Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,197,294 T→C E783G (GAA→GGA)  glnE ← fused deadenylyltransferase/adenylyltransferase for glutamine synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,197,2940TC98.3% 424.2 / ‑1.5 117E783G (GAA→GGA) glnEfused deadenylyltransferase/adenylyltransferase for glutamine synthetase
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (53/62);  total (54/63)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.80e-01

Only 100 of 120 total aligned reads displayed.

CGCATTTTCTCGCGCATTTCCCGCACTTCCGTTTGCAGAGTTTTACCTTCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTTCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCAAATGCTTCTGCGGATGTCACCAGCATTCCCGCCGCCCCGGACGGACGCAGTCGAGCATCCACTTCATACAAAATGCCGGAAGAGGTACGCGTACTGAACAGATGCATAATG  >  NC_000913/3197147‑3197440
                                                                                                                                                   |                                                                                                                                                  
cgcATTTTCTCGCGCATTTCCCGCACTTCCGTTTGCAGAGTTTTACCTTCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTcccc                                                                                                                                                 <  2:831342/151‑1 (MQ=255)
cgcATTTTCTCGCGCATTTCCCGCACTTCCGTTTGCAGAGTTTTACCTTCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTcccc                                                                                                                                                 <  4:709314/151‑1 (MQ=255)
 gcATTTTCTCGCGCATTTCCCGCACTTCCGTTTGCAGAGTTTTACCTTCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTCCCCa                                                                                                                                                <  2:527101/151‑1 (MQ=255)
  cATTTTCTCGCGCATTTCCCGCACTTCCGTTTGCAGAGTTTTACCTTCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTCCCCa                                                                                                                                                <  2:1151567/150‑1 (MQ=255)
  cATTTTCTCGCGCATTTCCCGCACTTCCGTTTGCAGAGTTTTACCTTCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTCCCCAc                                                                                                                                               >  4:648447/1‑151 (MQ=255)
    ttttCTCGCGCATTTCCCGCACTTCCGTTTGCAGAGTTTTACCTTCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTCCCCACGt                                                                                                                                             >  2:1663829/1‑151 (MQ=255)
     tttCTCGCGCATTTCCCGCACTTCCGTTTGCAGAGTTTTACCTTCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTCCCCACGTc                                                                                                                                            >  2:1381981/1‑151 (MQ=255)
       tctcGCGCATTTCCCGCACTTCCGTTTGCAGAGTTTTACCTTCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTCCCCACGTcc                                                                                                                                           <  2:804160/150‑1 (MQ=255)
          cgcgCATTTCCCGCACTTCCGTTTGCAGAGTTTTACCTTCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTCCCCACGTCCa                                                                                                                                          <  2:586179/148‑1 (MQ=255)
             gcATTTCCCGCACTTCCGTTTGCAGAGTTTTACCTTCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTCCCCACGTCCAGGCCTc                                                                                                                                    <  2:950943/151‑1 (MQ=255)
             gcATTTCCCGCACTTCCGTTTGCAGAGTTTTACCTTCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTCCCCACGTCCAGGCCTc                                                                                                                                    >  4:903438/1‑151 (MQ=255)
               aTTTCCCGCACTTCCGTTTGCAGAGTTTTACCTTCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTCCCCACGTCCAGGCCTCGt                                                                                                                                  >  2:1264755/1‑151 (MQ=255)
                  tCCCGCACTTCCGTTTGCAGAGTTTTACCTTCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTCCCCACGTCCAGGCCTCGttt                                                                                                                                <  4:506229/150‑1 (MQ=255)
                   cccGCACTTCCGTTTGCAGAGTTTTACCTTCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTCCCCACGTCCAGGCCTCGttttt                                                                                                                              <  1:187243/151‑1 (MQ=255)
                    ccGCACTTCCGTTTGCAGAGTTTTACCTTCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCCGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTCCCCACGTCCAGGCCTCGTTTTTc                                                                                                                             <  2:366461/151‑1 (MQ=255)
                       cACTTCCGTTTGCAGAGTTTTACCTTCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTCCCCACGTCCAGGCCTCGTTTTTCTg                                                                                                                           >  4:123652/1‑150 (MQ=255)
                          ttCCGTTTGCAGAGTTTTACCTTCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTCCCCACGTCCAGGCCTCGTTTTTCTGATa                                                                                                                        >  3:101156/1‑150 (MQ=255)
                          ttCCGTTTGCAGAGTTTTACCTTCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTCCCCACGTCCAGGCCTCGTTTTTCTGATa                                                                                                                        >  2:593931/1‑150 (MQ=255)
                           tCCGTTTGCAGAGTTTTACCTTCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTCCCCACGTCCAGGCCTCGTTTTTCTGATaa                                                                                                                       <  4:1554360/150‑1 (MQ=255)
                           tCCGTTTGCAGAGTTTTACCTTAACTCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGAGCGCCCTACACTACACGCGCACGCACCAGCGCCTGATGTCCCCAGGTCCAGGCCTCGTTTTTc                                                                                                                             >  1:52369/1‑144 (MQ=255)
                            ccGTTTGCAGAGTTTTACCTTCACGCGGCAGCGTCATGATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTCCCCACGTCCAGGCCTCGTTTTTCTGATAATc                                                                                                                     <  4:909769/151‑1 (MQ=255)
                            ccGTTTGCAGAGTTTTACCTTCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTCCCCACGTCCAGGCCTCGTTTTTCTGATAATc                                                                                                                     >  2:1449132/1‑151 (MQ=255)
                             cGTTTGCAGAGTTTTACCTTCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTCCCCACGTCCAGGCCTCGTTTTTCTGATAATCg                                                                                                                    <  2:319292/151‑1 (MQ=255)
                              gTTTGCAGAGTTTTACCTTCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTCCCCACGTCCAGGCCTCGTTTTTCTGATAATCgg                                                                                                                   <  4:1540948/151‑1 (MQ=255)
                              gTTTGCAGAGTTTTACCTTCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTCCCCACGTCCAGGCCTCGTTTTTCTGATAATCgg                                                                                                                   >  4:1644483/1‑151 (MQ=255)
                               tttGCAGAGTTTTACCTTCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTCCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGc                                                                                                                  >  4:825840/1‑151 (MQ=255)
                               tttGCAGAGTTTTACCTTCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTCCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGc                                                                                                                  >  4:114113/1‑151 (MQ=255)
                                ttGCAGAGTTTTACCTTCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTCCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGc                                                                                                                  >  4:1628068/1‑150 (MQ=255)
                                 tGCAGAGTTTTACCTTCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTCCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCa                                                                                                                 <  2:1469254/150‑1 (MQ=255)
                                 tGCAGAGTTTTACCTTCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTCCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCa                                                                                                                 <  1:33189/150‑1 (MQ=255)
                                   cAGAGTTTTACCTTCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTTCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCAAAt                                                                                                              >  2:634931/1‑151 (MQ=255)
                                   cAGAGTTTTACCTTCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTCCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCaaa                                                                                                               >  4:177602/1‑150 (MQ=255)
                                      agTTTTACCTTCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTCCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCAAATGCt                                                                                                           >  4:559000/1‑151 (MQ=255)
                                       gTTTTACCTTCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTCCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCAAATGCtt                                                                                                          >  2:1643085/1‑151 (MQ=255)
                                          ttACCTTCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTCCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCAAATGCTTCt                                                                                                        <  2:1000742/150‑1 (MQ=255)
                                            aCCTTCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTCCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCAAATGCTTCTGc                                                                                                      <  3:249594/150‑1 (MQ=255)
                                              cTTCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTCCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCAAATGCTTCTGCgg                                                                                                    >  2:1445396/1‑150 (MQ=255)
                                              cTTCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTCCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCAAATGCTTCTGCgg                                                                                                    >  4:452611/1‑150 (MQ=255)
                                                tCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTCCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCAAATGCTTCTGCGGATg                                                                                                 >  2:618746/1‑151 (MQ=255)
                                                 cACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTCCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCAAATGCTTCTGCGGATGt                                                                                                <  4:1268385/151‑1 (MQ=255)
                                                   cgcgGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTCCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCAAATGCTTCTGCGGATGTc                                                                                               >  4:1400760/1‑150 (MQ=255)
                                                    gcgGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTCCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCAAATGCTTCTGCGGATGTCAc                                                                                             >  2:1524900/1‑151 (MQ=255)
                                                     cgGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTCCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCAAATGCTTCTGCGGATGTCAcc                                                                                            >  4:1091499/1‑151 (MQ=255)
                                                      ggCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTCCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCAAATGCTTCTGCGGATGTCACCa                                                                                           <  2:834760/151‑1 (MQ=255)
                                                        cAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTCCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCAAATGCTTCTGCGGATGTCACCAGc                                                                                         <  1:174794/151‑1 (MQ=255)
                                                        cAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTCCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCAAATGCTTCTGCGGATGTCACCAGc                                                                                         <  1:154122/151‑1 (MQ=255)
                                                        cAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTCCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCAAATGCTTCTGCGGATGTCACCAGc                                                                                         >  1:440391/1‑151 (MQ=255)
                                                         aGCGTTATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTCCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCAAATGCTTCTGCGGATGTCACCAGc                                                                                         >  1:367666/1‑150 (MQ=255)
                                                          gCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTCCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCAAATGCTTCTGCGGATGTCACCAGCAt                                                                                       <  4:1413581/151‑1 (MQ=255)
                                                           cGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTCCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCAAATGCTTCTGCGGATGTCACCAGCAtt                                                                                      >  2:672728/1‑151 (MQ=255)
                                                             tCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTCCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCAAATGCTTCTGCGGATGTCACCAGCATTcc                                                                                    >  4:525840/1‑151 (MQ=255)
                                                              cATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTCCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCAAATGCTTCTGCGGATGTCACCAGCATTccc                                                                                   <  4:253840/151‑1 (MQ=255)
                                                              cATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTCCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCAAATGCTTCTGCGGATGTCACCAGCATTccc                                                                                   <  4:1254048/151‑1 (MQ=255)
                                                                  aTCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTCCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCAAATGCTTCTGCGGATGTCACCAGCATTCccgccg                                                                               >  2:439339/1‑151 (MQ=255)
                                                                  aTCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTCCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCAAATGCTTCTGCGGATGTCACCAGCATTCccgccg                                                                               >  2:369439/1‑151 (MQ=255)
                                                                      cgcgACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTCCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCAAATGCTTCTGCGGATGTCACCAGCATTCCCGCCGcccc                                                                           >  2:1146294/1‑151 (MQ=255)
                                                                       gcgACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTCCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCAAATGCTTCTGCGGATGTCACCAGCATTCCCGCCGCCCcg                                                                          <  4:1021153/151‑1 (MQ=255)
                                                                         gACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTCCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCAAATGCTTCTGCGGATGTCACCAGCATTCCCGCCGCCCcgg                                                                         >  4:1439606/1‑150 (MQ=255)
                                                                         gACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTCCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCAAATGCTTCTGCGGATGTCACCAGCATTCCCGCCGCCCcgg                                                                         >  3:335360/1‑150 (MQ=255)
                                                                          aCGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTCCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCAAATGCTTCTGCGGATGTCACCAGCATTCCCGCCGCCCcggac                                                                       <  4:1059240/151‑1 (MQ=255)
                                                                           cGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTCCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCAAATGCTTCTGCGGATGTCACCAGCATTCCCGCCGCCCcggact                                                                      <  2:43868/151‑2 (MQ=255)
                                                                           cGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTCCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCAAATGCTTCTGCGGATGTCACCAGCATTCCCGCCGCCCcggacg                                                                      >  2:1588022/1‑151 (MQ=255)
                                                                              aCTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTCCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCAAATGCTTCTGCGGATGTCACCAGCATTCCCGCCGCCCCggacggac                                                                   <  3:166812/151‑1 (MQ=255)
                                                                                 gCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTCCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCAAATGCTTCTGCGGATGTCACCAGCATTCCCGCCGCCCCGGACGGACGc                                                                 <  2:113926/150‑1 (MQ=255)
                                                                                 gCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTCCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCAAATGCTTCTGCGGATGTCACCAGCATTCCCGCCGCCCCGGACGGACGCa                                                                <  1:348350/151‑1 (MQ=255)
                                                                                  cGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTCCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCAAATGCTTCTGCGGATGTCACCAGCATTCCCGCCGCCCCGGACGGACGCAg                                                               <  4:1406040/151‑1 (MQ=255)
                                                                                  cGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTCCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCAAATGCTTCTGCGGATGTCACCAGCATTCCCGCCGCCCCGGACGGACGCAg                                                               <  2:1441933/151‑1 (MQ=255)
                                                                                    tCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTCCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCAAATGCTTCTGCGGATGTCACCAGCATTCCCGCCGCCCCGGACGGACGCAGTc                                                             >  4:322117/1‑151 (MQ=255)
                                                                                     cAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTCCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCAAATGCTTCTGCGGATGTCACCAGCATTCCCGCCGCCCCGGACGGACGCAGTc                                                             <  4:1073138/150‑1 (MQ=255)
                                                                                     cAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTCCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCAAATGCTTCTGCGGATGTCACCAGCATTCCCGCCGCCCCGGACGGACGCAGTc                                                             <  4:601190/150‑1 (MQ=255)
                                                                                     cAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTCCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCAAATGCTTCTGCGGATGTCACCAGCATTCCCGCCGCCCCGGACGGACGCAGTc                                                             <  2:1675576/150‑1 (MQ=255)
                                                                                     cAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTCCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCAAATGCTTCTGCGGATGTCACCAGCATTCCCGCCGCCCCGGACGGACGCAGTCg                                                            <  3:33254/151‑1 (MQ=255)
                                                                                      aaaGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTCCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCAAATGCTTCTGCGGATGTCACCAGCATTCCCGCCGCCCCGGACGGACGCAGTCg                                                            >  4:1219664/1‑150 (MQ=255)
                                                                                       aaGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTCCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCAAATGCTTCTGCGGATGTCACCAGCATTCCCGCCGCCCCGGACGGACGCAGTCg                                                            >  4:614459/1‑149 (MQ=255)
                                                                                       aaGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTCCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCAAATGCTTCTGCGGATGTCACCAGCATTCCCGCCGCCCCGGACGGACGCAGTCGAg                                                          <  2:1218130/151‑1 (MQ=255)
                                                                                                gTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTCCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCAAATGCTTCTGCGGATGTCACCAGCATTCCCGCCGCCCCGGACGGACGCAGTCGAGCATCCACtt                                                 <  2:578972/151‑1 (MQ=255)
                                                                                                   aGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTCCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCAAATGCTTCTGCGGATGTCACCAGCATTCCCGCCGCCCCGGACGGACGCAGTCGAGCATCCACTTCAt                                              >  2:796288/1‑151 (MQ=255)
                                                                                                         agATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTCCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCAAATGCTTCTGCGGATGTCACCAGCATTCCCGCCGCCCCGGACGGACGCAGTCGAGCATCCACTTCATACaaa                                         >  2:891510/2‑150 (MQ=255)
                                                                                                          gATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTCCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCAAATGCTTCTGCGGATGTCACCAGCATTCCCGCCGCCCCGGACGGACGCAGTCGAGCATCCACTTCATACAAAAt                                       <  2:1442672/151‑1 (MQ=255)
                                                                                                             cGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTCCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCAAATGCTTCTGCGGATGTCACCAGCATTCCCGCCGCCCCGGACGGACGCAGTCGAGCATCCACTTCATACAAAATGcc                                    >  4:45567/1‑151 (MQ=255)
                                                                                                              gCCGTACACTACACGCGCACGCACCAGCGCCTGATGTCCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCAAATGTTTCTGCGGATGTCACCAGCATTCCCGCCGCCCCGGACGGACGCAGTCGAGCATCCACTTCATACAAAATGCCg                                   <  3:253700/151‑1 (MQ=255)
                                                                                                               ccGTACACTACACGCGCACGCACCAGCGCCTGATGTCCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCAAATGCTTCTGCGGATGTCACCAGCATTCCCGCCGCCCCGGACGGACGCAGTCGAGCATCCACTTCATACAAAATGCCgg                                  >  2:419895/1‑151 (MQ=255)
                                                                                                               ccGTACACTACACGCGCACGCACCAGCGCCTGATGTCCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCAAATGCTTCTGCGGATGTCACCAGCATTCCCGCCGCCCCGGACGGACGCAGTCGAGCATCCACTTCATACAAAATGCCgg                                  <  2:2183/151‑1 (MQ=255)
                                                                                                                 gTACACTACACGCGCACGCACCAGCGCCTGATGTCCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCAAATGCTTCTGCGGATGTCACCAGCATTCCCGCCGCCCCGGACGGACGCAGTCGAGCATCCACTTCATACAAAATGCCGGaa                                <  2:711425/151‑1 (MQ=255)
                                                                                                                 gTACACTACACGCGCACGCACCAGCGCCTGATGTCCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCAAATGCTTCTGCGGATGTCACCAGCATTCCCGCCGCCCCGGACGGACGCAGTCGAGCATCCACTTCATACAAAATGCCGGa                                 >  4:1227714/1‑150 (MQ=255)
                                                                                                                      ctacacGCGCACGCACCAGCGCCTGATGTCCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCAAATGCTTCTGCGGATGTCACCAGCATTCCCGCCGCCCCGGACGGACGCAGTCGAGCATCCACTTCATACAAAATGCCGGAAGAGGt                           >  4:538497/1‑151 (MQ=255)
                                                                                                                        acacGCGCACGCACCAGCGCCTGATGTTCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCAAATGCTTCTGCGGATGTCACCAGCATTCCCGCCGCCCCGGACGGACGCAGTCGAGCATCCACTTCATACAAAATGCCGGAAGAGGTAc                         <  4:635789/151‑1 (MQ=255)
                                                                                                                        acacGCGCACGCACCAGCGCCTGATGTCCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCAAATGCTTCTGCGGATGTCACCAGCATTCCCGCCGCCCCGGACGGACGCAGTCGAGCATCCACTTCATACAAAATGCCGGAAGAGGTAc                         <  3:105856/151‑1 (MQ=255)
                                                                                                                         cacGCGCACGCACCAGCGCCTGATGTCCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCAAATGCTTCTGCGGATGTCACCAGCATTCCCGCCGCCCCGGACGGACGCAGTCGAGCATCCACTTCATACAAAATGCCGGAAGAGGTAc                         <  2:624361/150‑1 (MQ=255)
                                                                                                                           cgcgCACGCACCAGCGCCTGATGTCCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCAAATGCTTCTGCGGATGTCACCAGCATTCCCGCCGCCCCGGACGGACGCAGTCGAGCATCCACTTCATACAAAATGCCGGAAGAGGTAcgcg                      <  2:1579038/151‑1 (MQ=255)
                                                                                                                               cacgcacCAGCGCCTGATGTCCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCAAATGCTTCTGCGGATGTCACCAGCATTCCCGCCGCCCCGGACGGACGCAGTCGAGCATCCACTTCATACAAAATGCCGGAAGAGGTACGCGTAc                   >  2:773133/1‑150 (MQ=255)
                                                                                                                               cacgcacCAGCGCCTGATGTCCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCAAATGCTTCTGCGGATGTCACCAGCATTCCCGCCGCCCCGGACGGACGCAGTCGAGCATCCACTTCATACAAAATGCCGGAAGAGGTACGCGTACt                  >  4:1551616/1‑151 (MQ=255)
                                                                                                                                acgcacCAGCGCCTGATGTCCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCAAATGCTTCTGCGGATGTCACCAGCATTCCCGCCGCCCCGGACGGACGCAGTCGAGCATCCACTTCATACAAAATGCCGGAAGAGGTACGCGTACTg                 >  4:795289/1‑151 (MQ=255)
                                                                                                                                 cgcacCAGCGCCTGATGTCCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCAAATGCTTCTGCGGATGTCACCAGCATTCCCGCCGCCCCGGACGGACGCAGTCGAGCATCCACTTCATACAAAATGCCGGAAGAGGTACGCGTACTg                 <  2:1069551/150‑1 (MQ=255)
                                                                                                                                  gcacCAGCGCCTGATGTCCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCAAATGCTTCTGCGGATGTCACCAGCATTCCCGCCGCCCCGGACGGACGCAGTCGAGCATCCACTTCATACAAAATGCCGGAAGAGGTACGCGTACTGag               <  3:389696/151‑2 (MQ=255)
                                                                                                                                     ccAGCGCCTGATGTCCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCAAATGCTTCTGCGGATGTCACCAGCATTCCCGCCGCCCCGGACGGACGCAGTCGAGCATCCACTTCATACAAAATGCCGGAAGAGGTACGCGTACTGAACa             <  2:566430/150‑1 (MQ=255)
                                                                                                                                     ccAGCGCCTGATGTCCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCAAATGCTTCTGCGGATGTCACCAGCATTCCCGCCGCCCCGGACGGACGCAGTCGAGCATCCACTTCATACAAAATGCCGGAAGAGGTACGCGTACTGAACAg            <  4:1681340/151‑1 (MQ=255)
                                                                                                                                      cAGCGCCTGATGTCCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCAAATGCTTCTGCGGATGTCACCAGCATTCCCGCCGCCCCGGACGGACGCAGTCGAGCATCCACTTCATACAAAATGCCGGAAGAGGTACGCGTACTGAAc              >  4:200876/1‑148 (MQ=255)
                                                                                                                                       aGCGCCTGATGTCCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCAAATGCTTCTGCGGATGTCACCAGCATTCCCGCCGCCCCGGACGGACGCAGTCGAGCATCCACTTCATACAAAATGCCGGAAGAGGTACGCGTACTGAACAGAt          <  4:1447203/151‑1 (MQ=255)
                                                                                                                                               aTGTCCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCAAATGCTTCTGCGGATGTCACCAGCATTCCCGCCGCCCCGGACGGACGCAGTCGAGCATCCACTTCATACAAAATGCCGGAAGAGGTACGCGTACTGAACAGATGCATAATg  <  2:1593252/151‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CGCATTTTCTCGCGCATTTCCCGCACTTCCGTTTGCAGAGTTTTACCTTCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTTCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCAAATGCTTCTGCGGATGTCACCAGCATTCCCGCCGCCCCGGACGGACGCAGTCGAGCATCCACTTCATACAAAATGCCGGAAGAGGTACGCGTACTGAACAGATGCATAATG  >  NC_000913/3197147‑3197440

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