Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 741,041 G→T 5.5% intergenic (+116/+34) phr → / ← dtpD deoxyribodipyrimidine photolyase, FAD‑binding/dipeptide and tripeptide permease D

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913741,0410GT 5.5% 200.1 / 2.1 73intergenic (+116/+34)phr/dtpDdeoxyribodipyrimidine photolyase, FAD‑binding/dipeptide and tripeptide permease D
Reads supporting (aligned to +/- strand):  ref base G (35/34);  new base T (2/2);  total (37/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.42e-01

GCAGCGTATGAGGCGGCGCGGAAGGGGAAATAACAGACGCGTCAGGCAATCGAGCCCAGATGCCGGATGCGGCGTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACAGCCAAATCATCAGTACCAGTCCTACACACGCCAGTGCGCCCCAGGTAATTTGATCAAATACTTC  >  NC_000913/740893‑741188
                                                                                                                                                    |                                                                                                                                                   
gCAGCGTATGAGGCGGCGCGGAAGGGGAAATAACAGACGCGTCAGGCAATCGAGCCCAGATGCCGGATGCGGCGTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGAcgcg                                                                                                                                                     >  1:302262/1‑149 (MQ=255)
 cAGCGTATGAGGCGGCGCGGAAGGGGAAATAACAGACGCGTCAGGCAATCGAGCCCAGATGCCGGATGCGGCGTGAACGCCTTTTCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGt                                                                                                                                                    >  1:337822/1‑149 (MQ=255)
 cAGCGTATGAGGCGGCGCGGAAGGGGAAATAACAGACGCGTCAGGCAATCGAGCCCAGATGCCGGATGCGGCGTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGt                                                                                                                                                    >  1:455175/1‑149 (MQ=255)
 cAGCGTATGAGGCGGCGCGGAAGGGGAAATAACAGACGCGTCAGGCAATCGAGCCCAGATGCCGGATGCGGCGTGAACGCCTTATCCAGACTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGt                                                                                                                                                    >  1:385220/1‑149 (MQ=255)
  aGCGTATGAGGCGGCGCGGAAGGGGAAATAACAGACGCGTCAGGCAATCGAGCCCAGATGCCGGATGCGGCGTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTc                                                                                                                                                   <  1:1042020/149‑1 (MQ=255)
   gCGTATGAGGCGGCGCGGAAGGGGAAATAACAGACGCGTCAGGCAATCGAGCCCAGATGCCGGATGCGGCGTGCACGCCTTATCCAGCCTACAAATTTGGGCAAATTCAATAATTTGCAGGAAAGACGTAGGCCTGATAAGACGCGTCa                                                                                                                                                  >  1:242689/1‑149 (MQ=255)
              ggcgCGGAAGGGGAAATAACAGACGCGTCAGGCAATCGAGCCCAGATGCCGGATGCGGCGTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATc                                                                                                                                       >  1:51297/1‑149 (MQ=255)
               gcgcGGAAGGGGAAATAACAGACGCGTCAGGCAATCGAGCCCAGATGCCGGATGCGGCGTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCa                                                                                                                                      <  2:169341/149‑1 (MQ=255)
                cgcgGAAGGGGAAATAACAGACGCGTCAGGCAATCGAGCCCAGATGCCGGATGCGGCGTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAg                                                                                                                                     <  2:212843/149‑1 (MQ=255)
                cgcgGAAGGGGAAATAACAGACGCGGCAGGCAATCGAGCCGAGATGCCGGATGCGGCGTGAACGCCTTATCCAGCCTACAACTTTGTGCAAATTCAATAAATTGCAGGAAATACGTATGCCTGATGCGACGCTTGACGCGTCTTATCAg                                                                                                                                     <  2:484264/149‑1 (MQ=21)
                 gcgGAAGGGGAAATAACAGACGCGTCAGGCAATCGAGCCCAGATGCCGGATGCGGCGTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATGCGACGCTTGACGCGTCTTATCAgg                                                                                                                                    <  2:730154/149‑1 (MQ=25)
                    gAAGGGGAAATAACAGACGCGTCAGGCAATCGAGCCCAGATGCCGGATGCGGCGTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCAt                                                                                                                                 >  1:91555/1‑149 (MQ=255)
                       gaggAAATACCAGACGCGTCAGGCAATCGAGCCCAGATGCCGGAGGCGGCGTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCgg                                                                                                                              <  2:68162/147‑1 (MQ=255)
                          gAAATAACAGACGCGTCAGGCAATCGAGCCCAGATGCCGGATGCGGCGTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGc                                                                                                                           >  2:382985/1‑149 (MQ=255)
                             aTAACAGACGCGTCAGGCAATCGAGCCCAGATGCCGGATGCGGCGTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTgg                                                                                                                        >  1:398293/1‑149 (MQ=255)
                                 cAGACGCGTCAGGCAATCGAGCCCAGATGCCGGATGCGGCGTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCtt                                                                                                                    >  1:220095/1‑148 (MQ=255)
                                 cAGACGCGTCAGGCAATCGAGCCCAGATGCCGGATGCGGCGTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTa                                                                                                                    >  2:106378/1‑149 (MQ=255)
                                  aGACGCGTCAGGCAATCGAGCCCAGATGCCGGATGCGGCGTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCTCATCAGGCATCGCTGCTGGCCTAt                                                                                                                   >  1:1061590/1‑149 (MQ=255)
                                  aGACGCGTCAGGCAATCGAGCCCAGATGCCGGATGCGGCGTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTAt                                                                                                                   <  1:472714/149‑1 (MQ=255)
                                     cgcgTCAGGGAATCGAGCCACGATTCCGGATGCGGCGTGACAGCCTTATCCAGCCTACAAATTTGTGCAAATTCACTAAATTGCAGGAAATACTTAGGCATGATAAGACTCGTCATGCGTCGCATCAGGCATCGGTGCTGGCCTATTaa                                                                                                                <  1:773985/149‑1 (MQ=37)
                                      gcgTCAGGCAATCGAGCCCAGATGCCGGATGCGGCGTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAg                                                                                                               >  1:468341/1‑149 (MQ=255)
                                          cAGGCAATCGAGCCCAGATGCCGGATGCGGCGTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTc                                                                                                           <  2:156679/149‑1 (MQ=255)
                                            ggCAATCGAGCCCAGATGCCGGATGCGGCGTGAACGCCTTATCCAGCCTACAAATTTGTTCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCACGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTcca                                                                                                         >  2:106465/1‑149 (MQ=255)
                                            ggCAATCGAGCCCAGATGCCGGATGCGGCGTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTccc                                                                                                         >  2:674295/1‑148 (MQ=255)
                                                aTCGAGCCCAGATGCCGGATGCGACGTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTccagcgc                                                                                                     >  1:9263/1‑149 (MQ=255)
                                                   gaccccAGCTGCCGGACGCGGCGTGCACGCCTTATCCAGCCTACAAATTTGTGCAACTTCCATACATTGCCGGAAATACGTAGGCCTGCTCAGACGCGTCAAGCGTCGCATCAGGCAGCGGTGCGGGCCTATTAAGACTccagcgccag                                                                                                  <  2:504599/146‑1 (MQ=37)
                                                   gAGCCCAGATGCCGGATGCGGCGTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTccagcgccag                                                                                                  >  1:205633/1‑149 (MQ=255)
                                                   gAGCCCAGATGCCGGATGCGGCGTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTccagcgccag                                                                                                  >  1:157189/1‑149 (MQ=255)
                                                    aGCCCAGATGCCGGATGCGGCGTGAACGCCTTATCCAGCATACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGAAAAGCCGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTccagcgccagc                                                                                                 <  1:111299/149‑1 (MQ=255)
                                                         aGATGCCGGATGCGGCGTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGTTGCGGGCCTAGTAAGACTCCAGCGCCAGcgcgcg                                                                                            >  2:406749/1‑149 (MQ=255)
                                                         aGATGCCGGATGCGGCGTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGcgcgcg                                                                                            >  1:329349/1‑149 (MQ=255)
                                                             gCCGGATGCGGCGTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGttt                                                                                        <  2:601577/149‑1 (MQ=255)
                                                               cGGATGCGGCGTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCt                                                                                      <  1:706410/149‑1 (MQ=255)
                                                                 gATGCGGCGTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGTCGCGGCAAGCGTCGCATCAGGCATCGGTGCTGGCGTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGa                                                                                    <  2:367978/149‑1 (MQ=255)
                                                                  ataaggcgTTCACGCCGCATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGaa                                                                                   >  2:192180/5‑149 (MQ=38)
                                                                       gcgTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCa                                                                              >  2:971062/1‑149 (MQ=255)
                                                                        cgTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAg                                                                             <  1:971062/149‑1 (MQ=255)
                                                                         gTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAgc                                                                            >  1:1161981/1‑149 (MQ=255)
                                                                            aaCGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGCCTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGcc                                                                         >  1:108713/1‑149 (MQ=255)
                                                                            aaCGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGcc                                                                         <  2:33333/149‑1 (MQ=255)
                                                                               gcgttATGCAGCCTAGAAATTTGTGCTACATCAATAAATTGCAGGAAATTCGTTGGCATGATAAGTCGGGGCATGCGTCGCATCAGGCATCGGTGCTGGCCTATTTAGACTGGAGCGCCAGCGCGCGGTTTCTGAATTTGTGCGCCTGa                                                                      <  2:645990/146‑1 (MQ=255)
                                                                               gCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGa                                                                      >  1:814737/1‑149 (MQ=255)
                                                                               gCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGa                                                                      >  2:614665/1‑149 (MQ=255)
                                                                               gCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGa                                                                      <  1:743191/149‑1 (MQ=255)
                                                                                 cTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATa                                                                    <  1:106378/149‑1 (MQ=255)
                                                                                 cTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATa                                                                    <  1:697388/149‑1 (MQ=255)
                                                                                     tCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGGCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACAGc                                                                <  1:810556/149‑1 (MQ=255)
                                                                                      ccAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACAGcc                                                               >  2:107011/1‑149 (MQ=255)
                                                                                       cAGCCTACAAATTTGTGCAAATTCATTAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCGTATTAAGCCTCCAGCGCCAGCGCGAGGTTTCTGAGTTTCAGCGCCTGATACAGccc                                                              >  2:719476/1‑148 (MQ=39)
                                                                                        aGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACAGCCaa                                                             >  1:596673/1‑149 (MQ=255)
                                                                                             aCAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACAGCCAAatcat                                                        <  2:137598/149‑1 (MQ=255)
                                                                                             aCAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACAGCCAAatcat                                                        <  2:173836/149‑1 (MQ=255)
                                                                                              cAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGGGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACAGCCAAatcatc                                                       <  2:1115188/149‑1 (MQ=255)
                                                                                                   ttGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACAGCCAAATCATCAGTAc                                                  <  1:300762/149‑1 (MQ=255)
                                                                                                   ttGGGCAAATTCAATAAATAGCTGGAAATACGTAGGCCTGATAAGCCGCGTCAAGCGTCGCATCTGGCAGCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACAGCCAAATCATCAGTAc                                                  <  2:1054586/149‑1 (MQ=255)
                                                                                                         aaaTTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACAGCCAAATCATCAGTACCAGTcc                                            >  1:37372/1‑149 (MQ=255)
                                                                                                          caTTCAATACATGGCAGGAAATACGTAGGCCTGTTAAGACGCGTCAAGCGTAGAATCAGGCATCGGTGCTGGCGTCTTACGACGCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACCGCCAAATCATCAGTACCAGTCCt                                           <  2:786595/148‑1 (MQ=255)
                                                                                                                 tAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACAGCCAAATCATCAGTACCAGTCCTACACACg                                    <  1:883200/149‑1 (MQ=255)
                                                                                                                            aaaTACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACAGCCAAATCATCAGTACCAGTCCTACACACGCCAGTGCGccc                         <  2:302262/149‑1 (MQ=255)
                                                                                                                               tACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACAGCCAAATCATCAGTACCAGTCCTACACACGCCAGTGCGCCCCAg                      <  1:107011/149‑1 (MQ=255)
                                                                                                                                    atgCCTGATGCGACGCTTGACGCGTCTTATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACAGCCAAATCATCAGTACCAGTCCTACACACGCCAGTGCGCCCCAGGTAAt                 >  2:687385/3‑149 (MQ=38)
                                                                                                                                    atgCCTGATGCGACGCTTGACGCGTCTTATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACAGCCAAATCATCAGTACCAGTCCTACACACGCCAGTGCGCCCCAGGTAAt                 >  1:314888/3‑149 (MQ=38)
                                                                                                                                     ggCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACAGCCAATCATCAGTACCAGTCCTACACACGCCAGTGCGCCCCAGGTAAttt               <  2:455175/149‑1 (MQ=255)
                                                                                                                                        cTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACAGCCAAATCATCAGTACCAGTCCTACACACGCCAGTGCGCCCCAGGTAAtt                <  2:207856/146‑1 (MQ=255)
                                                                                                                                        cTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACAGCCAAATCATCAGTACCAGTCCTACACACGCCAGTGCGCCCCAGGTAAtt                >  1:207856/1‑146 (MQ=255)
                                                                                                                                        cTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACAGCCAAATCATCAGTACCAGTCCTACACACGCCAGTGCGCCCCAGGTAATTTGa             <  1:222337/149‑1 (MQ=255)
                                                                                                                                         tGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACAGCCAAATCATCAGTACCAGTCCTACACACGCCAGTGCGCCCCAGGTAATTTGAt            <  2:205633/149‑1 (MQ=255)
                                                                                                                                          gATAAGTCGCGTCAAGCGTCGCATCAGGCATGGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACAGCCAAATCATCAGTACCAGTCCTACACACGCCAGTGCGCCCCAGGTAATTTGATc           <  2:1161966/149‑1 (MQ=255)
                                                                                                                                          gATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACAGCCAAATCATCAGTACCAGTCCTACACACGCCAGTGCGCCCCAGGTAATTTGATc           >  1:448588/1‑149 (MQ=255)
                                                                                                                                          gATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACAGCCAAATCATCAGTACCAGTCATACACACGCCAGAGCGCCCCAGGTAATTTGATc           >  1:878806/1‑149 (MQ=255)
                                                                                                                                           aTAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACAGCCAAATCATCAGTACCAGTCCTACCCACGCCAGTGCGCCCCAGGTAATTTGATCa          >  1:914535/1‑149 (MQ=255)
                                                                                                                                             aaGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACAGCCAAATCATCAGTACCAGTCCTACACACGCCAGTGCGCCCCAGGTAATTTGATCaaa        <  1:614665/149‑1 (MQ=255)
                                                                                                                                                aCGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACAGCCAAATCATCAGTACCAGTCCTACACACGCCAGTGCGCCCCAGGTAATTTGATCAAATAc     >  2:611933/1‑149 (MQ=255)
                                                                                                                                                aCGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACAGCCAAATCATCAGTACCAGTCCTACACACGCCAGTGCGCCCCAGGGAATGTGAGCAAATAc     >  2:88445/1‑149 (MQ=255)
                                                                                                                                                  gcgTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACAGCCAAATCATCAGTACCAGTCCTACACACGCCAGTGCGCCCCAGGTAATTTGATCAAATACtt   <  1:621993/149‑1 (MQ=255)
                                                                                                                                                   cgTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACAGCCAAATCATCAGTACCAGTCCTACACACGCCAGTGCGCCCCAGGTAATTTGATCAAATACTTc  >  1:980197/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GCAGCGTATGAGGCGGCGCGGAAGGGGAAATAACAGACGCGTCAGGCAATCGAGCCCAGATGCCGGATGCGGCGTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACAGCCAAATCATCAGTACCAGTCCTACACACGCCAGTGCGCCCCAGGTAATTTGATCAAATACTTC  >  NC_000913/740893‑741188

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: