Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 741,034 A→G 5.7% intergenic (+109/+41) 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,0340AG 5.7% 201.2 / 2.8 71intergenic (+109/+41)phr/dtpDdeoxyribodipyrimidine photolyase, FAD‑binding/dipeptide and tripeptide permease D
Reads supporting (aligned to +/- strand):  ref base A (34/32);  new base G (2/2);  total (36/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.33e-01

GCAGCGTATGAGGCGGCGCGGAAGGGGAAATAACAGACGCGTCAGGCAATCGAGCCCAGATGCCGGATGCGGCGTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACAGCCAAATCATCAGTACCAGTCCTACACACGCCAGTGCGCCCCAGGTAATTTGATCAAA  >  NC_000913/740893‑741182
                                                                                                                                             |                                                                                                                                                    
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:157189/1‑149 (MQ=255)
                                                   gAGCCCAGATGCCGGATGCGGCGTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTccagcgccag                                                                                            >  1:205633/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                                                                >  2:614665/1‑149 (MQ=255)
                                                                               gCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGa                                                                >  1:814737/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           >  1:314888/3‑149 (MQ=38)
                                                                                                                                    atgCCTGATGCGACGCTTGACGCGTCTTATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACAGCCAAATCATCAGTACCAGTCCTACACACGCCAGTGCGCCCCAGGTAAt           >  2:687385/3‑149 (MQ=38)
                                                                                                                                     ggCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACAGCCAATCATCAGTACCAGTCCTACACACGCCAGTGCGCCCCAGGTAAttt         <  2:455175/149‑1 (MQ=255)
                                                                                                                                        cTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACAGCCAAATCATCAGTACCAGTCCTACACACGCCAGTGCGCCCCAGGTAAtt          >  1:207856/1‑146 (MQ=255)
                                                                                                                                        cTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACAGCCAAATCATCAGTACCAGTCCTACACACGCCAGTGCGCCCCAGGTAAtt          <  2:207856/146‑1 (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)
                                                                                                                                             |                                                                                                                                                    
GCAGCGTATGAGGCGGCGCGGAAGGGGAAATAACAGACGCGTCAGGCAATCGAGCCCAGATGCCGGATGCGGCGTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACAGCCAAATCATCAGTACCAGTCCTACACACGCCAGTGCGCCCCAGGTAATTTGATCAAA  >  NC_000913/740893‑741182

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: