Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,900,679 T→C 8.3% I303V (ATC→GTC)  yoaE ← putative membrane protein/conserved protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,900,6790TC 8.3% 197.1 / 8.0 73I303V (ATC→GTC) yoaEputative membrane protein/conserved protein
Reads supporting (aligned to +/- strand):  ref base T (28/37);  new base C (3/3);  total (32/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.11e-01

CTTTAGCACGTACAATACCGATGATTTCATCCAGTTCACCGCGACATACCGGGAACAGACTGTGCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAAC  >  NC_000913/1900532‑1900825
                                                                                                                                                   |                                                                                                                                                  
cTTTAGCACGTACAATACCGATGATTTCATCCAGTTCACCGCGACATACCGGGAACAGACTGTGCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATcc                                                                                                                                                  >  2:123109/1‑150 (MQ=255)
    aGCACGTACAATACCGATGATTTCATCCAGTTCACCGCGACATACCGGGAACAGACTGTGCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCcgcg                                                                                                                                              <  1:331670/150‑1 (MQ=255)
     gCACGTACAATACCGATGATTTCATCCAGTTCACCGCGACATACCGGGAACAGACTGTGCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGTCGTCATGATCCCgcgc                                                                                                                                             >  1:95345/1‑150 (MQ=255)
     gCACGTACAATACCGATGATTTCATCCAGTTCACCGCGACATACCGGGAACAGACTGTGCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCgcgc                                                                                                                                             <  1:395341/150‑1 (MQ=255)
        cGTACAATACCGATGATTTCATCCAGTTCACCGCGACATACCGGGAACAGACTGTGCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCaga                                                                                                                                          <  1:445535/150‑1 (MQ=255)
         gTACAATACCGATGATTTCATCCAGTTCACCGCGACATACCGGGAACAGACTGTGCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCagag                                                                                                                                         <  2:797580/150‑1 (MQ=255)
          tACAATACCGATGATTTCATCCAGTTCACCGCGACATACCGGGAACAGACTGTGCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGGATCATGACGCCgcgcggtga                                                                                                                                        <  1:1248457/150‑6 (MQ=255)
          tACAATACCGATGATTTCATCCAGTTCACCGCGACATACCGGGAACAGACTGTGCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGGATCATGACGCCgcgcggtga                                                                                                                                        <  2:601519/150‑6 (MQ=255)
           aCAATACCGATGATTTCATCCAGTTCACCGCGACATACCGGGAACAGACTGTGCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGGATCATGACGCCgcgcggt                                                                                                                                          >  1:1087027/1‑144 (MQ=255)
           aCAATACCGATGATTTCATCCAGTTCACCGCGACATACCGGGAACAGACTGTGCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGGATCATGACGCCgcgcggt                                                                                                                                          <  2:1087027/147‑4 (MQ=255)
              aTACCGATGATTTCATCCAGATCACCGCGACATACCGGGAACAGACTGTGCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAAcgc                                                                                                                                    <  1:461438/150‑1 (MQ=255)
                  cGATGATATCATCCAGTTCACCGCGACATACCGGGAACAGACTGTGCGGAGAAGAGAGCACTAGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCAAGCATAATACACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGAcc                                                                                                                                >  1:59624/1‑149 (MQ=255)
                        tttCATCCAGTTCACCGCGACATACCGGGAACAGACTGTGCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCg                                                                                                                          >  2:1251160/1‑150 (MQ=255)
                          tCATCCTGTTCACCGCGACAAACCGGGAACAGACTGAGCGGTGTAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTc                                                                                                                        <  2:1469472/150‑1 (MQ=255)
                           cATCCAGTTCACCGCGACATACCGGGAACAGACTGTGCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCa                                                                                                                       >  2:1683141/1‑150 (MQ=255)
                             tCCAGTTCACCGCGACATACCGGGAACAGACTGTGCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGt                                                                                                                     >  1:363618/1‑150 (MQ=255)
                                aGTTCACCGCGACATACCGGGAACAGACTGTGCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACg                                                                                                                  >  2:1145014/1‑150 (MQ=255)
                                aGTTCACCGCGACATACCGGGAACAGACTGTGCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACg                                                                                                                  <  1:1420020/150‑1 (MQ=255)
                                      ccGCGACATACCGGGAACAGACTGTGCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTa                                                                                                            <  2:23456/150‑1 (MQ=255)
                                      ccGCGACATACCGGGAACAGACTGTGCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTa                                                                                                            <  2:1243657/150‑1 (MQ=255)
                                          gACATACCGGGAACCGACTGTGCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATATCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCa                                                                                                        <  1:299194/150‑1 (MQ=255)
                                            cATACCGGGAACAGACTGTGCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATg                                                                                                      >  1:934852/1‑150 (MQ=255)
                                             aTACCGGGAACAGACTTTGCGGTGAAGAGAGCTGTTGCTCGCGGATTTCATCGACCCCGAGATTATAGTCAACCCAGCTTATTTCACCGCGCGTCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGt                                                                                                     <  1:1630345/150‑1 (MQ=255)
                                               aCCGGGAACAGACTGTGCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTaa                                                                                                   <  1:123109/150‑1 (MQ=255)
                                                ccGGGAACAGACTGTGCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAAc                                                                                                  >  2:1691794/1‑150 (MQ=255)
                                                   ggAACAGACTGTGCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGtt                                                                                               <  2:445518/150‑1 (MQ=255)
                                                    gAACAGACTGTGCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCACCCCAGCTTATTTCACCCCGCGGCGTCATGATCCACCGCAGAGAACGCGACGCCCGCGTCAGTACGCCGTTAATCAAGTAACGttc                                                                                              >  2:886225/1‑150 (MQ=255)
                                                    gAACAGACTGTGCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGttc                                                                                              <  2:363618/150‑1 (MQ=255)
                                                     aaCAGACTGTGCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATACCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGttct                                                                                             <  1:211790/150‑1 (MQ=255)
                                                     aaCAGACTGAGCTGTGAAGTGATCATTTTGTCGTGGATTTGATAGACCCCGAGATTAGCGTGAACGCTTCTTATTTCTCAGCGCGGAGTAATGAGCCCTCGCAGAGAAAGCGACTCGAGCGTCTGTACGCCGTTAATCATGTAACGttct                                                                                             <  1:292453/150‑1 (MQ=255)
                                                      aCAGACTGTGCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGttctt                                                                                            >  1:1714310/1‑150 (MQ=255)
                                                       cAGACTGTGCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGttcttc                                                                                           >  1:1266412/1‑150 (MQ=255)
                                                         gACTGTCCGGTGAAGAGAGCAGATGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTAATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACACCATCGTAAGTACGCCGTTAATCATGTAACGttcttctt                                                                                         >  1:657786/1‑150 (MQ=255)
                                                           cTGTGCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCg                                                                                       >  1:473018/1‑150 (MQ=255)
                                                           cTGTGCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCg                                                                                       <  2:1485117/150‑1 (MQ=255)
                                                           cTGTGCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCg                                                                                       <  2:384437/150‑1 (MQ=255)
                                                                cGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCaaa                                                                                  >  1:435437/1‑150 (MQ=255)
                                                                  ttGAAGAGAGCAGTTGCTCGAGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAAAGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACTTTCTTCTTCGGCAAATg                                                                                <  2:1094991/149‑1 (MQ=255)
                                                                    gAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCa                                                                              <  1:1683141/150‑1 (MQ=255)
                                                                      cgcgagCAGTTGCGCGCGGAGTTCAGCGACCCCGAGATTAGCGTCCACCCAGCTTATTTCACCGCGCGGCGTCATGAGCCCGCGCAGAGAACGCGACGCCAGCGTCAGTCCGCCGGTAATCAGGTAACGGTCTTCTGCGGCAAAGGCAcc                                                                            >  1:895766/4‑150 (MQ=255)
                                                                           gCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCg                                                                       <  2:502021/150‑1 (MQ=255)
                                                                            cAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCgg                                                                      <  2:857446/150‑1 (MQ=255)
                                                                                    cgcgGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTACTCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGatctgca                                                              >  2:102123/1‑150 (MQ=255)
                                                                                      cgGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGatcggcatc                                                            >  2:410152/1‑150 (MQ=255)
                                                                                        gATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGatcggcatcgg                                                          >  2:760320/1‑150 (MQ=255)
                                                                                        gATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGatcggcatcgg                                                          >  1:1029965/1‑150 (MQ=255)
                                                                                         aTTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCAGGatcggcatcggc                                                         >  1:1170546/1‑150 (MQ=255)
                                                                                            tCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCggcatcggcatc                                                      >  2:1325823/1‑150 (MQ=255)
                                                                                               tCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCggg                                                   <  1:754122/150‑1 (MQ=255)
                                                                                                cGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGt                                                  <  2:650932/150‑1 (MQ=255)
                                                                                                  aCCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTa                                                <  2:1266412/150‑1 (MQ=255)
                                                                                                    cccGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATc                                              <  2:1376441/150‑1 (MQ=255)
                                                                                                    cccGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATc                                              <  2:405667/150‑1 (MQ=255)
                                                                                                      cGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCgg                                            <  2:1029965/150‑1 (MQ=255)
                                                                                                          aTTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATc                                        >  1:775051/1‑150 (MQ=255)
                                                                                                           ttAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCg                                       <  2:934852/150‑1 (MQ=255)
                                                                                                            tAGCTTCAACCCAGCGTATTTCCCCGGGCGGCGTAATGATCCCGCGCCGAGAACGCGACGCCAGCGGCAGTACTCCGTTAATCATGTAACGTTCTTCTTCGTCAAATGCACCTTCCGTTAGCGGCATCGGCATCGGGTTATCGGAATCGt                                      <  1:88567/150‑1 (MQ=255)
                                                                                                               cGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCt                                   <  2:435437/150‑1 (MQ=255)
                                                                                                                gTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTg                                  >  2:213159/1‑150 (MQ=255)
                                                                                                                         gCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGc                         >  2:437779/1‑150 (MQ=255)
                                                                                                                              tttCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCCCCTTCCGGGATCGGCCTCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGAc                    >  1:230672/1‑150 (MQ=255)
                                                                                                                              tttCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGAc                    <  1:1691794/150‑1 (MQ=255)
                                                                                                                                tctctgcgcgGGATCATGACGCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCGTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGt                  >  1:558916/6‑150 (MQ=255)
                                                                                                                                tCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGt                  <  1:1145014/150‑1 (MQ=255)
                                                                                                                                tCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCATCGGCAAATGCACCTTACGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGt                  <  1:170238/150‑1 (MQ=255)
                                                                                                                                  aCCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCAGCGGGGTATCGGGATCGTGCTGAACATTGGCCTGACGttt                >  1:102426/1‑150 (MQ=255)
                                                                                                                                   ctgcgcgGGATCATGACGCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTc               >  1:1584967/3‑150 (MQ=255)
                                                                                                                                      cgcgGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCg            >  2:1141173/1‑150 (MQ=255)
                                                                                                                                         ggcgTCATGATCCCGCGCAGAGAACGCGCCGCCAGCGTCAGGACGCCGTTAATCATGTAACGGTCTTCTTCGGCACATGCACCTTCCGGGATCGGCATCGGCATCGGGTGATCGGCATCGTGCTGAACATTGGCCTGACGTTgcccgccc         <  1:267642/150‑1 (MQ=255)
                                                                                                                                             tCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCTTGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCa     <  2:59624/150‑1 (MQ=255)
                                                                                                                                             tCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCa     <  1:1325823/150‑1 (MQ=255)
                                                                                                                                               aTGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCaaa   >  1:870129/1‑150 (MQ=255)
                                                                                                                                                tGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAAc  <  2:473018/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CTTTAGCACGTACAATACCGATGATTTCATCCAGTTCACCGCGACATACCGGGAACAGACTGTGCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAAC  >  NC_000913/1900532‑1900825

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