Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,454,266 C→A 7.9% pseudogene (736/2646 nt) yfcU ← pseudogene, FimD fimbrial export usher family;putative membrane; Not classified; putative outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,454,2660CA 7.9% 1905.5 / 19.3 634pseudogene (736/2646 nt)yfcUpseudogene, FimD fimbrial export usher family;putative membrane; Not classified; putative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base C (287/272);  new base A (2/73);  total (289/345)
Fisher's exact test for biased strand distribution p-value = 9.86e-19
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 654 total aligned reads displayed.

TGGCGCGTAGCCGCGCAGATTGGGAGGCAACATTTGATCGTCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAGCGCGCATACGCCATGGCCCCAGGTTAACCCCGACCGTCCCGTTGCCGC  >  NC_000913/2454117‑2454415
                                                                                                                                                     |                                                                                                                                                     
tGGCGCGTAGCCGCGCAGATTGGGAGGCAACATTTGATCGTCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGc                                                                                                                                                       >  1:5971257/1‑150 (MQ=255)
tGGCGCGTAGCCGCGCAGATTGGGAGGCAACATTTGATCGTCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGc                                                                                                                                                       <  2:2915069/150‑1 (MQ=255)
tGGCGCGTAGCCGCGCAGATTGGGAGGCAACATTTGATCGTCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGc                                                                                                                                                       >  1:8558150/1‑150 (MQ=255)
gggCGCGTAGCCGCGCAGATTGGGAGGCAACATTTGATCGTCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCTGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGc                                                                                                                                                       <  1:7045358/149‑1 (MQ=255)
  gcgcGTAGCCGCGCAGATTGGGAGGCAACATTTGATCGTCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCat                                                                                                                                                     >  2:4789042/1‑150 (MQ=255)
   cgcgTAGCCGCGCAGATTGGGAGGCAACATTTGATCGTCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCata                                                                                                                                                    >  2:7658628/1‑150 (MQ=255)
    gcgTAGCCGCGCAGATTGGGAGGCAACATTTGATCGTCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCataa                                                                                                                                                   >  1:3429723/1‑150 (MQ=255)
    gcgTAGCCGCGCAGATTGGGAGGCAACATTTGATCGTCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCataa                                                                                                                                                   >  2:5497887/1‑150 (MQ=255)
    gcgTAGCCGCGCAGATTGGGAGGCAACATTTGATCGTCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCataa                                                                                                                                                   <  2:8794162/150‑1 (MQ=255)
    gcgTAGCCGCGCAGATTGGGAGGCAACATTTGATCGTCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCataa                                                                                                                                                   >  2:1334417/1‑150 (MQ=255)
     cgTAGCCGCGCAGATTGGGAGGCAACATTTGATCGTCAGTGCTGACACTGCCACCAACATAATTAACACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCataat                                                                                                                                                  >  1:3582633/1‑150 (MQ=255)
      gTAGCCGCGCAGATTGGGAGGCAACATTTGATCGTCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCataata                                                                                                                                                 <  2:6103135/150‑1 (MQ=255)
       tAGCCGCGCAGATTGGGAGGCAACATTTGATCGTCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAg                                                                                                                                                >  2:7207827/1‑150 (MQ=255)
           cgcgCAGATTGGGAGGCAACATTTGATCGTCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGac                                                                                                                                            <  2:2109840/150‑1 (MQ=255)
           cgcgCAGATTGGGAGGCAACATTTGATCGTCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGac                                                                                                                                            <  2:5689955/150‑1 (MQ=255)
            gcgcAGATTGGGAGGCAACATTTGATCGTCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGact                                                                                                                                           >  1:7486937/1‑150 (MQ=255)
             cgcAGTTTGGGAGGCAACATTTGATCGTCTGTGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCGGATTTGAGGTAATCTTCGCCCAGCGCCTGTTGGGCTTTTAATGAGGGTAACGCCCGCCCGGCATAATCGCGactc                                                                                                                                          <  2:4352955/150‑1 (MQ=255)
                  aTTGGGAGGCAACATTTGATCGTCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGactccactc                                                                                                                                     >  1:5578871/1‑150 (MQ=255)
                      ggAGGCAACATTTGATCGTCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCAt                                                                                                                                 >  2:5963127/1‑150 (MQ=255)
                       gAGGCAACATTTGATCGTCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCAtt                                                                                                                                >  1:1265039/1‑150 (MQ=255)
                           cAACATTTGATCGTCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCAtttttt                                                                                                                            >  1:8402891/1‑150 (MQ=255)
                            aaCATTTGATCGTCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCAttttttt                                                                                                                           >  1:3717021/1‑150 (MQ=255)
                             aCATTTGATCGTCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCAtttttttt                                                                                                                          >  1:2267257/1‑150 (MQ=255)
                              cATTTGATCGTCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTtg                                                                                                                         >  1:572740/1‑150 (MQ=255)
                                tttGATCGTCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCAtttttttgtgt                                                                                                                       >  2:3067048/1‑146 (MQ=255)
                                tttGATCGTCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTtgtg                                                                                                                       >  2:149390/1‑150 (MQ=255)
                                 tttaTCGTCTGTGCTGACACTGCCACCAACATAATTAAAACCTTCTAAAATATCGTAATTGAGTTAATCTTCGCCAAGCGCAATTTTGGATTTTAATGACGGTAACGCCCGCCTGGAATAATAGCGACTCCACTCCCATTTTTTTtgtgt                                                                                                                      <  2:3584813/147‑1 (MQ=255)
                                    aTCGTCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTc                                                                                                                   >  2:2536520/1‑150 (MQ=255)
                                    aTCGTCAGTGCTGACACTACCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTc                                                                                                                   >  2:2318379/1‑150 (MQ=255)
                                     tCGTCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGTGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCa                                                                                                                  <  2:5254572/150‑1 (MQ=255)
                                       gTCAGTGCTTACACTGCCACAAACATAATTAAAACCGTCTAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATc                                                                                                                <  2:9133343/150‑1 (MQ=255)
                                       gTCAGTGCTGACACTGCCACCAACATATTTAAAACGTTCTAAAATATCGGAATGGTGTTATTCTGCGCCCCGCGCCAGTTGGGATTTTATTGCCGTTAACGCCCGCCTGGCTTAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATc                                                                                                                <  2:3470364/150‑1 (MQ=255)
                                          tgTGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCGTATTTGAGGTAATCTTCGCCCCGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATcgcc                                                                                                             <  2:6302529/149‑1 (MQ=255)
                                            tGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATcgccgc                                                                                                           >  1:2749672/1‑150 (MQ=255)
                                               tGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCaa                                                                                                        >  1:8964914/1‑150 (MQ=255)
                                                 aaacccccacaaACAAAATAAAAACCTACCAAAATAACCAAATTCACCAAATCTACCCCCACCCCAAGTTTGTATTTTAATCACGCTAACCCCCCCCAGCAATAAAAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAAt                                                                                                      <  1:5204102/144‑1 (MQ=255)
                                                  cacTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAAtt                                                                                                     >  2:8406969/1‑150 (MQ=255)
                                                   acTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTc                                                                                                    >  2:6673817/1‑150 (MQ=255)
                                                     cgcaaccaACATAATTAAAACCGTCTAAAATATGGGATTTGTGTTAATCTTCGCCCCGCGCCTGTTTGGTTTTTATTGAGGTTACCGCCCGCCCGGAATAATCGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCAt                                                                                                  <  1:9201003/146‑1 (MQ=255)
                                                            aaCATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGtcatc                                                                                           >  1:1712523/1‑150 (MQ=255)
                                                            aaCATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGtcatc                                                                                           >  2:7254548/1‑150 (MQ=255)
                                                            aaCATAATTAAAACCGTCGAAAATAGCGTAATTGAGGTAATCTGCGCCCAGCGCCAGTTTGGCTTTTAATGACGTTAACGCCCGCCAGGAATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGtcatc                                                                                           <  2:4419503/150‑1 (MQ=255)
                                                              cATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGtcatcat                                                                                         >  1:4769217/1‑150 (MQ=255)
                                                                   ttAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTa                                                                                    >  1:7242330/1‑150 (MQ=255)
                                                                      aaaCCTTCTAATATATCTTATTTGAGTTATTCTTCGCCCAGCGCCTTTTTTTATTTTAATTAGGTTACCGCCCGCCTGGCAAAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTg                                                                                 <  1:149246/150‑1 (MQ=255)
                                                                       aaCCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTgc                                                                                >  1:7595117/1‑150 (MQ=255)
                                                                       aaCCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTgc                                                                                >  2:1963389/1‑150 (MQ=255)
                                                                        aCCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTgcg                                                                               <  2:2153803/150‑1 (MQ=255)
                                                                         ccGTCTAAAATATCGGATTTGAGGTAATCTCCGCCCCGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTgcgc                                                                              <  1:4610448/150‑1 (MQ=255)
                                                                         ccGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCTGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTgcgc                                                                              <  2:7975916/150‑1 (MQ=255)
                                                                         ccGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTgcgc                                                                              >  1:8048294/1‑150 (MQ=255)
                                                                             cGAAAATATCGTAATTGCGTTAATCTTCGCCCAGCGCAATTTTGGCTTTTAATGACGGTACCGCCCGCCAGGAATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTAt                                                                          <  1:2459218/150‑1 (MQ=255)
                                                                              gAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATg                                                                         >  2:6519132/1‑150 (MQ=255)
                                                                               aaaaTATCGGAATTGAGGTAATCTTCGCCCTGCGCCTTTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGt                                                                        <  1:6202922/150‑1 (MQ=255)
                                                                               aaaaTATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGt                                                                        <  2:6866497/150‑1 (MQ=255)
                                                                                aaaTATGGCATTTGTGTTAATCTGCGCCCAGCGCCTGTTTGTATTTTAATGACGTTAACGCCCGCCAGGAATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGtt                                                                       <  1:67345/150‑1 (MQ=255)
                                                                                aaaTATCGTAATTGAGGTAATCTTCGCCCAGCGCCAGTTGGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATCGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGtt                                                                       <  2:181163/150‑1 (MQ=255)
                                                                                  atatCGTATTTGAGGTAATCTTCGCCCAGCGCCTGTTTGGATTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGa                                                                     <  1:7155780/150‑1 (MQ=255)
                                                                                       gtaaTTGTGGAAATCTCCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGtt                                                                <  1:267053/148‑1 (MQ=255)
                                                                                       gtaaTTGAGGTAATCTGCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGAATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGtt                                                                <  1:2110229/148‑1 (MQ=255)
                                                                                       ggAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGtt                                                                <  1:1333736/150‑1 (MQ=255)
                                                                                       gcatttGAGGTAATCTTCGCCCAGCGCCTGTTTGGATTTTAATGACGTTACCGCCCGCCTGGAATAATCGCGACTCCACTCCCATTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTg                                                               <  2:3910207/146‑1 (MQ=255)
                                                                                          aTTGAGTTAATCTTCGCCCAGCGCCAGTTGGGCTTTTAATGACGGTAACGCCCGCCAGGAATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGt                                                             <  1:7415003/150‑1 (MQ=255)
                                                                                          aTTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGt                                                             <  1:5116368/150‑1 (MQ=255)
                                                                                          aTTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGAGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGt                                                             >  2:7021870/1‑150 (MQ=255)
                                                                                            tgtggTTATCTTCGCCCCGCGCCTGTTTGGATTTTAATGAGGGTAACGCCCGCCCGGCATAACCGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCt                                                           <  1:3136124/147‑1 (MQ=255)
                                                                                                gTAATCTTCGCCCAGCGCCTGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCa                                                       <  2:8267452/150‑1 (MQ=255)
                                                                                                 tAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAg                                                      >  1:2002438/1‑150 (MQ=255)
                                                                                                    tCTTCGCCCCGCGCCTTTTTGGCTTTTAATGACGGTAACGCCCGCCAGGAATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCa                                                   <  2:4154548/150‑1 (MQ=255)
                                                                                                    tCTTCGCCCAGCGCCTTTTTGGCTTTTAATGACGTTAACGCCCGCCCGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCa                                                   <  1:1082882/150‑1 (MQ=255)
                                                                                                       tCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAgcg                                                <  2:2750316/150‑1 (MQ=255)
                                                                                                        cGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAgcgc                                               <  1:6363199/150‑1 (MQ=255)
                                                                                                         gCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAgcgcg                                              <  1:1014731/150‑1 (MQ=255)
                                                                                                           ccAGCGCCTGTTTGGCTTTTACTGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAGCGCGCa                                            <  1:4957152/150‑1 (MQ=255)
                                                                                                           ccAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATATCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCAACATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAGCGCGCa                                            >  1:8671503/1‑150 (MQ=255)
                                                                                                           ccAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAGCGCGCa                                            >  1:4914909/1‑150 (MQ=255)
                                                                                                                gcCTGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAGCGCGCATACGc                                       <  1:4283650/150‑1 (MQ=255)
                                                                                                                 ccAGTTTGGCTTTTAATGACGGTAACGCCCACCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCAACATTACTGCGCGTATGTTGATAGTTGGTCTGCCCGTCAGCGCGCAAACGcc                                      >  1:4465897/1‑150 (MQ=255)
                                                                                                                  cattttGGCTTTTATTGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAGCGCGCATACGCCa                                     <  2:111368/147‑1 (MQ=255)
                                                                                                                  cAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAGCGCGCATACGCCa                                     <  1:8185989/150‑1 (MQ=255)
                                                                                                                      ttGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAGCGCGCATACGCCATGGc                                 <  1:6519132/150‑1 (MQ=255)
                                                                                                                      tgggCTTTTAATGACGGTACCGCCCGCCCGGCATAAGAGCGACTCTACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAGCGCGCATACGCCATGGc                                 <  2:5367434/148‑1 (MQ=255)
                                                                                                                           ttttAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAGCGCGCATACGCCATGGCCCCAg                            <  1:539777/150‑1 (MQ=255)
                                                                                                                           ttttAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAGCGCGCATACGCCATGGCCCCAg                            >  2:7778352/1‑150 (MQ=255)
                                                                                                                           ttttAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAGCGCGCATACGCCATGGCCCCAg                            >  1:8793834/1‑150 (MQ=255)
                                                                                                                                aTGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGtcatcat                                                                                         >  1:2943265/1‑84 (MQ=255)
                                                                                                                                aTGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAGCGCGCATACGCCATGGCCCCAGGTTaa                       <  2:7372745/150‑1 (MQ=255)
                                                                                                                                 tGAGGGTAACGCCCGCCAGGAATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAGCGCGCATACGCCATGGCCCCAGGTTAAc                      <  2:3347654/150‑1 (MQ=255)
                                                                                                                                   aCGGTAACGCCCGCCCGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAGCGCGCATACGCCATGGCCCCAGGTTAAccc                    <  2:8653327/150‑1 (MQ=255)
                                                                                                                                   aCGGTAACGCCCGCCAGGAATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAGCGCGCATACGCCATGGCCCCAGGTTAAccc                    <  2:1932779/150‑1 (MQ=255)
                                                                                                                                      gTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAGCGCGCATACGCCATGGCCCCAGGTTAACCCCGa                 <  2:7521840/150‑1 (MQ=255)
                                                                                                                                      gTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAGCGCGCATACGCCATGGCCCCAGGTTAACCCCGa                 >  1:1919870/1‑150 (MQ=255)
                                                                                                                                       tAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAGCGCGCATACGCCATGGCCCCAGGTTAACCCCGAc                >  2:1658420/1‑150 (MQ=255)
                                                                                                                                           gcccgccAGGCATAATCGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAGCGCGCATACGCCATGGCCCCAGGTTAACCCCGACCGTc            <  2:3755012/150‑1 (MQ=255)
                                                                                                                                           gcccgccAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAGCGCGCATACGCCATGGCCCCAGGTTAACCCCGACCGTc            <  1:4235914/150‑1 (MQ=255)
                                                                                                                                           gcccgccAGGAATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAGCGCGCATACGCCATGGCCCCAGGTTAACCCCGACCGTc            <  1:7768669/150‑1 (MQ=255)
                                                                                                                                             ccgccAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAGCGCGCATACGCCATGGCCCCAGGTTAACCCCGACCGTccc          <  1:3864676/150‑1 (MQ=255)
                                                                                                                                             ccgccAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAGCGCGCATACGCCATGGCCCCAGGTTAACCCCGACCGTccc          >  1:3643611/1‑150 (MQ=255)
                                                                                                                                             ccgccAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAGCGCGCATACGCCATGGCCCCAGGTTAACCCCGACCGTccc          <  1:3297255/150‑1 (MQ=255)
                                                                                                                                                  aGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAGCGCGCATACGCCATGGCCCCAGGTTAACCCCGACCGTCCCGTTGc     <  2:5134115/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGGCGCGTAGCCGCGCAGATTGGGAGGCAACATTTGATCGTCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAGCGCGCATACGCCATGGCCCCAGGTTAACCCCGACCGTCCCGTTGCCGC  >  NC_000913/2454117‑2454415

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: