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.

ACTCGCTGTTTTGTTCGGTGGCTTACAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCG  >  NC_000913/4229289‑4229586
                                                                                                                                                      |                                                                                                                                                    
  tCGCTGTTTTGTTCGGTGGCTTACAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCAt                                                                                                                                                    <  2:1134971/150‑1 (MQ=255)
   cGCTGTTTTGTTCGGTGGCCTACAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATc                                                                                                                                                   >  2:500356/1‑150 (MQ=255)
     cTGTTTTGTACGGTGGCTTACAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCTTCGCATCAGGCACTATGAAAAAGATGCCGGAT‑CAATCATCCg                                                                                                                                                 >  2:594676/1‑150 (MQ=39)
       gTTTTGTTCGGTGGCTTACAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATGCATG‑ATCCGGc                                                                                                                                               <  1:1176351/150‑1 (MQ=255)
       gTTTTGTTCGGTGGCTTACAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGc                                                                                                                                               <  1:1411967/150‑1 (MQ=255)
       gTTTTGTTCGGTGGCTTACAGGGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGc                                                                                                                                               <  2:91779/150‑1 (MQ=255)
         tttGTTCGGTGGCTTACAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATGCATG‑ATCCGGCat                                                                                                                                             <  1:1582794/150‑2 (MQ=255)
         tttGTTCGGTGGCTTACAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATGCATG‑ATCCGGCat                                                                                                                                             <  2:1255211/150‑2 (MQ=255)
          ttGTTCGGTGGCTTACAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATGCATG‑ATCCGGCatc                                                                                                                                            <  1:88382/150‑3 (MQ=255)
          ttGTTCGGTGGCTTACAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAAt                                                                                                                                            <  1:1398145/150‑1 (MQ=255)
           tGTTCGGTGGCTTACAGTGTTCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTACGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATGCATG‑ATCCGGCatct                                                                                                                                           <  1:1320554/150‑4 (MQ=38)
           tGTTCGGTGGCTTACAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGAATCAGGCACAATGAACAAGATGCCGGAT‑CATGCATCCGGCAAtt                                                                                                                                           >  1:759824/1‑150 (MQ=255)
            gTTCGGTGGCTTACAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTc                                                                                                                                          <  2:1250227/150‑1 (MQ=255)
              tCGTTGGCTTACAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATGCATG‑ATCCGGCatcttgt                                                                                                                                        <  1:1649034/150‑7 (MQ=38)
              tCGGTGGCTTACAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCAc                                                                                                                                        <  1:529270/150‑1 (MQ=255)
               cGGTGGCTTTCAGTTTGCTGGAACAGCCGGATGAAGATGAAGGACGCGACGAGTATAAATATTGTGAATCCTGTAGGCCTGATACGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACt                                                                                                                                       <  1:1074585/150‑1 (MQ=39)
               cGGTGGCTTACAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACt                                                                                                                                       <  2:378505/150‑1 (MQ=255)
                ggTGGCTTACAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACtt                                                                                                                                      >  2:853723/1‑150 (MQ=255)
                 gTGGCTTACAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTa                                                                                                                                     >  1:277731/1‑150 (MQ=255)
                 gTGGCTTACAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTa                                                                                                                                     >  2:1633843/1‑150 (MQ=255)
                  tGGCTTACAGTGTTCTGGAACAGCCGGATGAATATGTAGGACGGGACGAGTATTAATTTTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTaa                                                                                                                                    <  2:150557/150‑1 (MQ=255)
                      ttACAGTGTGCTGGAACAGCCGGATGAAGATGAATGACGTGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAAGAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAAc                                                                                                                                <  1:249537/150‑1 (MQ=255)
                        aCAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACgg                                                                                                                              >  2:728136/1‑150 (MQ=255)
                          aGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTg                                                                                                                            <  1:1022796/150‑1 (MQ=255)
                            tgtgCTGGAACAGCCGGATGAAGATGAAGGACGGGACTAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGAc                                                                                                                          <  1:1108737/150‑1 (MQ=255)
                             gtgCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGAcc                                                                                                                         <  1:754815/150‑1 (MQ=255)
                               gcgggAAAAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCAGTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAg                                                                                                                       <  1:1510058/147‑1 (MQ=255)
                                cTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATAGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGc                                                                                                                      >  1:481312/1‑150 (MQ=255)
                                    aaCAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTcc                                                                                                                  <  1:1545593/150‑1 (MQ=255)
                                       aGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGc                                                                                                               <  1:1728300/150‑1 (MQ=255)
                                        gCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCg                                                                                                              >  1:1171469/1‑150 (MQ=255)
                                          cGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGa                                                                                                            <  1:988609/150‑1 (MQ=255)
                                           ggATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGaa                                                                                                           >  1:820212/1‑150 (MQ=255)
                                              tgaagctgaagGACGGGACGAGTATTAATCTTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGTAT‑CATGCATCCGGCAATTCACTTAAAAACTGTGACCAGCTTCCAGCGGAACCg                                                                                                        <  2:216427/150‑1 (MQ=255)
                                                 agatgaagGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGctt                                                                                                     >  1:1035429/1‑150 (MQ=255)
                                                     gaagGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCttcttc                                                                                                 <  2:811509/150‑1 (MQ=255)
                                                      aagGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CTTGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCa                                                                                                <  1:695397/150‑1 (MQ=255)
                                                       agGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCAGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCAc                                                                                               >  2:921133/1‑150 (MQ=255)
                                                          acgggacgAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAg                                                                                            <  2:1212357/150‑1 (MQ=255)
                                                           cgggacgAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAGCGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGc                                                                                           >  1:1252041/1‑150 (MQ=255)
                                                           cgggacgAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGc                                                                                           >  2:1406877/1‑150 (MQ=255)
                                                             ggacgAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGTCACTATGAACAAGATGCCGGAT‑CATGCATCCGTCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCTCCAGCCt                                                                                         >  2:1370987/1‑150 (MQ=255)
                                                             ggacgAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCt                                                                                         <  1:1692296/150‑1 (MQ=255)
                                                                         atatTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAg                                                                             >  1:86919/1‑150 (MQ=255)
                                                                           atTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAgtg                                                                           >  1:1690878/1‑150 (MQ=255)
                                                                             tgtgAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAgtggt                                                                         <  1:1594887/150‑1 (MQ=255)
                                                                                    ccTGTATGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCggg                                                                  <  1:460700/150‑1 (MQ=255)
                                                                                    ccTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCAGTTAAAAACGGTGACCAGCTTCCAGCGGCCCCGCTTCTTCCCCAGCCATAAACACATAAGTGGTGTTCggg                                                                  >  1:1091662/1‑150 (MQ=38)
                                                                                    ccTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCggg                                                                  <  1:1283009/150‑1 (MQ=255)
                                                                                     cTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGc                                                                 <  2:1276285/150‑1 (MQ=255)
                                                                                      tGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGcc                                                                <  1:360213/150‑1 (MQ=255)
                                                                                       gTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCa                                                               >  1:638964/1‑150 (MQ=255)
                                                                                         aGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCAcc                                                             <  1:1524162/150‑1 (MQ=255)
                                                                                           gCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGACCCGCTTCTTCACCAGCCTTAAACCCATAAGTGGTGTTCGGTCCACCCa                                                           >  2:1094204/1‑150 (MQ=255)
                                                                                            ccTGATGCGACGCTTGACGCGTCTTATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCAACAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCag                                                          >  1:877975/1‑150 (MQ=255)
                                                                                             cTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCagc                                                         >  2:140940/1‑150 (MQ=255)
                                                                                             cTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCagc                                                         >  1:917708/1‑150 (MQ=255)
                                                                                              tGATAAGACGCGTCAAGCGTCGCATCAGGCACTTTGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCagcg                                                        >  2:911665/1‑150 (MQ=255)
                                                                                              tGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCagcg                                                        <  1:1226684/150‑1 (MQ=255)
                                                                                               gATAATACGCGTCAAGAGTCGCATCAGGAACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCAATTAAAAACGGTGACCAGCTTCAATCGGAACCCCTAGATCACCAGCCTTAAACCCATAAGTGGTGTTCGGGCCCCCCcgcgt                                                       >  2:157090/1‑150 (MQ=255)
                                                                                                aTAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCagcgta                                                      >  2:544031/1‑150 (MQ=255)
                                                                                                  aaGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGTAT‑CATGCATACGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCagcgtagc                                                    <  2:1675747/150‑1 (MQ=255)
                                                                                                  aaGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCagcgtagc                                                    >  2:768222/1‑150 (MQ=255)
                                                                                                        cgTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGAcc                                              <  1:154710/150‑1 (MQ=255)
                                                                                                          tCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCtt                                            >  1:695824/1‑150 (MQ=255)
                                                                                                          tCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCtt                                            >  1:382157/1‑150 (MQ=255)
                                                                                                          tCAAGCGTCGCATCAGGCACTATGAACAAGATGCAGGAT‑CATGCATCCGGCAATTCACTTAAAACCTGTGACCAGCTTCCCGCGGAACCGATTCTACAACAGCCTTAAACACATAACAGGTTAACGGGCCACCCAGCGAAGCGTCTCCtt                                            >  1:75558/1‑150 (MQ=255)
                                                                                                              gCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACt                                        <  2:890704/150‑1 (MQ=255)
                                                                                                                 tCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCAc                                     <  2:1429208/150‑1 (MQ=255)
                                                                                                                   gCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACtg                                   >  1:547730/1‑150 (MQ=255)
                                                                                                                   gCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACtg                                   <  1:1159848/150‑1 (MQ=255)
                                                                                                                     aTCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACtgtg                                 >  1:1154226/1‑150 (MQ=255)
                                                                                                                        aGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATc                              <  2:1035429/150‑1 (MQ=255)
                                                                                                                        aGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATc                              <  2:438346/150‑1 (MQ=255)
                                                                                                                         ggCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATcc                             <  2:312567/150‑1 (MQ=255)
                                                                                                                         ggCACAATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATcc                             <  2:759824/150‑1 (MQ=255)
                                                                                                                           cACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAg                           <  2:583879/150‑1 (MQ=255)
                                                                                                                           cACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACAGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAg                           <  1:131809/150‑1 (MQ=255)
                                                                                                                            aCTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGt                          >  2:794349/1‑150 (MQ=255)
                                                                                                                             cTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGtt                         >  2:358307/1‑150 (MQ=255)
                                                                                                                              tATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATTAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTa                        >  2:193874/1‑150 (MQ=255)
                                                                                                                                tGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCCGCGGAACCGCTTCTTCCCCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCCCTGTGATCCAGTTAcc                      >  1:1522415/1‑150 (MQ=255)
                                                                                                                                tGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTAcc                      <  1:1504143/150‑1 (MQ=255)
                                                                                                                                 gAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCt                     >  1:183984/1‑150 (MQ=255)
                                                                                                                                 gAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCt                     <  1:887198/150‑1 (MQ=255)
                                                                                                                                  aaCAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTTGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCtt                    <  2:1356370/150‑1 (MQ=255)
                                                                                                                                   aCAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCAGCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTATGATCCAGTTACCTTc                   >  2:1114149/1‑150 (MQ=255)
                                                                                                                                     gaattGCCGGATGCATG‑ATCCGGCAATTCACTTAAAACCGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGCCCTTTACTCACTGTGATCAAGTTACCTTCCg                 >  1:394490/5‑150 (MQ=255)
                                                                                                                                      aattGCCGGATGCATG‑ATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCCGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCgg                >  2:1360320/4‑150 (MQ=255)
                                                                                                                                      aattGCCGGATGCATG‑ATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGCCCTTTACTCACTGTGATCCAGTTACCTTCCgg                >  2:198561/4‑150 (MQ=255)
                                                                                                                                      aGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGATTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAACTGGTGTTCGGGCCACCCAGCGGAGCGTGCCCTGTTCTCCCTGTGATCCAGTTACCTTCCgg                >  2:697467/1‑150 (MQ=255)
                                                                                                                                       attGCCGGATGCATG‑ATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGt               >  1:809265/3‑150 (MQ=255)
                                                                                                                                         tGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAg             >  1:1612055/1‑150 (MQ=255)
                                                                                                                                          gCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGa            >  2:143614/1‑150 (MQ=255)
                                                                                                                                           ccGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGAc           <  1:1633843/150‑1 (MQ=255)
                                                                                                                                            cGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCATAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGAcc          <  2:845671/150‑1 (MQ=255)
                                                                                                                                             ggAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCa         >  2:513662/1‑150 (MQ=255)
                                                                                                                                             ggAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCa         >  1:150277/1‑150 (MQ=255)
                                                                                                                                                     gCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCg  <  2:789971/150‑1 (MQ=255)
                                                                                                                                                     gCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCg  >  1:793513/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
ACTCGCTGTTTTGTTCGGTGGCTTACAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCG  >  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: ‑