Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,208,622 T→G L219* (TTA→TGA)  yehP → VMA domain putative YehL ATPase stimulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,208,6220TG96.2% 174.3 / ‑3.1 53L219* (TTA→TGA) yehPVMA domain putative YehL ATPase stimulator
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base G (24/27);  total (24/29)
Fisher's exact test for biased strand distribution p-value = 4.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.33e-01

TTGAAGAAATTATGGCTCGACTGGCAAAGGAAGTTCGTCAGGCTTTTTCTGGTGTCCGCGATCGCCGTCGCCGTTCATTTATTCCACTGGCGCGAAACTTTGATTTCAAAAGTACTCTGCGCGCCAACCTGCAACACTGGCACCCGCAACACGGCAAGTTGTATATCGAATCCCCCCGCTTTAACAGCCGCATTAAACGCCAAAGCGAACAATGGCAACTGGTCTTACTGGTTGATCAAAGCGGATCGATGGTCGATTCGGTGATCCACTCTGCGGTGATGGCGGCCTGTTTGTGGCAGTTACCCGGCATTCGTACCCATCTGGTGGCGTTTGACACAAGCGTCGTTGATCTCACGGCAGACGTTGCCGATCCGGTAGAGTTATTAATGAAAGTACAGTTGGGCGGCGGGACCAATATCGCCAGTGCCGTGGAGTATGGTCGGCAACTTATTGAA  >  NC_000913/2208397‑2208851
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                     
ttgAAGAAATTATGGCCCGACTGGCAAAGGAAGTTCGTCAGGCTTTTTCTGGTGTCCGCGATCGCCGTCGCCGTTCATTTATTCCACTGGCGCGAAACTTTGATTTCAAAAGTACTCTGCGCGCCAACCTGCAACACTGGCACCCGCAACACGGCAAGTTGTATATCGAATCCCCCCGCTTTAACAGCCGCATTAAACGCCAAAGCGAACAATGGCAACTGGTCTGACTgg                                                                                                                                                                                                                                  <  1:140376/231‑1 (MQ=255)
     taaaTTATGGCTCGACTGGCAAAGGAATTTCGTCCGGCTTTTTCTGGTGTCCGCGTTCCACGTCCCCCTTCCTTTATTCCCCTGGCCCGAAACTTTGATTTCAAAAGTACTCTGCCCGCCAACCTGCCACCCTGGCACCCGCAACACGGCAAGTTGTATATCGAATCCCCCCGCTTTAACAGCCGCATTAAACGCCAAAGCGAACAATGGCAACTGGTCTGACTGGTTGATCAAAg                                                                                                                                                                                                                        <  2:167659/235‑1 (MQ=255)
              gCTCGACTGGAAAGGGAGTTTCGTCAGGCTTTTTCTGGTGTCCGCGATCGCCGTCGCCTTTCATTTATTCCACTGGCGCGAACCTTTGATTCCAAAAGTACTCTGCGCGCCAACCTGCAACACTGGCACCCGCAACACGGCAAGTTGTATATCGAATCCCCCCGCTTTAACAGCCGCATTAAACGCCAAAGCGAACAATGGCAACTGGTCTGACTGGTTGATCAAAGCg                                                                                                                                                                                                                      <  2:509598/229‑1 (MQ=255)
                  gACTGGCAAAGGAAGTTCGTCAGGCTTTTTCTGGTGTCCGCGATCGCCGTCGCCGTTCATTTATTCCACTGGCGCGAAACTTTGATTTCAAAAGTACTCTGCGCGCCAACCTGCAACACTGGCACCCGCAACACGGCAAGTTGTATATCGAATCCCCCCGCTTTAACAGCCGCATTAAACGCCAAAGCGAACAATGGCAACTGGTCTGACTGGTTGATCAAAGCGGATCGATGGt                                                                                                                                                                                                            <  1:641597/235‑1 (MQ=255)
                   aCTGGCAAAAGAAGTTCGTCAGGCTTTTTCTGGTGTCCGCGATCGCCGTCGCCGTTCATTTATTCCACTGGCGCGAAACTTTGATTTCAAAAGTACTCTGCGCGCCAAACTGCAACACTGGCACCCGCAACACGGCAAGTTGTATATCCAATCCCCCCGCGTTAACAGCCCCATTTAAAGCCCAAGAGCACAAAGGAAACAGGGCTGACTGGTTTGTCAAAACGGacagaacgtcg                                                                                                                                                                                                          >  2:703982/1‑226 (MQ=255)
                                 tgcGTCAGGCTTTGTCTGGTGTCCGCGATCGCCGTCGCCGTTCATTTATTCCACTGGCGCGAAACTTTGATTTCAAAAGTACTCTGCGCGCAAACCTGCAACACTGGCACCCGCAACACGGCAAGTTGTATATCGAATCCCCCCGCTTTAACAGCCGCATTAAACGCCAAAGCGAACAATGGCAACTGGTCTGACTGGTTGATCAAAGCGGATCGATGGTCGATTCGGTGAt                                                                                                                                                                                                <  1:353634/230‑1 (MQ=255)
                                  tCGTCAGGCTTTTTCTGGTGTCCGCGATCGCCGTCGCCGTTCATTTATTCCACTGGCGCGAAACTTTGATTTCAAAAGTACTCTGCGCGCCAACCTGCAACACTGGCACCCGCAACACGGCAAGTTGTATATCGAATCCCCCCGCTTTAACAGCCGCATTAAACGCCAAAGCGAACAATGGCAACTGGTCTGACTGGTTGATCAAAGCGGATCGATGGTCGATTCGGTGATc                                                                                                                                                                                               >  1:329335/1‑232 (MQ=255)
                                  tCGTCAGGCTTTTTCTGGGTTCCGAGATCGCCGTCGCCGTTCATTTATTCCACTGGCGCGAAACTTTGATTTCAAAAGTACTCTGCGCGCCAACCTGCAACACTGGCACCCGCAACACGGCAAGTTGTATATCGAATCCCCCCGCTTTAACAGCCGCATTAAACGCCAAAGCGAACAATGGCAACTGGTCTGACTGGTTGATCAAAGCGGATCGATGGTCGATTCGGTGATCCa                                                                                                                                                                                             <  1:368063/234‑1 (MQ=255)
                                     tCAGGCTTTTTCTGGTGTCCGCGATCGCCGTCGCCGTTCATTTATTCCACTGGCGCGAAACTTTGATTTCAAAAGTACTCTGCGCGCCAACCTGCAACACTGGCACCCGCAACACGGCAAGTTGTATATCGAATCCCCCCGCTTTAACAGCCGCATTAAACGCCAAAGCGAACAATGGCAACTCGTCTGACTCGTTGATCCAAACGGATCGATGGTCGATTCAGTGATCCAcccta                                                                                                                                                                                        >  2:301611/1‑235 (MQ=255)
                                              ttCTGGTGTCCGCGATCGCCGTCGCCGTTCATTTATTCCACTGGCGCGAAACTTTGATTTCAAAAGTACTCTGCGCGCCAACCTGCAACACTGGCACCCGCAACACGGCAAGTTGTATATCGAATCCCCCCGCTTTAACAGCCGCATTAAACGCCAAAGCGAACAATGGCAACTGGTCTGACTGGTTGATCAAAGCGGATCGATGGTCGATTCGGTGATCCACTCTGCGGTGATg                                                                                                                                                                                <  1:543120/235‑1 (MQ=255)
                                              gtCTGGTGTCCGCGATCGCTGTCGCCTTTCATTTATTCCACTGGCGCGAAACTTTGATTTCAAAAGTACTCTGCGCGCCAACCTGCAACACTGGCACCCGCAACACGGCAAGTTGTATATCGAATCCCCCCGCTTTAACAGCCGCATTAAACGCCAAAGCGAACAATGGCAACTGGTCTGACTGGTTGATCAAAGCGGATCGATGGTCGATTCGGTGATCCACTCTGCGGTGATgg                                                                                                                                                                               <  2:383010/235‑1 (MQ=255)
                                                  ggTTTCCGAGATCGCGGTCACCGTTCATTTATCCCACGGGCGCGAAACTTTGATTTCAAAAGTACTCTGCGCGCCAACCTGCAACACTGGCACCCGCAACACGGCAAGTTGTATATCGAATCCCCCCGCTTTAACAGCCGCATTAAACGCCAAAGCGAACAATGGCAACTGGTCTGACTGGTTGATCAAAGCGGATCGATGGTCGATTCGGTGATCCACTCTGCGGTGATggcggc                                                                                                                                                                           <  2:355863/236‑1 (MQ=255)
                                                  ggTGTCCGCAATCGCCGTCGCCGTTCATTTATTCCACTGGCGCGAAACTTTGATTTCAAAAGTACTCTGCGCGCCAACCTGCAACACTGGCACCCGCAACACGGCAAGTTGTATATCGAATCCCCCCGCTTTAACAGCCGCATTAAACGCCAAAGCGAACAATGGCAACTGGTCTTACTGGTTGATCAAAGCGGATCGATGGTCGATTCGGTGATCCACTCTGCGGTGATggcggc                                                                                                                                                                           <  1:306942/236‑1 (MQ=255)
                                                   gtgtCCGCGATCGCCGTCGCCGTTCATTTATTCCACTGGCGCGAAACTTTGATTTCAAAAGTACTCTGCGCGCCAACCTGCAACACTGGCACCCGCAACACGGCAAGTTGTATATCGAATCCCCCCGCTTTAACAGCCGCATTAAACGCCAAAGCGAACAATGGCAACTGGTCTGACTGGTTGATCAAAGCGGATCGATGGTCGATTCGGTGATCCACTCTGCGGGGATgg                                                                                                                                                                               >  1:376007/1‑231 (MQ=255)
                                                    tgtCCGCGATCGCCGTCGCCGTTCATTTATTCCACTGGCGCGAAACTTTGATTTCAAAAGTACTCTGCGCGCCAACCTGCAACACTGGCACCCGCAACACGGCAAGTTGTATATCGAATCCCCCCGCTTTAACAGCCGCATTAAACGCCAAAGCGAACAATGGCAACTGGTCTGACTGGTTGATCAAAGCGGATCGATGGTCGATTCGGTGATCCACTCTGCGGTGATGGCGGCCt                                                                                                                                                                         <  1:429129/236‑1 (MQ=255)
                                                                     gccgtTCATTTATTCCACTGGCGCGAAACTTTGATTTCAAAAGTACTCTGCGCGCCAACCTGCAACACTGGCACCCGCAACACGGCAAGTTGTATATCGAATCCCCCCGCTTTAACAGCCGCATTAAACGCCAAAGCGAACAATGGCAACTGGTCTGACTGGTTGATCAAAGCGGATCGATGGTCGATTCGGTGATCCACTCTGCGGTGATggcggc                                                                                                                                                                           <  1:349489/217‑1 (MQ=255)
                                                                     gccgtTCATTTATTCCACTGGCGCGAAACTTTGATTTCAAAAGTACTCTGCGCGCCAACCTGCAACACTGGCACCCGCAACACGGCAAGTTGTATATCGAATCCCCCCGCTTTAACAGCCGCATTAAACGCCAAAGCGAACAATGGCAACTGGTCTGACTGGTTGATCAAAGCGGATCGATGGTCGATTCGGTGATCCACTCTGCGGTGATggcggc                                                                                                                                                                           >  2:349489/1‑217 (MQ=255)
                                                                       cgtTCATTTATTCCCCTGGCGCGAAACTTAGATTTCAAAAGTGCTCTGCGCACCAACCTGCCACACTGGCCCGCGCAACACGGCAAGTTGTATATCGAATCCCCCCGCTTTAACAGCCGCATTAAACGCCAAAGCGAACAATGGCAACTGGTCTGACTGGTTGATCAAAGCGGATCGATGGTCGATTCGGTGATCCACTCTGCGGTGATGGCGGCCTGTTTGTGGCAGTTACCCg                                                                                                                                                       <  2:329335/235‑1 (MQ=255)
                                                                                      cTGGCGCGAAACTTTGATTTCAAAAGTACTCTGCGCGCCAACCTGCAACACTGGCTCCCGCAACACGGCAAGTTGTATATCGAATCCCCCCGCTTTAACAGCCGCATTAAACGCCAAAGCGAACAATGGCAACTGGTCTGACTGGTTGATCAAAGCGGATCGATGGTCGATTCGGTGAt                                                                                                                                                                                                >  2:707029/1‑179 (MQ=255)
                                                                                      cTGGCGCGAAACTTTGATTTCAAAAGTACTCTGCGCGCCAACCTGCAACACTGGCACCCGCAACACGGCAAGTTGTATATCGAATCCCCCCGCTTTAACAGCCGCATTAAACGCCAAAGCGAACAATGGCAACTGGTCTGACTGGTTGATCAAAGCGGATCGATGGTCGATTCGGTGAt                                                                                                                                                                                                <  1:707029/179‑1 (MQ=255)
                                                                                            cgAAACTTTGATTTCAAAAGTACTCTGCGCGCCAACCTGCAACACTGGCACCCGCAACACGGCAAGTTGTATACCGAATCCCCCCGCTTTAACAGCCGCATTAAACGCCAAAGCGAACAATGGCAACTGGTCTGACTGGTTGATCAAAGCGGATCGATGGTCGATTCGGTGATCCACTCTGCGGTGATGGCGGCCTGTTTGTGGCAGTTACCCGGCATTCGTACCCATCTGGTGGc                                                                                                                                 >  1:448675/1‑236 (MQ=255)
                                                                                                aCTTTGATTTCAAAAGTACTCTGCGCGCCAACCTGCAACACTGGCACCCGCAACACGGCAAGTTGTATATCGAATCCCCCCGCTTTAACAGCCGCATTAAACGCCAAAGCGAACAATGGCAACTGGTCTGACTGGTTGATCAAAGCGGATCGATGGTCGATTCGGTGATCCACTCTGCGGTGATGGCGGCCTGTTTGTGGCAGTTACCCGGCATTCGTACCCATCTGGTGGCGtt                                                                                                                              <  2:212279/235‑1 (MQ=255)
                                                                                                                     tGCGCGCCAACCTGCAACACTGGCACCCGCAACACGGCAAGTTGTATATCGAATCCCCCCGCTTTAACAGCCGCATTAAACGCCAAAGCGAACAATGGCAACTGGTCTGACTGGTTGATCAAAGCGGATCGATGGTCGATTCg                                                                                                                                                                                                     >  2:465667/1‑143 (MQ=255)
                                                                                                                     tGCGCGCCAACCTGCAACACTGGCACCCGCAACACGGCAAGTTGTATATCGAATCCCCCCGCTTTAACAGCCGCATTAAACGCCAAAGCGAACAATGGCAACTGGTCTGACTGGTTGATCAAAGCGGATCGATGGTCGATTCg                                                                                                                                                                                                     <  1:465667/143‑1 (MQ=255)
                                                                                                                               ccTGCAACCCTGGCCCCCGCAACCCGGCAAGTTGTATATCGAATCCCCCCGCTTTAACAGCCGCATTAAACGCCAAAGCGAACAATGGCAACTGGTCTGACTGGTTGATCAAAGCGGATCGATGGTCGATTCGGTGATCCACTCTGCGGTGATGGCGGCCTGTTTGTGGCAg                                                                                                                                                              <  2:120584/172‑1 (MQ=255)
                                                                                                                                      cacTGGCACCCGCAACACGGCAAGTTGTATATCGAATCCCCCCGCTTTAACAGCCGCATTAAACGCCAAAGCGAACAATGGCAACTGGTCTGACTGGTTGATCAAAGCGGATCGATGGTCGATTCGGTGATCCACTCTGCGGTGATGGCGGCCTGTTTGTGGCAGTTACCCg                                                                                                                                                       >  1:120584/1‑172 (MQ=255)
                                                                                                                                        cTGGCACCCGCAACACGGCAAGTTGTATATCGAATCCCCCCGCTTTAACAGCCGCATTAAACGCCAAAGCGAACAATGGCAACTGGTCTGACTGGTTGATCAAAGCGGATCGATGGTCGATTCGGTGATCCACTCTGCGGTGATGGCGGCCTGTTTGTGGCAGTTACCCGGCATTCGTACCCATCTGGTGGCGTTTGACACAAGCGTCGTTGATCTCa                                                                                                       >  1:555238/1‑218 (MQ=255)
                                                                                                                                           gcccccGCAACCCGGCAAGTTGTATATCGAATCCCCCCGCTTTAACAGCCGCATTAAACGCCAAAGCGAACAATGGCAACTGGTCTGACTGGTTGATCAAAGCGGATCGATGGTCGATTCGGTGATCCACTCTGCGGTGATGGCGGCCTGTTTGTGGCAGTTACCCGGCATTCGTACCCATCTGGTGGCGTTTGACACAAGCGTCGTTGATCTCACGGCAGACGTTGCCGATCCgg                                                                                  <  1:431745/233‑1 (MQ=255)
                                                                                                                                                cGCAACACGGCAAGTTGTATATCGAATCCCCCCGCTTTAACAGCCGCATTAAACGCCAAAGCGAACAATGGCAACTGGTCTGACTGGTTGATCAAAGCGGATCGATGGTCGATTCGGTGATCCACTCTGCGGTGATGGCGGCCTGTTTGTGGCAGTTACCCGGCATTCGTACCCATCTGGTGGCGTTTGACACAAGCGTCGTTGATCTCACGGCAGACGTTTCCGATCCGGTagag                                                                             >  1:718549/1‑236 (MQ=255)
                                                                                                                                                cGCAACACGGCAAGTTGTATATCGAATCCCCCCGCTTTAACAGCCGCATTAAACGCCAAAGCGAACAATGGCAACTGGTCTGACTGGTTGATCAAAGCGGATCGATGGTCGATTCGGTGATCCACTCTGCGGTGATGGCGGCCTGTTTGTGGCAGTTACCCGGCATTCGTACCCATCTGGTGGCGTGTGACTCAAACGTCGGTGATCTCACTGCAGACGGTGCCGATCCGGTagag                                                                             >  2:690296/1‑236 (MQ=255)
                                                                                                                                                 gCAACACGGCAAGTTGTATATCGAATCCCCCCGCTTTAACAGCCGCATTAAACGCCAAAGCGAACAATGGCAACTGGTCTGACTGGTTGATCAAAGCGGATCGATGGTCGATTCGGTGATCCACTCTGCGGTGATGGCGGCCTGTTTGTGGCAGTTACCCGGCATTCGTACCCATCTGGTGGCGTTTGACACAAGCGTCGTTGATCTCACGGCAGACGTTGCCGATTCGGTAGAGt                                                                            >  1:399869/1‑236 (MQ=255)
                                                                                                                                                  cAACACGGCAAGTTGTATATCGAATCCCCCCGCTTTAACAGCCGCATTAAACGCCAAAGCGAACAATGGCAACTGGTCTGACTGGTTGATCAAAGCGGATCGATGGTCGATTCGGTGATCCACTCTGCGGTGATGGCGGCCTGTTTGTGGCAGTTACCCGGCATTCGTACCCATCTGGTGGCGTTTGACGCAAGCGTCGTTGATCTCACGGCAGACGTGTGCGATCCGGTagagcg                                                                           >  2:550398/1‑234 (MQ=255)
                                                                                                                                                           aaGTTGTATATCGAATCCCCCCGCTTTAACAGCCGCATTAAACGCCAAAGCGAACAATGGCAACTGGTCTGACTGGTTGATCAAAGCGGATCGATGGTCGATTCGGTGATCCACTCTGCGGTGATGGCGGCCTGTTTGTGGCAGTTACCCGGCATTCGTACCCATCTGGTGGCGTTTGACACAAGCGTCGTTGATCTCACGGCAGACGTTGCCGATCCGGTa                                                                                <  1:533134/222‑1 (MQ=255)
                                                                                                                                                               tGTATATCGAATCCCCCCGCTTTAACAGCCGCATTAAACGCCAAAGCGAACAATGGCAACTGGTCTGAATGGTTGATCAAAGCGGATCGATGGTCGATTCGGTGATCCACTCTGCGGTGATGGCGGCCTGTTTGTGGCAGTTACCCGGCATTCGTACCCATCTGGTGGCGTTTGACACAAGCGTCGTTGATCTCACGGCAGACGTTGCCGATCCGGTAGAGTTATTAATGAAAGTa                                                              <  1:349225/236‑1 (MQ=255)
                                                                                                                                                               tGTATATCGAACCCACCCGCTTTAACAGCCGCATTAAACGCCAAAGCGAACAAGGGCAACTGGTCTGACTGGTTGATCAAAGCGGATCGATGGTCGATTCGGTGATCCACTCTGCGGTGATGGCGGCCTGTTTGTGGCAGTTACCCGGCGTTCGTACCCATCTGGTGGCGTTTGACACAAGCGTCGTTGATCTCACGGCAGACGTTGCCGATCCGGTAGAGTTATTAATGAAAGTa                                                              <  1:601421/236‑1 (MQ=255)
                                                                                                                                                                 tataTCGAATCCCCCCGCTTTAACATCCGCATTAAACGCCAAAGCGAACAATGGCAACTGGTCTGACTGGTTGATCAAAGCGGATCGATGGTCGATTCGGTGATCCACTCTGCGGTGATGGCGGCCTGTTTGTGGCAGTTACCCGGCATTCGTACCCATCTGGTGGCGTTTGACACAAGCGTCGTTGATCTCACGGCAGACGTTGCCGATCCGGTAGAGTTATTAATGAAAg                                                                <  1:718269/232‑1 (MQ=255)
                                                                                                                                                                  atatCGAATCCCCCCGCTTTAACAGCCGCATTAAACGCCAAAGCGAACAATGGCAACTGGTCTGACTGGTTGATCAAAGCGGATCGATGGTCGATTCGGTGATCCACTCTGCGGTGATGGCGGCCTGTTTGTGGCAGTTACCCGGCATTCGTACCCATCTGGTGGCGTTTGACACAAGCGTCGTTGATCTCACGGAAGACGGTGCCGGTCCGGGagagctgt                                                                         >  2:533134/1‑218 (MQ=255)
                                                                                                                                                                   tatCGAATCCCCCCGCTTTAACAGCCGCATTAAACGCCAAAGCGAACAATGGCAACTGGTCTGACTGGTTGATCAAAGCGGATCGATGGTCGATTCGGTGATCCACTCTGCGGTGATGGCGGCCTGTTTGTGGCAGTTACCCGGCATTCGTACCCATCTGGTGGCGTTTGACACAAGCGTCGTTGATCTCACGGCAGACGTTGCCGATCCGGTAGAGTTATTAATGAAAGTACAGt                                                          >  2:108679/1‑236 (MQ=255)
                                                                                                                                                                        aaTCCCCCCGCTTTAACAGCCGCATTAAACGCCAAAGCGAACAATGGCAACTGGTCTGACTGGTTGATCAAAGCGGATCGATGGTCGATTCGGTGATCCACTCTGCGGTGATGGCGGCCTGTTTGTGGCAGTTACCCGGCATTCGTACCCATCTGGTGGCGTTTGACACAAGCGTCGTTGACCTCACTGCAGTCGTTTCAGAGCCGGTAGAGGTATATATGTAAGTACAGtt                                                         >  2:718269/1‑232 (MQ=255)
                                                                                                                                                                        aaTCCCCCCGCTTTAACAGCCGCATTAAACGCCAAAGCGAACAATGGCAACTGGTCTGACTGGTTGATCAAAGCGGATCGATGGTCGATTCGGTGATCCACTCTGCGGTGATGGCGGCCTGTTTGTGGCAGTTACCCGGCATTCGTACCCATCTGGTGGCGTTTGACACAAGAGTCGTTGATCTCACGGCAGACGTTGCCGATCCGGTAGAGTTATTAATGAAAGTaacgttcgcg                                                     >  2:728449/1‑227 (MQ=255)
                                                                                                                                                                          tCCCCCCGCTTTAACAGCCGCATTAAACGCCAAAGCGAACAATGGCAACTGGTCTGACTGGTTGATCAAAGCGGATCGATGGTCGATTCGGTGATCCACTCTGCGGTGATGGCGGCCTGTTTGTGGCAGTTACCCGGCATTCGTACCCATCTGGTGGCGTTTGACACAAGCGTCGTTGATCTCACGGCAGACGTTGCCGATCCGGTAGAGTTATTAATGAAAg                                                                >  2:693657/1‑223 (MQ=255)
                                                                                                                                                                          tCCCCCCGCTTTAACAGCCGCATTAAACGCCAAAGCGAACAATGGCAACTGGTCTGACTGGTTGATCAAAGCGGATCGATGGTCGATTCGGTGATCCACTCTGCGGTGATGGCGGCCTGTTTGTGGCAGTTACCCGGCATTCGTACCCATCTGGTGGCGTTTGACACAAGCGTCGTTGATCTCACGGCAGACGTTGCCGATCCGGTAGAGTTATTAATGAAAg                                                                <  1:693657/223‑1 (MQ=255)
                                                                                                                                                                             ccccGTTTGAACAGCCGCATTAAACGCCAATGCGAACAATGGCAACTGGTCTGACTGGTTGATCAAAGCGGATCGATGGTCGATTCGGTGATCCACTCTGCGGTGATGGCGGCCTGTTTGTGGCAGTTACCCGGCATTCGTACCCATCTGGTGGCGTTTGACACAAGCGTCGTTGATCTCACGGCAGACGTTGCCGATCCGGTAGAGTTATTAATGAAAGTACAGTTGGGcggcgg                                                <  2:399869/236‑1 (MQ=255)
                                                                                                                                                                                             gCATTAAACGCCAAAGCGAACAATGGCAACTGGTCTGACTGGTTGATCAAAGCGGATCGATGGTCGATTCGGTGATCCACTCTGCGGTGATGGCGGCCTGTTTGTGGCAGTTACCCGGCATTCGTACCCATCTGGTGGCGTTTGACACAAGCGGCGTTGATCTCACGGCAGACGGTGCCGATCCGGTAGAGGTATTAATGAAAGGACAGATGGGCGGCCGGTCCAATATTGCCAGt                                >  2:226022/1‑236 (MQ=255)
                                                                                                                                                                                               aTTAAACGCCAAAGCGAACAATGGCAACTGGTCTGACTGGTTGATCAAAGCGGATCGATGGTCGATTCGGTGATCCACTCTGCGGTGATGGCGGCCTGTTTGTGGCAGTTACCCGGCATTCGTACCCATCTGGTGGCGTTTGACACAAGCGTCGTTGATCTCACGGCAGACGTTGCCGATCCGGTAGAGTTATTAATGAAAGTACAGTTGGGCGGCGGGACCAATATCGCCa                                  <  1:703982/232‑1 (MQ=255)
                                                                                                                                                                                                  aaaCGCCAAAGCGAACAATGGCAACTGGTCTGACTGGTTGATCAAAGCGGATCGATGGTCGATTCGGTGATCCACTCTGCGGTGATGGCGGCCTGTTTGTGGCAGTTACCCGGCATTCGTACCCATCTGGTGGCGTTTGACACAAGCGTCGTTGATCTCACGGCAGACGTTGCCGATCCGGTAGAGTTATTAATGAAAGTACAGTTGGGCGGCGGGACCAATATCGCCAGTGCCg                            <  1:550398/235‑1 (MQ=255)
                                                                                                                                                                                                           aGCGAACAATGGCAACTGGTCTGACTGGTTGATCAAAGCGGATCGATGGTCGATTCGGTg                                                                                                                                                                                                  >  1:315635/1‑60 (MQ=255)
                                                                                                                                                                                                           aGCGAACAATGGCAACTGGTCTGACTGGTTGATCAAAGCGGATCGATGGTCGATTCGGTg                                                                                                                                                                                                  <  2:315635/60‑1 (MQ=255)
                                                                                                                                                                                                                aCAATGGCAACTGGTCTGACTGGTTGGTCAAAGCGGATCGATGGTCGATTCGGTGATCCACTCTGCGGTGATGGCGGCCTGTTTGTGGCAGTTACCCGGCATTCGTACCCATCTGGTGGCGTTTGACACAAGCGTCGTTGATCTCACGGCAGACGTTGCCGATCCGGTAGAGTTATTAATGAAAGTACAGTTGGGCGGCGGGACCAATATCGCCAGTGCCGTGGAGTATGGTCg               >  1:641136/1‑234 (MQ=255)
                                                                                                                                                                                                                  attGGCAAGTGGACTTGCTGGTGAGTCAAAGCGGTTCGATGGTCGATTCGTTGATCCAGTCTGCGGTGATGGTGGCCTTTTTGTGGCAGTTCCCCGGCATTCGTACCCATCTGGTGGCGTTTGACGCAAGCGTCGTTGATCTCACGGCGGACGTTGACGATCCGGTAGAGTTATTAATGAAAGTACAGTTGGGCGGCGGGACCAATATCGCCAGTGCCGTGGAGTATGGTCGGCa            <  2:376007/233‑1 (MQ=255)
                                                                                                                                                                                                                            ggTCTGACTGGTTGATCAAAGCGGATCGATGGTCGATTCGGTGATCCACTCTGCGGTGATGGCGGCCTGTTTGTGGCAGTTACCCGGCATTCGTACCCATCTGGTGGCGTTTGACACAAGCGTCGTTGATCTCACGGCAGACGGTGCCGATCCGGTAGAGTTATTAATGAAAGTACAGTTGGGCGGCGGGACCAATATCGCCAGGGCCCGGGCGATTGGGTGGCAACTTTTTGaa  >  1:20764/1‑235 (MQ=255)
                                                                                                                                                                                                                             gTCTGACTGGTTGATCAAAGCGGATCGATGGTCGATTCGGTGATCCACTCTGCGGTGATGGCGGCCTGTTTGTGGCAGTTACCCGGCATTCGTACCCATCTGGTGGCGTTTGACACAAGCGTCGTTGATCTCACGGCAGACGTTGCCGATCCGGTAGAGTTATTAATGAAAGTACAGTTGGGCGGCGGGACCAATATCGCCAGTGCCGTGGAGTa                     <  1:338990/215‑1 (MQ=255)
                                                                                                                                                                                                                             gTCTGACTGGTTGATCAAAGCGGATCGATGGTCGATTCGGTGATCCACTCTGCGGTGATGGCGGCCTGTTTGTGGCAGTTACCCGGCATTCGTACCCATCTGGTGGCGTTTGACACAAGCGTCGTTGATCTCACGGCAGACGTTGCCGATCCGGTAGAGTTATTAATGAAAGTACAGTTGGGCGGCGGGACCAATATCGCCAGTGCCGTGGAGTa                     >  2:338990/1‑215 (MQ=255)
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                     
TTGAAGAAATTATGGCTCGACTGGCAAAGGAAGTTCGTCAGGCTTTTTCTGGTGTCCGCGATCGCCGTCGCCGTTCATTTATTCCACTGGCGCGAAACTTTGATTTCAAAAGTACTCTGCGCGCCAACCTGCAACACTGGCACCCGCAACACGGCAAGTTGTATATCGAATCCCCCCGCTTTAACAGCCGCATTAAACGCCAAAGCGAACAATGGCAACTGGTCTTACTGGTTGATCAAAGCGGATCGATGGTCGATTCGGTGATCCACTCTGCGGTGATGGCGGCCTGTTTGTGGCAGTTACCCGGCATTCGTACCCATCTGGTGGCGTTTGACACAAGCGTCGTTGATCTCACGGCAGACGTTGCCGATCCGGTAGAGTTATTAATGAAAGTACAGTTGGGCGGCGGGACCAATATCGCCAGTGCCGTGGAGTATGGTCGGCAACTTATTGAA  >  NC_000913/2208397‑2208851

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