Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 580,020 G→A 100% R47R (AGG→AGA ybcW → DLP12 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913580,0200GA99.1% 12782.6 / 21.2 4405R47R (AGG→AGAybcWDLP12 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (26/19);  new base A (2113/2218);  total (2155/2250)
Fisher's exact test for biased strand distribution p-value = 2.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 4405 total aligned reads displayed.

GACAAAATGAATAAAGAACAATCTGCTGATGATCCCTCCGTGGATCTGATTCGTGTAAAAAATATGCTTAATAGCACCATTTCTATGAGTTACCCTGATGTTGTAATTGCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTA  >  NC_000913/579874‑580166
                                                                                                                                                  |                                                                                                                                                  
 aCAAAATGAATAAAGAACAATCTGCTGATGATCCCTCCGTGGATCTGATTCGTGTAAAAAATATGCTTAATAGCACCATTTCTATGAGTTACCCTGATGTTGTAATTGCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGAGc                                                                                                                                                  >  6:253538/1‑148 (MQ=255)
      aTGAATAAAGAACAATCTGCTGATGATCCCTCCGTGGATCTGATTCGTGTAAAAAATATGCTTAATAGCACCATTTCTATGAGTTACCCTGATGTTGTAATTGCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGa                                                                                                                                            <  2:41348/149‑1 (MQ=255)
       ttaaTAAACAACAATCTGCTCATGATCCCTCCGTTGATCTGATTCGTTTAAAAAATATGCTTAATAGCACCATTTCTATGAGTTACCCTGTTGTTGTAATTGCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGAg                                                                                                                                           <  6:112840/147‑1 (MQ=255)
          aTAAAGAACAATCTGCTGATGATCCCTCCGTGGATCTGATTCGTGTAAAAAATATGCTTAATAGCACCATTTCTATGAGTTACCCTGATGTTGTAATTGCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGAgg                                                                                                                                          >  1:17686/1‑147 (MQ=255)
          aTAAAGAACAATCTGCTGATGATCCCTCCGTGGATCTGATTCGTGTAAAAAATATGCTTAATAGCACCACTTCTATGAGTTACCCTGATGTTGTAATTGCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGAgg                                                                                                                                          >  7:328311/1‑147 (MQ=255)
           tAAAGAACAATCTGCTGATGATCCCTCCGTGGATCTGATTCGTGTAAAAAATATGCTTAATAGCACCATTTCTATGAGTTACCCTGATGTTGTAATTGCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGAGGCAg                                                                                                                                       >  3:266421/1‑149 (MQ=255)
            aaaGAACAATCTGCTGATGATCCCTCCGTGGATCTGATTCGTGTAAAAAATATGCTTAATAGCACCATTTCTATGAGTTACCCTGATGTTGTAATTGCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGAGGCAGc                                                                                                                                      >  7:111281/1‑149 (MQ=255)
                aaCAATCTGCTGATGATCCCTCCGTGGATCTGATTCGTGTAAAAAATATGCTTAATAGCACCATTTCTATGAGTTACCCTGATGTTGTAATTGCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGAGGCAGCGTTg                                                                                                                                  <  2:259939/149‑1 (MQ=255)
                   aaTCTGCTGATGATCCCTCCGTGGATCTGATTCGTGTAAAAAATATGCTTAATAGCACCATTTCTATGAGTTACCCTGATGTTGTAATTGCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGAgg                                                                                                                                          <  5:37717/138‑1 (MQ=255)
                       tgctgATGATCCCTCCGTGGATCTGATTCGTGTAAAAAATATGCTTAATAGCCCCATTTCTATGAGTTACCCTGATGTTGTAATTGCATGTATAGAACATATGGTGTCTCTGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAg                                                                                                                            >  3:8155/1‑148 (MQ=255)
                         ctgATGATCCCTCCGTGGATCTGATTCGTGTAAAAAATATGCTTAATAGCACCATTTCTATGAGTTACCCTGATGTTGTAATTGCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAg                                                                                                                            <  2:313800/146‑1 (MQ=255)
                         ctgATGATCCCTCCGTGGATCTGATTCGTGTAAAAAATATGCTTAATAGCACCATTTCTATGAGTTACCCTGATGTTGTAATTGCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCa                                                                                                                          <  1:47095/148‑1 (MQ=255)
                         ctgATGATCCCTCCGTGGATCTGATTCGTGTAAAAAATATGCTTAATAGCACCATTTCTATGAGTTACCCTGATGTTGTAATTGCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCAc                                                                                                                         <  1:25637/149‑1 (MQ=255)
                         ctgATGATCCCTCCGTGGATCTGATTCGTGTAAAAAATATGCTTAATAGCACCATTTCTATGAGTTACCCTGATGTTGTAATTGCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCAc                                                                                                                         <  1:312894/149‑1 (MQ=255)
                          tgatgaTCCCTCCGTGGATCTGATTCGTGTAAAAAATATGCTTAATAGCACCATTTCTATGAGTTACCCTGATGTTGTAATTGCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGa                                                                                                                                            <  8:318788/129‑1 (MQ=255)
                           gatgatCCCTCCGTGGATCTGATTCGTGTAAAAAATATGCTTAATAGCACCATTTCTATGAGTTACCCTGATGTTGTAATTGCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGAGGCAGCGTTg                                                                                                                                  <  3:307418/138‑1 (MQ=255)
                           gatgatCCCTCCGTGGATCTGATTCGTGTAAAAAATATGCTTAATAGCACCATTTCTATGAGTTACCCTGATGTTGTAATTGCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGa                                                                                                                       >  5:212564/1‑149 (MQ=255)
                            atgatCCCTCCGTGGATCTGATTCGTGTAAAAAATATGCTTAATAGCACCATTTCTATGAGTTACCCTGATGTTGTAATTGCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGAGCAATTg                                                                                                                                             <  8:183853/126‑1 (MQ=255)
                            atgatCCCTCCGTGGATCTGATTCGTGTAAAAAATATGCTTAATAGCACCATTTCTATGAGTTACCCTGATGTTGTAATTGCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGa                                                                                                                       >  5:319853/1‑148 (MQ=255)
                               atCCCTCCGTGGATCTGATTCGTGTAAAAAATATGCTTAATAGCACCATTTCTATGAGTTACCCTGATGTTGTAATTGCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGataa                                                                                                                    <  4:84121/148‑1 (MQ=255)
                                tCCCTCCGTGGATCTGATTCGTGTAAAAAATATGCTTAATAGCACCATTTCTATGAGTTACCCTGATGTTGTAATTGCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGataat                                                                                                                   <  8:160247/148‑1 (MQ=255)
                                 cccTCCGTGGATCTGATTCGTGTAAAAAATATGCTTAATAGCACCATTTCTATGAGTTACCCTGATGTTGTAATTGCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGataataa                                                                                                                 >  6:137759/1‑149 (MQ=255)
                                 cccTCCGTGGATCTGATTCGTGTAAAAAATATGCTTAATAGCACCATTTCTATGAGTTACCCTGATGTTGTAATTGCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGAACGataataa                                                                                                                 >  2:139626/1‑149 (MQ=255)
                                   ctacGTTGATCTGATTCGTGTAAAAAATATGCTTAATAGCACCATTTCTATGAGTTACCCTGATGTTGTAATTGCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGataataata                                                                                                               <  4:359911/146‑1 (MQ=255)
                                   cTCCGTGGATCTGATTCGTGTAAAAAATATGCTTAATAGCACCATTTCTATGAGTTACCCTGATGTTGTAATTGCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGAGGCAGCGTTg                                                                                                                                  <  5:155111/130‑1 (MQ=255)
                                   cTCCGTGGATCTGATTCGTGTAAAAAATATGCTTAATAGCACCATTTCTATGAGTTACCCTGATGTTGTAATTGCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGataataata                                                                                                               >  4:214944/1‑149 (MQ=255)
                                   cTCCGTGGATCTGATTCGTGTAAAAAATATGCTTAATAGCACCATTTCTATGAGTTACCCTGATGTTGTAATTGCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGataataata                                                                                                               <  8:45950/149‑1 (MQ=255)
                                    tCCGTGGATCTGATTCGTGTAAAAAATATGCTTAATAGCACCATTTCTATGAGTTACCCTGATGTTGTAATTGCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAatat                                                                                                              >  7:269868/1‑149 (MQ=255)
                                      cGTGGATCTGATTCGTGTAAAAAATATGCTTAATAGCACCATTTCTATGAGTTACCCTGATGTTGTAATTGCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAatat                                                                                                              <  4:141053/147‑1 (MQ=255)
                                          gATCTGATTCGTGTAAAAAATATGCTTAATAGCACCATTTCTATGAGTTACCCTGATGTTGTAATTGCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGa                                                                                                        <  7:44259/149‑1 (MQ=255)
                                                 ttCGTGTAAAAAATATGCTTAATAGCACCATTTCTATGAGTTACCCTGATGTTGTAATTGCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGa                                                                                                        <  1:371028/142‑1 (MQ=255)
                                                  tCGTGTAAAAAATATGCTTAATAGCACCATTTCTATGAGTTACCCTGATGTTGTAATTGCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTccc                                                                                                >  8:349912/1‑149 (MQ=255)
                                                    gtgtAAAAAATTTGCTTAATAGCACCATTTCTATGAGTTACCCTGATGTTGTAATTGCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCt                                                                                               <  5:10805/148‑1 (MQ=255)
                                                    gtgtAAAAAATATGCTTAATAGCACCATTTCTATGAGTTACCCTGATGTTGTAATTGCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCtg                                                                                              >  1:23241/1‑149 (MQ=255)
                                                    gtgtAAAAAATATGCTTAATAGCACCATTTCTATGAGTTACCCTGATGTTGTAATTGCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCtg                                                                                              <  4:166847/149‑1 (MQ=255)
                                                    gggtAAAAAATATGCTTAATAGCACCATTTCTATGAGTTACCCTGATGTTGTAATTGCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCt                                                                                               <  5:268420/146‑1 (MQ=255)
                                                     tgtAAAAAATATGCTTAATAGCACCATTTCTATGAGTTACCCTGATGTTGTAATTGCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCtgg                                                                                             <  3:140664/149‑1 (MQ=255)
                                                     tgtAAAAAATATGCTTAATAGCACCATTTCTATGAGTTACCCTGATGTTGTAATTGCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCtgg                                                                                             >  7:142817/1‑149 (MQ=255)
                                                      gtAAAAAATATGCTTAATAGCACCATTTCTATGAGTTACACGGATGTTGTAATTGCATGTATAGAACATAAGGTGTCACTGGAAGCAATCAGAGCAATTGAGGCAGCGTGGGTGAAGCACGATAATAATATGAAGGATTATTCCCtggt                                                                                            >  6:279225/1‑149 (MQ=255)
                                                       tAAAAAATATGCTTAATAGCACCATTTCTATGAGTTACCCTGATGTTGTAATTGCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTccc                                                                                                <  2:23715/144‑1 (MQ=255)
                                                       tAAAAAATATGCTTAATAGCACCATTTCTATGAGTTACCCTGATGTTGTAATTGCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTcc                                                                                                 >  6:369763/1‑143 (MQ=255)
                                                          aaaaTATGCTTAATAGCACCATTTCTAGGAGTTACCCTGATGTTGTAATTGCATGGATAGAACATAAGGTGTCTCTGGAAGCATTCAGAGCAAATGAGGCAGCGTTGGTGAAGCACGATAATAATATGAACGATTAATCCCTGGTGGtt                                                                                        >  5:347077/1‑149 (MQ=255)
                                                             atatGCTTAATAGCACCATTTCTATGAGTTACCCTGATGTTGTAATTGCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGAc                                                                                     <  3:167523/149‑1 (MQ=255)
                                                                tGCTTAATAGCACCATTTCTATGAGTTACCCTGATGTTGTAATTGCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGa                                                                                  <  6:307537/149‑1 (MQ=255)
                                                                   ttAATAGCACCATTTCTATGAGTTACCCTGATGTTGTAATTGCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGataataa                                                                                                                 <  5:134781/115‑1 (MQ=255)
                                                                   ttAATAGCACCATTTCTATGAGTTACCCTGATGTTGTAATTGCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCa                                                                               >  2:78074/1‑149 (MQ=255)
                                                                      aTAGCACCATTTCTATGAGTTACCCTGATGTTGTAATTGCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCAcc                                                                             <  7:8038/148‑1 (MQ=255)
                                                                      aTAGCACCATTTCTATGAGTTACCCTGATGTTGTAATTGCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCAcc                                                                             <  6:86792/148‑1 (MQ=255)
                                                                      aTAGCACCATTTCTATGAGTTACCCTGATGTTGTAATTGCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCa                                                                            <  8:29166/149‑1 (MQ=255)
                                                                       tAGCACCATTTCTATGAGTTACCCTGATGTTGTAATTGCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGataa                                                                                                                    <  1:201947/108‑1 (MQ=255)
                                                                         gCACCATTTCTATGAGTTACCCTGATGTTGTAATTGCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATaa                                                                         >  8:382293/1‑149 (MQ=255)
                                                                          cACCATTTCTATGAGTTACCCTGATGTTGTAATTGCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAAc                                                                        >  4:272258/1‑149 (MQ=255)
                                                                          cACCATTTCTATGAGTTACCCTGATGTTGTAATTGCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAAc                                                                        <  5:245140/149‑1 (MQ=255)
                                                                           aCCATTTCTATGAGTTACCCTTATGTTGTAATTGCATTTATAGAACATAATGTGTCTCTGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAAc                                                                        <  1:124486/148‑1 (MQ=255)
                                                                               tttCTATGAGTTACCCTGATGTTGTAATTGCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCt                                                                    >  1:81103/1‑148 (MQ=255)
                                                                                   tATGAGTTACCCTGATGTTGTAATTGCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGAt                                                                                 <  6:335534/131‑1 (MQ=255)
                                                                                   tATGAGTTACCCTGATGTTGTAATTGCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCa                                                               >  3:104611/1‑149 (MQ=255)
                                                                                     tGAGTTACCCTGATGTTGTAATTGCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGAGGCGGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCa                                                               <  8:263912/147‑1 (MQ=255)
                                                                                       aGTTACCCTGATGTTGTAATTGCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCa                                                           <  7:349193/149‑1 (MQ=255)
                                                                                       aGTTACCCTGATGTTGTAATTGCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCa                                                           <  8:105843/149‑1 (MQ=255)
                                                                                           aCCCTGATGTTGTAATTGCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAac                                                        <  7:24425/148‑1 (MQ=255)
                                                                                             ccTGATGTTGTAATTGCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAactat                                                     >  2:198030/1‑148 (MQ=255)
                                                                                              cTGATGTTGTAATTGCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAactatt                                                    >  8:238360/1‑147 (MQ=255)
                                                                                               tGATGTTGTAATTGCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCAtt                                                             >  3:193984/1‑139 (MQ=255)
                                                                                               tGATGTTGTAATTGCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAactatt                                                    >  6:161654/1‑146 (MQ=255)
                                                                                               tGATGTTGTAATTGCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCACAACTGCTAATCATTCAAactatt                                                    >  3:361895/1‑146 (MQ=255)
                                                                                                 aTGTTGTAATTGCATGTATAGAACATAAGGTGTCTCTGGAATCATTCAGAGCAATTGAGGCAGCGTTGGTGTAGCACTATAATAATATGAAGGATTATTCCCTGGTGTTTGACTGATCACCATAACTTCTAATCATTCAAACTATTTa                                                  <  4:108404/148‑1 (MQ=255)
                                                                                                 aTGTTGTAATTGCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTa                                                  <  8:70026/148‑1 (MQ=255)
                                                                                                    ttgtAATTGCATGTATATAACATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTagtct                                              <  2:343557/149‑5 (MQ=255)
                                                                                                      gtAATTGCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTagtctg                                             <  4:19058/148‑6 (MQ=255)
                                                                                                       tAATTGCATGTATAGACCATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGAGGCCGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGGCTTATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGT                                            <  6:189080‑M1/148‑1 (MQ=255)
                                                                                                       tAATTGCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGa                                                                                                                                            >  2:284292/1‑52 (MQ=255)
                                                                                                       tAATTGCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGT                                            <  2:109316‑M1/148‑1 (MQ=255)
                                                                                                          ttGCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCtggt                                                                                            >  8:208661/1‑97 (MQ=255)
                                                                                                          ttGCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTa                                                  <  3:333515/139‑1 (MQ=255)
                                                                                                            gCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAG                                       <  5:9063‑M1/148‑1 (MQ=255)
                                                                                                              aTGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCtggt                                                                                            <  4:58434/93‑1 (MQ=255)
                                                                                                              aTGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGc                                    >  7:300434/1‑149 (MQ=255)
                                                                                                              aTGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGc                                    <  8:138524/149‑1 (MQ=255)
                                                                                                               tGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGCGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGc                                    >  3:218223/1‑148 (MQ=255)
                                                                                                               tGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGcc                                   <  5:299197/149‑1 (MQ=255)
                                                                                                                 tataGAACATAAGGTGGCTCTGGAAGCTTTCAGAGCAATTGAGGCAGCGTTGTGGAAGCACGATAATATTATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCaa                                 <  6:56735/149‑1 (MQ=255)
                                                                                                                   taGAACATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGAGGCAGCGTTGCTGAAGAACGAGAATAATATGAAGGATTATTCCCTGGTGGTTGTCTGATCCCCATAACTGCTAAGCATTAAAACTATTTAGTCTGTGACAGAGCCAacc                               >  4:246024/1‑148 (MQ=255)
                                                                                                                     gAACATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACg                             <  4:332260/149‑1 (MQ=255)
                                                                                                                     gAACATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACg                             >  4:308361/1‑149 (MQ=255)
                                                                                                                      aaCATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCATACTATTTAGTCTGTGACAGAGCCAACACg                             >  1:286204/1‑148 (MQ=255)
                                                                                                                      aaCATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACg                             >  7:169817/1‑148 (MQ=255)
                                                                                                                      aaCATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACg                             >  3:5852/1‑148 (MQ=255)
                                                                                                                            aGGTGTCTCTGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGtct                       >  2:171953‑M2/1‑145 (MQ=255)
                                                                                                                             ggTGTCTCTGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTACCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGtctgt                     <  5:185653‑M2/149‑6 (MQ=255)
                                                                                                                                 tctcTGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTa                                                  >  6:224792/1‑116 (MQ=255)
                                                                                                                                   tctGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGTctgtcactgt               >  7:377712/1‑149 (MQ=255)
                                                                                                                                   tctGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGTctgtcactgt               >  2:366415/1‑149 (MQ=255)
                                                                                                                                   tctGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGTctgtcactgt               >  7:218823/1‑149 (MQ=255)
                                                                                                                                   tctGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGTctgtcactgt               >  1:150242/1‑149 (MQ=255)
                                                                                                                                   tctGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGTctgtcactgt               >  7:193041/1‑149 (MQ=255)
                                                                                                                                      ggAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAg            >  1:33972/1‑149 (MQ=255)
                                                                                                                                       gAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGTctgtcactgtca             <  5:61166/147‑1 (MQ=255)
                                                                                                                                            aTTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCAGTGACAGAGCCAACCCGCAGTCTGTCACTGTCAGGAAAGt      >  1:101526/1‑149 (MQ=255)
                                                                                                                                             ttCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTg     >  3:244855/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GACAAAATGAATAAAGAACAATCTGCTGATGATCCCTCCGTGGATCTGATTCGTGTAAAAAATATGCTTAATAGCACCATTTCTATGAGTTACCCTGATGTTGTAATTGCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTA  >  NC_000913/579874‑580166

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