Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,229,438 Δ1 bp 7.0% intergenic (+76/+15) yjbB → / ← pepE putative Na+/Pi‑cotransporter/peptidase E, alpha‑aspartyl dipeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,229,4380C. 7.0% 457.2 / 27.3 143intergenic (+76/+15)yjbB/pepEputative Na+/Pi‑cotransporter/peptidase E, alpha‑aspartyl dipeptidase
Reads supporting (aligned to +/- strand):  ref base C (66/67);  new base . (4/6);  total (70/73)
Fisher's exact test for biased strand distribution p-value = 7.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.56e-01

Only 100 of 144 total aligned reads displayed.

ACTCGCTGTTTTGTTCGGTGGCTTACAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCG  >  NC_000913/4229289‑4229586
                                                                                                                                                      |                                                                                                                                                    
  tCGCTGTTTTGTTCGGTGGCTTACAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCAt                                                                                                                                                    <  2:1134971/150‑1 (MQ=255)
   cGCTGTTTTGTTCGGTGGCCTACAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATc                                                                                                                                                   >  2:500356/1‑150 (MQ=255)
     cTGTTTTGTACGGTGGCTTACAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCTTCGCATCAGGCACTATGAAAAAGATGCCGGATCAATCATCCg                                                                                                                                                 >  2:594676/1‑150 (MQ=39)
       gTTTTGTTCGGTGGCTTACAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATGCATGATCCGGc                                                                                                                                               <  1:1176351/150‑1 (MQ=255)
       gTTTTGTTCGGTGGCTTACAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGc                                                                                                                                               <  1:1411967/150‑1 (MQ=255)
       gTTTTGTTCGGTGGCTTACAGGGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGc                                                                                                                                               <  2:91779/150‑1 (MQ=255)
         tttGTTCGGTGGCTTACAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATGCATGATCCGGCat                                                                                                                                             <  1:1582794/150‑2 (MQ=255)
         tttGTTCGGTGGCTTACAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATGCATGATCCGGCat                                                                                                                                             <  2:1255211/150‑2 (MQ=255)
          ttGTTCGGTGGCTTACAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATGCATGATCCGGCatc                                                                                                                                            <  1:88382/150‑3 (MQ=255)
          ttGTTCGGTGGCTTACAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAAt                                                                                                                                            <  1:1398145/150‑1 (MQ=255)
           tGTTCGGTGGCTTACAGTGTTCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTACGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATGCATGATCCGGCatct                                                                                                                                           <  1:1320554/150‑4 (MQ=38)
           tGTTCGGTGGCTTACAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGAATCAGGCACAATGAACAAGATGCCGGATCATGCATCCGGCAAtt                                                                                                                                           >  1:759824/1‑150 (MQ=255)
            gTTCGGTGGCTTACAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTc                                                                                                                                          <  2:1250227/150‑1 (MQ=255)
              tCGTTGGCTTACAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATGCATGATCCGGCatcttgt                                                                                                                                        <  1:1649034/150‑7 (MQ=38)
              tCGGTGGCTTACAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCAc                                                                                                                                        <  1:529270/150‑1 (MQ=255)
               cGGTGGCTTTCAGTTTGCTGGAACAGCCGGATGAAGATGAAGGACGCGACGAGTATAAATATTGTGAATCCTGTAGGCCTGATACGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACt                                                                                                                                       <  1:1074585/150‑1 (MQ=39)
               cGGTGGCTTACAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACt                                                                                                                                       <  2:378505/150‑1 (MQ=255)
                ggTGGCTTACAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACtt                                                                                                                                      >  2:853723/1‑150 (MQ=255)
                 gTGGCTTACAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTa                                                                                                                                     >  1:277731/1‑150 (MQ=255)
                 gTGGCTTACAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTa                                                                                                                                     >  2:1633843/1‑150 (MQ=255)
                  tGGCTTACAGTGTTCTGGAACAGCCGGATGAATATGTAGGACGGGACGAGTATTAATTTTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTaa                                                                                                                                    <  2:150557/150‑1 (MQ=255)
                      ttACAGTGTGCTGGAACAGCCGGATGAAGATGAATGACGTGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAAGAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAAc                                                                                                                                <  1:249537/150‑1 (MQ=255)
                        aCAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACgg                                                                                                                              >  2:728136/1‑150 (MQ=255)
                          aGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTg                                                                                                                            <  1:1022796/150‑1 (MQ=255)
                            tgtgCTGGAACAGCCGGATGAAGATGAAGGACGGGACTAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGAc                                                                                                                          <  1:1108737/150‑1 (MQ=255)
                             gtgCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGAcc                                                                                                                         <  1:754815/150‑1 (MQ=255)
                               gcgggAAAAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCAGTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAg                                                                                                                       <  1:1510058/147‑1 (MQ=255)
                                cTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATAGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGc                                                                                                                      >  1:481312/1‑150 (MQ=255)
                                    aaCAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTcc                                                                                                                  <  1:1545593/150‑1 (MQ=255)
                                       aGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGc                                                                                                               <  1:1728300/150‑1 (MQ=255)
                                        gCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCg                                                                                                              >  1:1171469/1‑150 (MQ=255)
                                          cGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGa                                                                                                            <  1:988609/150‑1 (MQ=255)
                                           ggATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGaa                                                                                                           >  1:820212/1‑150 (MQ=255)
                                              tgaagctgaagGACGGGACGAGTATTAATCTTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGTATCATGCATCCGGCAATTCACTTAAAAACTGTGACCAGCTTCCAGCGGAACCg                                                                                                        <  2:216427/150‑1 (MQ=255)
                                                 agatgaagGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGctt                                                                                                     >  1:1035429/1‑150 (MQ=255)
                                                     gaagGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCttcttc                                                                                                 <  2:811509/150‑1 (MQ=255)
                                                      aagGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCTTGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCa                                                                                                <  1:695397/150‑1 (MQ=255)
                                                       agGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCAGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCAc                                                                                               >  2:921133/1‑150 (MQ=255)
                                                          acgggacgAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAg                                                                                            <  2:1212357/150‑1 (MQ=255)
                                                           cgggacgAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAGCGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGc                                                                                           >  1:1252041/1‑150 (MQ=255)
                                                           cgggacgAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGc                                                                                           >  2:1406877/1‑150 (MQ=255)
                                                             ggacgAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGTCACTATGAACAAGATGCCGGATCATGCATCCGTCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCTCCAGCCt                                                                                         >  2:1370987/1‑150 (MQ=255)
                                                             ggacgAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCt                                                                                         <  1:1692296/150‑1 (MQ=255)
                                                                         atatTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAg                                                                             >  1:86919/1‑150 (MQ=255)
                                                                           atTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAgtg                                                                           >  1:1690878/1‑150 (MQ=255)
                                                                             tgtgAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAgtggt                                                                         <  1:1594887/150‑1 (MQ=255)
                                                                                    ccTGTATGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCggg                                                                  <  1:460700/150‑1 (MQ=255)
                                                                                    ccTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCAGTTAAAAACGGTGACCAGCTTCCAGCGGCCCCGCTTCTTCCCCAGCCATAAACACATAAGTGGTGTTCggg                                                                  >  1:1091662/1‑150 (MQ=38)
                                                                                    ccTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCggg                                                                  <  1:1283009/150‑1 (MQ=255)
                                                                                     cTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGc                                                                 <  2:1276285/150‑1 (MQ=255)
                                                                                      tGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGcc                                                                <  1:360213/150‑1 (MQ=255)
                                                                                       gTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCa                                                               >  1:638964/1‑150 (MQ=255)
                                                                                         aGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCAcc                                                             <  1:1524162/150‑1 (MQ=255)
                                                                                           gCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGACCCGCTTCTTCACCAGCCTTAAACCCATAAGTGGTGTTCGGTCCACCCa                                                           >  2:1094204/1‑150 (MQ=255)
                                                                                            ccTGATGCGACGCTTGACGCGTCTTATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCAACAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCag                                                          >  1:877975/1‑150 (MQ=255)
                                                                                             cTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCagc                                                         >  2:140940/1‑150 (MQ=255)
                                                                                             cTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCagc                                                         >  1:917708/1‑150 (MQ=255)
                                                                                              tGATAAGACGCGTCAAGCGTCGCATCAGGCACTTTGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCagcg                                                        >  2:911665/1‑150 (MQ=255)
                                                                                              tGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCagcg                                                        <  1:1226684/150‑1 (MQ=255)
                                                                                               gATAATACGCGTCAAGAGTCGCATCAGGAACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCAATTAAAAACGGTGACCAGCTTCAATCGGAACCCCTAGATCACCAGCCTTAAACCCATAAGTGGTGTTCGGGCCCCCCcgcgt                                                       >  2:157090/1‑150 (MQ=255)
                                                                                                aTAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCagcgta                                                      >  2:544031/1‑150 (MQ=255)
                                                                                                  aaGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGTATCATGCATACGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCagcgtagc                                                    <  2:1675747/150‑1 (MQ=255)
                                                                                                  aaGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCagcgtagc                                                    >  2:768222/1‑150 (MQ=255)
                                                                                                        cgTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGAcc                                              <  1:154710/150‑1 (MQ=255)
                                                                                                          tCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCtt                                            >  1:695824/1‑150 (MQ=255)
                                                                                                          tCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCtt                                            >  1:382157/1‑150 (MQ=255)
                                                                                                          tCAAGCGTCGCATCAGGCACTATGAACAAGATGCAGGATCATGCATCCGGCAATTCACTTAAAACCTGTGACCAGCTTCCCGCGGAACCGATTCTACAACAGCCTTAAACACATAACAGGTTAACGGGCCACCCAGCGAAGCGTCTCCtt                                            >  1:75558/1‑150 (MQ=255)
                                                                                                              gCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACt                                        <  2:890704/150‑1 (MQ=255)
                                                                                                                 tCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCAc                                     <  2:1429208/150‑1 (MQ=255)
                                                                                                                   gCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACtg                                   >  1:547730/1‑150 (MQ=255)
                                                                                                                   gCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACtg                                   <  1:1159848/150‑1 (MQ=255)
                                                                                                                     aTCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACtgtg                                 >  1:1154226/1‑150 (MQ=255)
                                                                                                                        aGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATc                              <  2:1035429/150‑1 (MQ=255)
                                                                                                                        aGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATc                              <  2:438346/150‑1 (MQ=255)
                                                                                                                         ggCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATcc                             <  2:312567/150‑1 (MQ=255)
                                                                                                                         ggCACAATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATcc                             <  2:759824/150‑1 (MQ=255)
                                                                                                                           cACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAg                           <  2:583879/150‑1 (MQ=255)
                                                                                                                           cACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACAGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAg                           <  1:131809/150‑1 (MQ=255)
                                                                                                                            aCTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGt                          >  2:794349/1‑150 (MQ=255)
                                                                                                                             cTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGtt                         >  2:358307/1‑150 (MQ=255)
                                                                                                                              tATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATTAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTa                        >  2:193874/1‑150 (MQ=255)
                                                                                                                                tGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCCGCGGAACCGCTTCTTCCCCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCCCTGTGATCCAGTTAcc                      >  1:1522415/1‑150 (MQ=255)
                                                                                                                                tGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTAcc                      <  1:1504143/150‑1 (MQ=255)
                                                                                                                                 gAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCt                     >  1:183984/1‑150 (MQ=255)
                                                                                                                                 gAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCt                     <  1:887198/150‑1 (MQ=255)
                                                                                                                                  aaCAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTTGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCtt                    <  2:1356370/150‑1 (MQ=255)
                                                                                                                                   aCAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCAGCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTATGATCCAGTTACCTTc                   >  2:1114149/1‑150 (MQ=255)
                                                                                                                                     gaattGCCGGATGCATGATCCGGCAATTCACTTAAAACCGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGCCCTTTACTCACTGTGATCAAGTTACCTTCCg                 >  1:394490/5‑150 (MQ=255)
                                                                                                                                      aattGCCGGATGCATGATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCCGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCgg                >  2:1360320/4‑150 (MQ=255)
                                                                                                                                      aattGCCGGATGCATGATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGCCCTTTACTCACTGTGATCCAGTTACCTTCCgg                >  2:198561/4‑150 (MQ=255)
                                                                                                                                      aGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGATTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAACTGGTGTTCGGGCCACCCAGCGGAGCGTGCCCTGTTCTCCCTGTGATCCAGTTACCTTCCgg                >  2:697467/1‑150 (MQ=255)
                                                                                                                                       attGCCGGATGCATGATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGt               >  1:809265/3‑150 (MQ=255)
                                                                                                                                         tGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAg             >  1:1612055/1‑150 (MQ=255)
                                                                                                                                          gCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGa            >  2:143614/1‑150 (MQ=255)
                                                                                                                                           ccGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGAc           <  1:1633843/150‑1 (MQ=255)
                                                                                                                                            cGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCATAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGAcc          <  2:845671/150‑1 (MQ=255)
                                                                                                                                             ggATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCa         >  2:513662/1‑150 (MQ=255)
                                                                                                                                             ggATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCa         >  1:150277/1‑150 (MQ=255)
                                                                                                                                                     gCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCg  <  2:789971/150‑1 (MQ=255)
                                                                                                                                                     gCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCg  >  1:793513/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
ACTCGCTGTTTTGTTCGGTGGCTTACAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCG  >  NC_000913/4229289‑4229586

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