Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 572,454 C→T C29C (TGC→TGT ylcH → uncharacterized protein, DLP12 prophage

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913572,4540CT99.0% 320.9 / ‑4.0 98coding (87/102 nt)ylcHuncharacterized protein, DLP12 prophage
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (49/48);  minor base . (1/0);  total (50/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.03e-01

CCGGAGCCTGAACAACTCCGGTGACTTCTGCGCTAAACGGGGACGTTTATGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTGGAGGAGGGCAAATGCTTTCGCATGGTGCTTAAACCGTGGCGT  >  NC_000913/572320‑572585
                                                                                                                                      |                                                                                                                                   
ccGGAGCCTGAACAACTCCGGTGACTTCTGCGCTAAACGGGACGTTTATGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAgg                                                                                                                                >  1:811880/1‑139 (MQ=255)
 cGGAGCCTGAACAACTCCGGTGACTTCTGCGCTAAACGGGGACGTTTATGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAgg                                                                                                                                >  2:1719173/1‑139 (MQ=255)
  ggAGCCTGAACAACTCCGGTGACTTCTGCGCTAAACGGGGACGTTTATGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAggg                                                                                                                               <  2:703665/139‑1 (MQ=255)
  ggAGCCTGAACAACTCCGGTGACTTCTGCGCTAAACGGGGACGTTTATGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAggg                                                                                                                               <  2:294360/139‑1 (MQ=255)
     gCCTGAACAACTCCGGTGACTTCTGCGCTAAACGGGGACGTTTATGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGaa                                                                                                                            >  1:1342297/1‑139 (MQ=255)
            caacTCCGGTGACTTCTGCGCTAAACGGGGACGTTTATGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCATTTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGaa                                                                                                                     <  1:1334936/139‑1 (MQ=255)
              acTCCGGTGACTTCTGCGCTAAACGGGGACGTTTATGCGCACATACAATCCTAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGaa                                                                                                                     <  2:863454/137‑1 (MQ=255)
              acTCCGGTGACTTCTGCGCTAAACGGGGACGTTTATGCGCACATACAATCCTAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGaa                                                                                                                     >  1:863454/1‑137 (MQ=255)
              acTCCGGTGACTTCTGCGCTAAACGGGGACGTTTATGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAAcc                                                                                                                   <  1:1203774/139‑1 (MQ=255)
              acTCCGGTGACTTCTGCGCTAAACGGGGACGTTTATGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAAcc                                                                                                                   <  1:902234/139‑1 (MQ=255)
              acTCCGGTGACTTCTGCGCTAAACGGGGACGTTTATGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAAcc                                                                                                                   >  2:232388/1‑139 (MQ=255)
                    gTGACTTCTGCGCTAAACGGGGACGTTTATGCGCACATACAATCCAAACTCTCTTCTCCCTTCACGGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTca                                                                                                              <  2:1058466/138‑1 (MQ=255)
                    gTGACTTCTGCGCTAAACGGGGACGTTTATGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTcac                                                                                                             >  1:1815244/1‑139 (MQ=255)
                    gTGACTTCTGCGCTAAACGGGGACGTTTATGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTcac                                                                                                             >  2:1476353/1‑139 (MQ=255)
                    gTGACTTCTGCGCTAAACGGGGACGTTTATGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTcac                                                                                                             >  1:436350/1‑139 (MQ=255)
                    gTGACTTCTGCGCTAAACGGGGACGTTTATGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTca                                                                                                              >  1:1058466/1‑138 (MQ=255)
                         ttCTGCGCTAAACGGGGACGTTTATGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGAc                                                                                                        >  2:200199/1‑139 (MQ=255)
                            tGCGCTAAACGGGGACGTTTATGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAggg                                                                                                                               <  2:1234704/113‑1 (MQ=255)
                            tGCGCTAAACGGGGACGTTTATGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAggg                                                                                                                               >  1:1234704/1‑113 (MQ=255)
                            tGCGCTAAACGGGGACGTTTATGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGc                                                                                                     >  2:1626637/1‑139 (MQ=255)
                            tGCGCTAAACGGGGACGTTTATGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGc                                                                                                     >  2:22493/1‑139 (MQ=255)
                             gcgcTAAACGGGGACGTTTATGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCa                                                                                                    <  2:487859/139‑1 (MQ=255)
                             gcgcTAAACGGGGACGTTTATGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCa                                                                                                    <  1:1719173/139‑1 (MQ=255)
                             gcgcTAAACGGGACGTTTATGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCAt                                                                                                   <  2:811880/139‑1 (MQ=255)
                                        ggACGTTTATGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACAt                                                                                         <  2:241533/139‑1 (MQ=255)
                                              ttATGCGCACATTCAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGc                                                                                   <  2:1559068/139‑1 (MQ=255)
                                                      acaTACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACctct                                                                                                                <  2:1265541/102‑1 (MQ=255)
                                                      acaTACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACctct                                                                                                                >  1:1265541/1‑102 (MQ=255)
                                                        aTACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCg                                                                         >  2:173514/1‑139 (MQ=255)
                                                             aTCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAAcc                                                                                                                   >  1:995423/1‑92 (MQ=255)
                                                             aTCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAAcc                                                                                                                   <  2:995423/92‑1 (MQ=255)
                                                              tCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACtt                                                                              <  2:1224629/128‑1 (MQ=255)
                                                              tCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACtt                                                                              >  1:1224629/1‑128 (MQ=255)
                                                              tCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTAtcg                                                                   >  1:907636/1‑139 (MQ=255)
                                                              tCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTAtcg                                                                   >  2:1605900/1‑139 (MQ=255)
                                                               ccAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTAtcgg                                                                  <  1:173514/139‑1 (MQ=255)
                                                               ccAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTAtcgg                                                                  <  1:22493/139‑1 (MQ=255)
                                                                  aaCTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTAtcggtcg                                                               >  2:1215136/1‑139 (MQ=255)
                                                                   aCTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTAtcggtcgg                                                              <  2:1188303/139‑1 (MQ=255)
                                                                      ctctTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCa                                                                                                    >  2:1669523/1‑98 (MQ=255)
                                                                      ctctTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCa                                                                                                    <  1:1669523/98‑1 (MQ=255)
                                                                        cttctCCCTTCACAGATGCAGAAATGCACCTGTAATTCTTTGCATCTAGCGTTTGACCTCT‑‑‑GTGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCa                                                      >  1:1116845/1‑139 (MQ=255)
                                                                         ttctCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATc                                                        <  1:383409/139‑1 (MQ=255)
                                                                              ccTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGcc                                                   <  1:350951/139‑1 (MQ=255)
                                                                                 tCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGaa                                                                                                                     >  1:280645/1‑70 (MQ=255)
                                                                                 tCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGaa                                                                                                                     <  2:280645/70‑1 (MQ=255)
                                                                                 tCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCa                                                                                                    >  1:1111749/1‑87 (MQ=255)
                                                                                 tCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCa                                                                                                    <  2:1111749/87‑1 (MQ=255)
                                                                                   acaGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTc                                                                             >  1:1677902/1‑108 (MQ=255)
                                                                                   acaGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTc                                                                             <  2:1677902/108‑1 (MQ=255)
                                                                                     aaaTGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGc                                                             <  2:1143658/120‑1 (MQ=255)
                                                                                     aaaTGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGc                                                             >  1:1143658/3‑122 (MQ=255)
                                                                                      gATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTg                                                                                                                       >  1:1394328/1‑63 (MQ=255)
                                                                                      gATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTg                                                                                                                       <  2:1394328/63‑1 (MQ=255)
                                                                                      gATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGaa                                                                                                                     >  1:1371325/1‑65 (MQ=255)
                                                                                      gATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGaa                                                                                                                     <  2:1371325/65‑1 (MQ=255)
                                                                                      gATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATcctgtctcttatac                                           >  2:895988/1‑127 (MQ=255)
                                                                                        tGCAGAAATGCACCTGCAATTCTTTGCATCTAGTGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACgg                                                                                                      <  1:1334220/78‑1 (MQ=255)
                                                                                        tGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGa                                                                                                         >  1:850305/1‑75 (MQ=255)
                                                                                        tGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGa                                                                                                         <  2:850305/75‑1 (MQ=255)
                                                                                        tGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACgg                                                                                                      >  2:1334220/1‑78 (MQ=255)
                                                                                        tGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGcc                                                   <  1:1249359/129‑1 (MQ=255)
                                                                                        tGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGcc                                                   >  2:1249359/1‑129 (MQ=255)
                                                                                         gCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTggagg                                        <  2:436350/139‑1 (MQ=255)
                                                                                         gCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTggagg                                        <  1:1215136/139‑1 (MQ=255)
                                                                                          cAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGaa                                                                                                                     <  1:1571407/61‑1 (MQ=255)
                                                                                          cAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGaa                                                                                                                     >  2:1571407/1‑61 (MQ=255)
                                                                                          cAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTggagga                                       >  1:1596251/1‑139 (MQ=255)
                                                                                            gAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTGGaggagg                                     >  1:261502/1‑139 (MQ=255)
                                                                                                 gCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTGGAGGAGGGCaaa                                >  2:1490861/1‑139 (MQ=255)
                                                                                                  cACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTGGAGGAGGGCAAAt                               >  2:1550213/1‑139 (MQ=255)
                                                                                                   aCCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACctc                                                                                                                 >  2:182512/1‑56 (MQ=255)
                                                                                                   aCCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACctc                                                                                                                 <  1:182512/56‑1 (MQ=255)
                                                                                                    ccTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATc                                              >  1:1508992/1‑122 (MQ=255)
                                                                                                    ccTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATc                                              <  2:1508992/122‑1 (MQ=255)
                                                                                                       gCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTGGAGGAGGGCAAATGCttt                          >  2:414405/1‑139 (MQ=255)
                                                                                                       gCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTGGAGGAGGGCAAATGCtt                           >  1:1298397/1‑138 (MQ=255)
                                                                                                       gCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTGGAGGAGGGCAAATGCtt                           <  2:1298397/138‑1 (MQ=255)
                                                                                                            tCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTGGAGGAGGGCAAATGCTTTCGCAt                     <  2:261502/139‑1 (MQ=255)
                                                                                                            tCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTGGAGGAGGGCAAATGCTTTCGCAt                     <  1:414405/139‑1 (MQ=255)
                                                                                                             cTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTGGAGGAGGGCAAATGCTTTCGCATg                    <  1:1626637/139‑1 (MQ=255)
                                                                                                             cTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTGGAGGAGGGCAAATGCTTTCGCATg                    <  1:1490861/139‑1 (MQ=255)
                                                                                                               ttGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGa                                                                                                         <  1:1379904/52‑1 (MQ=255)
                                                                                                               ttGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGa                                                                                                         >  2:1379904/1‑52 (MQ=255)
                                                                                                                tGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAAc                                                                                                                    <  2:1226328/40‑1 (MQ=255)
                                                                                                                tGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAAc                                                                                                                    >  1:1226328/1‑40 (MQ=255)
                                                                                                                tGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTGGAGGAGGGCAAATGCTTTCGCATGGTg                 >  2:1362401/1‑139 (MQ=255)
                                                                                                                  cATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTGGAGGAGGGCAAATGCTTTCGCATGGTGCt               >  2:1389835/1‑139 (MQ=255)
                                                                                                                       aGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACgg                                                                                                      >  2:1370105/1‑47 (MQ=255)
                                                                                                                       aGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACgg                                                                                                      <  1:1370105/47‑1 (MQ=255)
                                                                                                                       aGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGc                                                             >  2:1597860/1‑88 (MQ=255)
                                                                                                                       aGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGc                                                             <  1:1597860/88‑1 (MQ=255)
                                                                                                                            ttGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTGGAGGAGGGCAAATGCTTTCGCATGGTGCTTAAACCGTgg     >  2:1795605/1‑139 (MQ=255)
                                                                                                                             tGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGc                                                             <  2:1340612/82‑1 (MQ=255)
                                                                                                                             tGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGc                                                             >  1:1340612/1‑82 (MQ=255)
                                                                                                                              gACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTAt                                                                     <  2:1589141/73‑1 (MQ=255)
                                                                                                                              gACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTAt                                                                     >  1:1589141/1‑73 (MQ=255)
                                                                                                                               aCCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTGGAGGAGGGCAAATGCTTTCGCATGGTGCTTAAACCGTGGCGt  <  2:685328/139‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                   
CCGGAGCCTGAACAACTCCGGTGACTTCTGCGCTAAACGGGGACGTTTATGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTGGAGGAGGGCAAATGCTTTCGCATGGTGCTTAAACCGTGGCGT  >  NC_000913/572320‑572585

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