Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,306,090 G→T 6.1% Q170K (CAA→AAA)  yjcS ← metallo‑beta‑lactamase superfamily protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,306,0900GT 6.1% 311.9 / 4.5 109Q170K (CAA→AAA) yjcSmetallo‑beta‑lactamase superfamily protein
Reads supporting (aligned to +/- strand):  ref base G (39/62);  new base T (5/3);  total (44/65)
Fisher's exact test for biased strand distribution p-value = 2.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 109 total aligned reads displayed.

TTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATCCTGGCCGCGCACCTGATACATTTTATCGGTGAC  >  NC_000913/4305942‑4306237
                                                                                                                                                    |                                                                                                                                                   
ttttCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGc                                                                                                                                                    <  2:234925/150‑1 (MQ=255)
 tttCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCg                                                                                                                                                   <  1:1738882/150‑1 (MQ=255)
 tttCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCg                                                                                                                                                   <  1:1341950/150‑1 (MQ=255)
   tCGCTGATGTCTTCGTACATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGCTTTCACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGGCGGTGGGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGa                                                                                                                                                 <  2:1943851/150‑1 (MQ=255)
    cGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGAc                                                                                                                                                <  1:1007118/150‑1 (MQ=255)
       tGATGGCTTCGTCCGTAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGATTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGAAAAAATCGGTTTTTGCGGACGAt                                                                                                                                             >  1:1055922/1‑150 (MQ=255)
          tGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGAAGCt                                                                                                                                          >  2:284646/1‑150 (MQ=255)
            gCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTgg                                                                                                                                        <  2:571405/150‑1 (MQ=255)
              ttCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGaa                                                                                                                                      <  2:248166/150‑1 (MQ=255)
               tCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGTTTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAg                                                                                                                                     <  1:1507380/150‑1 (MQ=255)
               tCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTATTGAGGACGATGCTGGAAg                                                                                                                                     >  2:1618799/1‑150 (MQ=255)
                 gTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTa                                                                                                                                   >  1:97805/1‑150 (MQ=255)
                    cATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGATTTTGCGGACGATGCTGGAAGTAAAg                                                                                                                                >  1:235703/1‑150 (MQ=255)
                     aTAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAgg                                                                                                                               <  2:1201184/150‑1 (MQ=255)
                       aaaGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTc                                                                                                                             >  1:1083957/1‑150 (MQ=255)
                             tGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGc                                                                                                                       >  1:631250/1‑150 (MQ=255)
                             tGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGATGGAAGTAAAGGTCAAGTGc                                                                                                                       >  2:1629283/1‑150 (MQ=255)
                               cAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCgg                                                                                                                     <  1:425444/150‑1 (MQ=255)
                               cAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGGGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCgg                                                                                                                     <  2:646374/150‑1 (MQ=255)
                                aGGGGCAATCACCTGAACTTTGCAGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGc                                                                                                                    <  1:1497970/150‑1 (MQ=255)
                                      aaTCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGCTAATGTCTTTCACGCCACCATAGTGGTCGGTGTGGGTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGCTGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGc                                                                                                              >  1:71843/1‑150 (MQ=255)
                                      aaTCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGc                                                                                                              <  1:536648/150‑1 (MQ=255)
                                       aTCACCTGAACTTTGCCGGATTTAAGATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCg                                                                                                             <  2:175040/150‑1 (MQ=255)
                                       aTCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGATTTTGAATACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCg                                                                                                             >  2:423911/1‑150 (MQ=255)
                                          aCCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCa                                                                                                          <  1:881576/150‑1 (MQ=255)
                                          aCCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGGCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCa                                                                                                          <  1:1843652/150‑1 (MQ=255)
                                                 cTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGgcggcg                                                                                                   <  2:1678946/150‑1 (MQ=255)
                                                 cTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCAAAGTGGGCGGTGGGGCTGTGAGTGTAGATAACGGCATCACTCGGTTTTTGCGGACGATGCTGGAAGTAATGGTCAAGTGCGGCTTTTGCGGCAGgcggcg                                                                                                   <  1:332703/150‑1 (MQ=255)
                                                     gCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCAc                                                                                               >  1:1732127/1‑150 (MQ=255)
                                                      ccGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCAcc                                                                                              <  2:908801/150‑1 (MQ=255)
                                                      ccGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCAcc                                                                                              <  1:1374992/150‑1 (MQ=255)
                                                         gATTTAACATCGGCTTCAGAGATAATGCCTTGCCCGCCACCATAGGGGTCGGTGGTGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGc                                                                                           <  2:132509/150‑1 (MQ=255)
                                                          aTTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCg                                                                                          <  2:1093202/150‑1 (MQ=255)
                                                          aTTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCg                                                                                          >  1:1378786/1‑150 (MQ=255)
                                                           tttAACATCGGGTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGTACAAGTGCGGCTTTTGCTGCAGGCGTCGTCACCAGCgg                                                                                         >  1:261551/1‑150 (MQ=255)
                                                             tAACATCGGCTGCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGt                                                                                       <  2:953375/150‑1 (MQ=255)
                                                              aaCATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTc                                                                                      <  1:1230960/150‑1 (MQ=255)
                                                                  tCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATCTCGGAAAAAAACGTATAATGAAAAAAAACATGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGCCGGCACCAGCGGGTCGATa                                                                                  >  2:199141/1‑150 (MQ=255)
                                                                  tCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATa                                                                                  <  2:865886/150‑1 (MQ=255)
                                                                  tCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATa                                                                                  <  1:750451/150‑1 (MQ=255)
                                                                  tCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATa                                                                                  <  2:909216/150‑1 (MQ=255)
                                                                  tCGGCTGCAGAGATAATGCCTATCACGCCACCATAGTGGGCGGTGTGGCTGTGAGTGTCGATAACGGCAACACTCGGTGTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATa                                                                                  <  1:353157/150‑1 (MQ=255)
                                                                   cGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATaa                                                                                 >  2:707397/1‑150 (MQ=255)
                                                                   cGGCTTCAGAGAGAATGCCTTTCTCGCCACCATAGTGGGCGGGGTGGCTGTGAGGGGAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGATGTAAATGTCAAGGGCGGCTCTTGCGGCAGGCGGCGTCACCAGCGGGTCGATaa                                                                                 <  1:984397/150‑1 (MQ=255)
                                                                    ggCTTCAGAGATAATGCCTTTCCCGCCACCATAGTGGGCGGTGTGGCTGTGAGTGGAGATAACGGCAACAATCGGGTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAAc                                                                                <  1:184588/150‑1 (MQ=255)
                                                                    ggCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAAc                                                                                <  1:151940/150‑1 (MQ=255)
                                                                    ggCTTCAGAGATAATGCCTTTCACGCCACCATACCGGGCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGAGGACGATGCTGGAAGTAAAGGTCCAGTGCGGCTTTTGCGGCAGGCGTCGTCACCAGCGGGTCGATAAc                                                                                <  2:1794546/150‑1 (MQ=255)
                                                                      cTTCAGAGATAATGCCTTTAACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGAAAAAAACGTATAATTAAAAAAAAACTGGAAGTAAAGGTCAAGTGAGCCTTTTGCGGCAGGCGGCGTCACCTGCGGGTCGATAACaa                                                                              >  2:1214024/1‑150 (MQ=255)
                                                                         cAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACaataa                                                                           >  1:1346029/1‑150 (MQ=255)
                                                                         cAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACaataa                                                                           >  1:1442752/1‑150 (MQ=255)
                                                                         cAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGAAGATAACGGCAACAATCGGTTATTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTAGCGGCAGGCGGCGTCACCAGCGGGTCGCTAACaataa                                                                           >  1:1601022/1‑150 (MQ=255)
                                                                             gaTAATGCCTTTCAAGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTTAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGcc                                                                       <  2:1463944/150‑1 (MQ=255)
                                                                               tAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCtt                                                                     >  1:867556/1‑150 (MQ=255)
                                                                                  tGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTtct                                                                  <  2:1429733/150‑1 (MQ=255)
                                                                                  tGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTtct                                                                  <  2:32344/150‑1 (MQ=255)
                                                                                  tGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTtct                                                                  >  1:33885/1‑150 (MQ=255)
                                                                                   gCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTtctc                                                                 >  1:571374/1‑150 (MQ=255)
                                                                                    ccTTTCACGCCACCATAGTGGTCTGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTTCGGACGATGCTGGAAGTAAATGTCATGTGCGGCTTTTGTGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCg                                                                <  1:1226539/150‑1 (MQ=255)
                                                                                    ccTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAAAAAACGTATTTTTAAGACTATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCg                                                                >  2:262594/1‑150 (MQ=255)
                                                                                    ccTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGAAAACAACGTTTAATTAAAAAAATAATGGAAGAAAAGGTCAAGTGCGGCTTTTGCGGCAGGCTGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCg                                                                >  2:1384287/1‑150 (MQ=255)
                                                                                     ctggcACGCCAGCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGc                                                               <  2:1803464/146‑1 (MQ=255)
                                                                                     cTTTCACGCCACCATAGTGGTCGGTGTGTCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGAGGCTTTTGCGGCAGGCGGCTTCACCAGCGGGGAGATAACAATATTGCCTTTCTCGc                                                               >  1:346620/1‑150 (MQ=255)
                                                                                     cTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGc                                                               <  2:723674/150‑1 (MQ=255)
                                                                                      tttCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGcc                                                              >  1:1615392/1‑150 (MQ=255)
                                                                                        tCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGATTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCt                                                            >  2:1683818/1‑150 (MQ=255)
                                                                                         cACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTATAGGTCATGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTc                                                           <  2:1329613/150‑1 (MQ=255)
                                                                                         cACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAAGGCCTTTCTCGCCCTc                                                           <  2:631250/150‑1 (MQ=255)
                                                                                               accaTAGTGGTCTGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCaacgaa                                                     <  2:164652/150‑1 (MQ=255)
                                                                                                ccaTAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCTACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCaacgaac                                                    <  1:768602/150‑1 (MQ=255)
                                                                                                ccaTAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCaacgaac                                                    <  2:1436656/150‑1 (MQ=255)
                                                                                                   tAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAACAATGCCTTTCTCGCCCTCAACGAACGTa                                                 <  2:2034575/150‑1 (MQ=255)
                                                                                                        gTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGtt                                            <  1:1930953/150‑1 (MQ=255)
                                                                                                         tCGGTGTGGCTGGGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTGTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTa                                           <  1:1618800/150‑1 (MQ=255)
                                                                                                          ctgtgtGGCTGTGAGTGTAGATATCGGCATCAATCGGTTTTTGCGGACGATTCTGGAAGTAAATGTCAAGTTCGGCTTTTGCGGGAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTag                                          <  1:211564/148‑1 (MQ=255)
                                                                                                          cGGTGTGGCTGTGAGTGTAGATAACGGAAAACAACGTATAATTAAAAAAAAAATGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTag                                          >  2:1555778/1‑150 (MQ=255)
                                                                                                            gtgtGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTagag                                        <  2:2007447/150‑1 (MQ=255)
                                                                                                            gtgtGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTagag                                        <  2:936484/150‑1 (MQ=255)
                                                                                                                  cTGTGAGTGTAGATACCGGCAACAATCGGTTTTTGCGGACGATGGTGGTAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATcc                                  <  1:1285568/150‑1 (MQ=255)
                                                                                                                  cTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATcc                                  >  1:2010189/1‑150 (MQ=255)
                                                                                                                    gtgAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGCTAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATCCTg                                >  1:1309059/1‑150 (MQ=255)
                                                                                                                     tgAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATATCAATAATGCCTTTCTCGCACTCATCGAACGTAATGTTAGAGATATCCTgg                               >  1:388759/1‑150 (MQ=255)
                                                                                                                        gtgtAGATAACGGCAACAATCGGTTTTTGCTGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATCCTGGCcg                            <  1:423911/150‑1 (MQ=255)
                                                                                                                           tAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATATTGCCGTTCTCGCCCTCAACGAACGTAATGTGAGCTATATCCTGGCCgcgc                         >  2:1937376/1‑150 (MQ=255)
                                                                                                                           tAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATCCTGGCCgcgc                         >  1:1116521/1‑150 (MQ=255)
                                                                                                                            aGATAACGGCAACAATCGGTTTTTGCGGACGTTGCTGGAAGTAAAGGTCAAGTGCGGCGTTTGCGGCAGGCGGCGTCACCAGCGGGGCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATCCTGGCCGCGCa                        <  1:1769994/150‑1 (MQ=255)
                                                                                                                            aGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATCCTGGCCGCGCa                        >  2:1130280/1‑150 (MQ=255)
                                                                                                                            aGAAAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATCCTGGCCGCGCa                        <  1:167034/150‑1 (MQ=255)
                                                                                                                             gATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATCCTGGCCGCGCAc                       >  1:1485583/1‑150 (MQ=255)
                                                                                                                               tAACGGCAACAATCGGTTTTTGCGGACGCTGCTGGAAGTAAAGGTCCAGTGCGGCTTTTGCGGCAGGCGGCGTCCCCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCCACGACCGGAATGTTAGAGATATCCTGGCCGCGCACCt                     >  2:1569285/1‑150 (MQ=255)
                                                                                                                               tAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAAAATTAATGCCTTTCTCGCCCTCAACGAACGTAATGTTTGAGATATCATGGCCGCGCACCt                     >  2:836133/1‑150 (MQ=255)
                                                                                                                                 aCGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATCCTGGCCGCGCACCTGa                   >  1:650044/1‑150 (MQ=255)
                                                                                                                                  cGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATCCTGGCCGCGCACCTGAt                  <  1:1343889/150‑1 (MQ=255)
                                                                                                                                    gCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATCCTGGCCGCGCACCTGATAc                >  1:1459009/1‑150 (MQ=255)
                                                                                                                                       acaaTCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTTCGGCTTTTGCGGCAGGCTGCGTCACCAGCGGTTCGATAACAATAATGCATTTCTCTCACTCAACGAACGTAATGTTAGAGATATCCTGGCCTCGCACCTGAAACAtt             >  1:1737316/1‑150 (MQ=255)
                                                                                                                                         caTCGGTTTTTGCGGCCGATGCTGGAAGTAAAGGGCAAGTGCGGCTTTTGCGGCAGGCGGCGTCCCCAGCGGGTCGCTAAAAATCATGCCTTTCTCGCCCTCCCCGAACGTAATGGTAGAGAGATCCTGGCCGCGCCCCTGCTACAtgtt           >  1:178733/2‑150 (MQ=255)
                                                                                                                                         aaTCGGTTTTTTGGGTCGATGCTGGAAGTAAAGGTCACGTGCGGAGATTGCGGCAGGCGGCGTCACCAGCGGGTCGAAAACAATACTGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATCCTGGCCTCGCACCTGATACAtttt           <  1:336451/150‑1 (MQ=255)
                                                                                                                                            cGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATCCTGGCCGCGCACCTGATACATTTTATc        <  1:284646/150‑1 (MQ=255)
                                                                                                                                                ttttTCGGACGATGCTGGAAGTTAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATCCTGGCCGCGCACCTGATACATTTTATCGGta    <  2:1795529/150‑2 (MQ=255)
                                                                                                                                                ttttGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATCCTGGCCGCGCACCTGATACATTTTATCGGTg    <  2:1346028/150‑1 (MQ=255)
                                                                                                                                                  ttGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATCCTGGCCGCGCACCTGATACATTTTATCGGTGAc  <  1:1092620/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATCCTGGCCGCGCACCTGATACATTTTATCGGTGAC  >  NC_000913/4305942‑4306237

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: