Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 741,045 A→C 6.4% intergenic (+120/+30) 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,0450AC 6.4% 191.6 / 2.2 72intergenic (+120/+30)phr/dtpDdeoxyribodipyrimidine photolyase, FAD‑binding/dipeptide and tripeptide permease D
Reads supporting (aligned to +/- strand):  ref base A (31/34);  new base C (3/2);  total (34/38)
Fisher's exact test for biased strand distribution p-value = 6.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

GGCGCGGAAGGGGAAATAACAGACGCGTCAGGCAATCGAGCCCAGATGCCGGATGCGGCGTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACAGCCAAATCATCAGTACCAGTCCTACACACGCCAGTGCGCCCCAGGTAATTTGATCAAATACTTCAAT  >  NC_000913/740907‑741191
                                                                                                                                          |                                                                                                                                                  
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                                                                         <  1:743191/149‑1 (MQ=255)
                                                                 gCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGa                                                                         >  1:814737/1‑149 (MQ=255)
                                                                 gCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGa                                                                         >  2:614665/1‑149 (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                   >  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)
                                                                                                                                  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)
                                                                                                                                       tCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACAGCCAAATCATCAGTACCAGTCCTACACACGCCAGTGCGCCCCAGGTAATTTGATCAAATACTTCaa   <  2:385220/149‑1 (MQ=255)
                                                                                                                                        cAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACAGCCAAATCATCAGTACCAGTCCTACACACGCCAGTGCGCCCCAGGTAATTTGATCAAATACTTCAat  <  1:916157/149‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                                  
GGCGCGGAAGGGGAAATAACAGACGCGTCAGGCAATCGAGCCCAGATGCCGGATGCGGCGTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACAGCCAAATCATCAGTACCAGTCCTACACACGCCAGTGCGCCCCAGGTAATTTGATCAAATACTTCAAT  >  NC_000913/740907‑741191

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: