Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,696,909 G→A G151E (GGA→GAA)  bcsE → c‑di‑GMP‑binding protein BcsE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,696,9090GA95.3% 187.1 / ‑0.9 64G151E (GGA→GAA) bcsEc‑di‑GMP‑binding protein BcsE
Reads supporting (aligned to +/- strand):  ref base G (1/2);  new base A (34/27);  total (35/29)
Fisher's exact test for biased strand distribution p-value = 5.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

TAGATGATTCTCAAGGGCCGGAAAAAATAAAATTATTTTCAATGCTAAATCATGAAAAAGGTCTATACTATTTGACCCGTGATTTGCAGTGTTCTATTGATCCCCATAATTACCTTTTTATTCTTGTTTGCGCAAATAACGCATGGCAAAACATTCCTGCCGAGCGGCTTCGCTCATGGTTGGATAAAATGAATAAATGGAGCAGGTTAAACCATTGTTCGCTTTTGGTAATTAATCCCGGAAATAATAACGATAAACAATTTTCATTGTTGCTTGAGGAATACCGTTCACTTTTTGGTCTTGCCAGTTTGCGTTTTCAGGGTGACCAACATTTGCTGGATATTGCCTTCTGGTGCAACGAAAAAGGGGTCAGCGCCCGTCAGCAGCTTAGCGTTCAGCAACAAAATGGTATCTGGACATTAGTTCAAAGCGAAGAGGCGGAGATCCAACCACGCAGCGACGAAAAACGCATTCTGAGT  >  NC_000913/3696669‑3697147
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                              
tAGATGATTCTCAAGGGCCGGAAAAAATAAAATTATTTTCAATGCTAAATCATGAAAAAGGTCTATACTATTTGACCCGTGATTTGCAGTGTTCTATTGATCCCCATAATTACCTTTTTATTCTTGTTTGCGCAAATAACGCATGGCAAAACATTCCTGCCGAGCGGCTTCGCTCATGGTTGGATAAAATGAATAAATGGAGCAGGTTAAACCATTGTTCGCTTTTGGTAATTAATCCCGAAAAtaataa                                                                                                                                                                                                                                       <  1:299173/250‑1 (MQ=255)
 aGATGATTCTCAAGGGCCGGAAAAAATAAAATTATTTTCAATGCTAAATCATGAAAAAGGTCTATACTATTTGACCCGTGATTTGCAGTGTTCTATTGATCCCCATAATTACCTTTTTATTCTTGTTTGCGCAAATAACGCATGGCAAAACATTCCTGCCGAGCGGCTTCGCTCATGGTTGGATAAAATGAATAAATGGAGCAGGTTAAACCATTGTTCGCTTTTGGTAATTAATCCCGAAAATAATAAc                                                                                                                                                                                                                                      <  1:132604/250‑1 (MQ=255)
 aGATGATTCTCAAGGGCCGGAAAAAATAAAATTATTTTCAATGCTAAATCATGAAAAAGGTCTATACTATTTGACCCGTGATTTGCAGTGTTCTATTGATCCCCATAATTACCTTTTTATTCTTGTTTGCGCAAATAACGCATGGCAAAACATTCCTGCCGAGCGGCTTCGCTCATGGTTGGATAAAATGAATAAATGGAGCAGGTTAAACCATTGTTCGCTTTTGGTAATTAATCCCGAAAATAATAAc                                                                                                                                                                                                                                      <  1:376499/250‑1 (MQ=255)
           cAAGGGCCGGAAAAAATAAAATTATTTTCAATGCTAAATCATGAAAAATGTCTATACTATATGACCCGTGATTTGCAGTGTTCTATTGATCCCCATAATTACCTTTTTATTCTTGTTTGCGCAAATAACGCATGGCAAAACATTCCTGCCGAGCGGCTTCGCTCATGGTTGGATAAAATGAATAAATGGAGCAGGTTAAACCATTGTTCGCTTTTGGTAATTAATCCCGAAAATAATAACGATAAACAAtt                                                                                                                                                                                                                           <  2:35373/251‑1 (MQ=255)
           cAAGGGCCGGAAAAAATAAAATTATTTTCAATGCTAAATCATGAAAAAGGTCTATACTATTTGACCCGTGATTTGCAGTGTTCTATTGATCCCCATAATTACCTTTTTATTCTTGTTTGCGCAAATAACGCATGGCAAAACATTCCTGCCGAGCGGCTTCGCTCATGGTTGGATAAAATGAATAAATGGAGCAGGTTAAACCATTGTTCGCTTTTGGTAATTAATCCCGAAAATAATAACGATAAACAAtt                                                                                                                                                                                                                           >  2:265629/1‑251 (MQ=255)
             aGGGCCGGAACAAATAAAATTATTTTCAATGCTAAATCATGAAAAGGGTCTATACTATTTGACCCGTGATTTGCAGTGTTCTATTAATCCCCATAATTACCTTTTTATTCTTGTTTGCGCAAATAACGCATGGCAAAACATTCCTGCCGAGCGGCTTCGATCATGGTTGGATAAAATGAATAAATGGAGCAGGTTACACCATTGTTCGCTTTTGGTAATTAATCCCGAAAATAATAACGATaacacatttt                                                                                                                                                                                                                         >  2:347581/1‑243 (MQ=255)
              gggCCGGAAAAAATAAAATTATTTTCAATGCTAAATCATGAAAAAGGTCTATACTATTTGACCCGTGATTTGCAGTGTTCTATTGATCCCCATAATTACCTTTTTATTCTTGTTTGCGCAAATAACGCATGGCAAAACATTCCTGCCGAGCGGCTTCGCTCATGGTTGGATAAAATGAATAAATGGAGCAGGTTAAACCATTGTTCGCTTTTGGTAATTAATCCCGAAAATAATAACGATAAACAATTTTc                                                                                                                                                                                                                        >  1:4599/1‑251 (MQ=255)
               ggCCGGAAAAAATAAAATTATTTTCAATGCTAAATCATGAAAAAGGTCTATACTATTTGACCCGTGATTTGCAGTGTTCTATTGATCCCCATAATTACCTTTTTATTCTTGTTTGCGCAAATAACGCATGGCAAAACATTCCTGCCGAGCGGCTTCGCTCATGGTTGGATAAAATGAATAAATGGAGCAGGTTAAACCATTGTTCGCTTTTGGTAATTAATCCCGAAAATAATAACGATAAACAATTTTc                                                                                                                                                                                                                        >  1:420932/1‑250 (MQ=255)
                    gAAAAAATAAAATTATTTTCAATGCTAAATCATGAAAAAGGTCTATACTATTTGACCCGTGATTTGCAGTGTTCTATTGATCCCCATAATTACCTTTTTATTCTTGTTTGCGCAAATAACGCATGGCAAAACATTCCTGCCGAGCGGCTTCGCTCATGGTTGGATAAAATGAATAAATGGAGCAGGTTAAACCATTGTTCGCTTTTGGTAATTAATCCCGAAAATAATAACGATAAACAATTTTCAttgtt                                                                                                                                                                                                                  <  2:352709/251‑1 (MQ=255)
                     aaaaaaTAAAATTATTTTCAATGCTAAATCATGAAAAAGGTCTATACTATTTGACCCGTGATTTGCAGTGTTCTATTGATCCCCATAATTACCTTTTTATTCTTGTTTGCGCAAATAACGCATGGCAAAACATTCCTGCCGAGCGGCTTCGCTCATGGTTGGATAAAATGAATAAATGGAGCAGGTTAAACCATTGTTCGCTTTTGGTAATTAATCCCGAAAATAATAACGATAAACAATTTTCAttgttg                                                                                                                                                                                                                 >  1:64122/1‑251 (MQ=255)
                               attattTTCAATGCTAAATCATGAAAAAGGTCTATACGATTTGACCCGGGATTTGCAGTGTTATATTGATCCACATAATTACCTTTTTATTCTTGTTTGCGCAAATAACGCATGGCAAAACATTCCTGCCGAGCGGCTTCGCTCATGGTTGGATAAAATGAATAAATGGAGCAGGTTAAACCATTGTTCGCTTTTGGTAATTAATCCCGAAAATAATAACGATAAACAATTTTCATTGTTGCTTGAGGAAt                                                                                                                                                                                                       <  2:64122/251‑1 (MQ=255)
                                      tcaattcgaaaTCATGAAAATGCTCTATACTATTTGACCCGTTAGTTGCAGTGAGCTATAGATCCCCATAATTACCTTTTTATTCTTGTTTGCGCAAATAACACATGGCAAAACCTTCCTGCCGAGCTGCTTCGCTCATGGTTGGATAAAATTAATAAATGGAGCAGGTTAAACCATTGTGCGCTTTTGGTAATGAATCCCGAAAATAATAACGAGAAACAAGCGTCATTTTTGCTTGAGTAATAccatt                                                                                                                                                                                                 <  2:107982/242‑4 (MQ=255)
                                       aaagtcTAAACCATGAAAAAGGTCTTTACTATTTGACCCGTGATTTGCAGTGTTCTATTGAGCCCCATAATTACCTTTTTATTCTGGTTTGCGAAAATAACGCATGGCAAAACATTCCTGCCGAGCGGCTTCGCTCATGGTTGGATAAAATGAATGAATGGAGCAGGTTAAACCATTGTTCGCTTTTGGTAATTAAGCCCGGAAATAATAACGATAAACAATTTTCATTGTGGCTTGAGGAATACCGTTCa                                                                                                                                                                                               <  2:767811/246‑1 (MQ=255)
                                                aTCATGAAAAAGGTCTATACTATTTGACCCGTGATTTGCAGTGTTCTATTGATCCCCATAATTACCTTTTTATTCTTGTTTGCGCAAATAACGCATGGCAAAACATTCCTGCCGAGCGGCTTCGCTCATGGTTGGATAAAATGAATAAATGGAGCAGGTTAAACCATTGTTCGCTTTTGGTAATTAATCCCGAAAATAATAACGATAAACAATTTTCATTGTTGCTTGAGGAATACCGTTCACTTTTTGGt                                                                                                                                                                                      <  2:420932/251‑1 (MQ=255)
                                                    tGAAAAAGGTCTATACTATTTGACCCGTGATTTGCAGTGTTCTATTGATCCCCATAATTACCTTTTTATTCTTGTTTGCGCAAATAACGCATGGCAAAACATTCCTGCCGAGCGGCTTCGCTCATGGTTGGATAAAATGAATAAATGGAGCAGGTTAAACCATTGTTCGCTTTTGGTAATTAATCCCGAAAATAATAACGATAAACAATTTTCATTGTTGCTTGAGGAATACCGTTCACTTTTTGGTCTTg                                                                                                                                                                                  >  2:294086/1‑251 (MQ=255)
                                                       aaaaGGTCTAATCTATTTGTCCCGTGATTTGCAGTGTTCTATTGATCCCCATAATTACCTTTTTATTCTTGTTTGCGCAAATAACGCATGGCAAAACATTCCTGCCGAGCGGCTTCGCTCATGGTTGGATAAAATGAATAAATGGAGCAGGTTAAACCATTGTTCGCTTTTGGTAATTAATCCCGAAAATAATAACGATAAACAATTTTCATTGTTGCTTGAGGAATACCGTTCACTTTTTGGTCTTGCCa                                                                                                                                                                               <  2:4307/251‑1 (MQ=255)
                                                          aGGTCTATACTATTTGACCCGTGATTTGCAGTGATCTATTGATCCACCTAATTACCTTTTTATTCTTGTTTGCGCAAATAACGCATGGCAAAACATTCCTGCCGAGCGGCTTCGCTCATGGTTGGATAAAATGAATAAATGGAGCAGGTTAAACCATTGTTCGCTTTTGGTAATTAATCCCGAAAATAATAACGATAAACAATTTTCATTGTTGCTTGAGGAATACCGTTCACTTTTTGGTCTTGCCCGtt                                                                                                                                                                            >  1:505265/1‑251 (MQ=255)
                                                                    tATTTGACCCGTGATTTGCAGTGTTCTATTGATCCCCATAATTACCTTTTTATTCTTGTTTGCGCAAATAACGCATGGCAAAACATTCCTGCCGAGCGGCTTCGCTCATGGTTGGATAAAATGAATAAATGGAGCAGGTTAAACCATTGTTCGCTTTTGGTAATTAATCCCGAAAATAATAACGATAAACAATTTTCATTGTTGCTTGAGGAATACCGTTCACTTTTTGGTCTTGCCAGTTTGCGTTTTc                                                                                                                                                                   >  1:618426/1‑250 (MQ=255)
                                                                     agttgatccttGATTTGAAGTTTTCTATTGATCCCCATAATTACCTTTTTATTCTTGTTTGCACAAATAACTCATGGCAAAACATTCCTGCCGAGCGGCTTCGTTCATGGTTGGATAAAATGAATAAATGGAGCAGGTTAAACCATTGTTCGCTTTTGGTTATTAATCCCGAAAATAATAACGATAAACAATTGTCATTGTTGCTTGAGGAATTCCGTTCACTTTTTGGGCTTGCCAGTTTGCGTTTTCa                                                                                                                                                                  <  2:456818/240‑1 (MQ=255)
                                                                     aTTTGACCCGTGATTTGCAGTGTTCTATTGATCCCCATAATTACCTTTTTATTCTTGTTTGCGCAAATAACGCATGGCAAAACATTCCTGCCGAGCGGCTTCGCTCATGGTTGGATAAAATGAATAAATGGAGCAGGTTAAACCATTGTTCGCTTTTGGTAATTAATCCCGAAAATAATAACGATAAACAATTTTCATTGTTGCTTGAGGAATACCGTTCACTTTTTGGTCTTGCCAGTTTGCGTTTTCa                                                                                                                                                                  <  1:294086/250‑1 (MQ=255)
                                                                      tttGACCCGTGATTTGCAGTGTTCTATTGATCCCCATAATTACCTTTTTATTCTTGTTTGCGCAAATAACGCATGGCAAAACATTCCTGCCGAGCGGCTTCGCTCATGGTTGGATAAAATGAATAAATGGAGCAGGTTAAACCATTGTTCGCTTTTGGTAATTAATCCCGAAAATAATAACGATAAACAATTTTCATTGTTGCTTGAGGAATACCGTTCACTTTTTGGTCTTGCCAGTTTGCGTTTTCAgg                                                                                                                                                                >  1:529130/1‑251 (MQ=255)
                                                                       ttGACCCGTGATTTGCAGGGGTCTATTGATCCCCATAATTACCTTGTTATTCTTGTTTTCGCAAATAACGCATGGCAAAACATTCCTGCCGAGCGGCTTCGCTCATGGTTGGATAAAATGAATAAATGGAGCAGGTTAAACCATTGTTCGCTTTTGGTAATTAATCCCGAAAATAATAACGATAAACAATTTTCATTGTTGCTTGAGGAATACCGTTCACTTTTTGGTCTTGCCAGTTTGCGTTTTCAgg                                                                                                                                                                <  2:51062/250‑1 (MQ=255)
                                                                                       aGTGTTCTATTGATCCCCATAATTACCTTTTTATTCTTGTTTGCGCAAATAACGCATGGCAAAACATTCCTGCCGAGCGGCTTCGCTCATGGTTGGATAAAATGAATAAATGGAGCAGGTTAAACCATTGTTCGCTTTTGGTAATTAATCCCGAAAATAATAACGATAAACAATTTTCATTGTTGCTTGAGGAATACCGTTCACTTTTTGGTCTTGCCAGTTTGCGTTTTCAGGGTGACCAACATGTGCTg                                                                                                                                               <  1:717819/251‑1 (MQ=255)
                                                                                             cTATTGATCCCCATAATTACCTTATTATTCTTGTTTGCGCAAATAACGCATGGCAAAACATTCCTGCCGAGCGGCTTCGCTCATGGTTGGATAAAATGAATAAATGGAGCAGGTTAAACCATTGTTCGCTTTTGGTAATTAATCCCGAAAATAATAACGATAAACAATTTTCATTGTTGCTTGAGGAATACCGTTCACTTTTTGGTCTTGCCAGTTTGCGTTTTCAGTGTGACCAACATTTGCTGGATAtt                                                                                                                                         >  2:800917/1‑251 (MQ=255)
                                                                                              tATTGATCCCCATAATTACCTTTTTATTCTTGTTTGCGCAAATAACGCATGGCAAAACATTCCTGCCGAGCGGCTTCGCTCATGGTTGGATAAAATGAATAAATGGAGCAGGTTAAACCATTGTTCGCTTTTGGTAATTAATCCCGGAaat                                                                                                                                                                                                                                            <  1:19126/151‑1 (MQ=255)
                                                                                              tATTGATCCCCATAATTACCTTTTTATTCTTGTTTGCGCAAATAACGCATGGCAAAACATTCCTGCCGAGCGGCTTCGCTCATGGTTGGATAAAATGAATAAATGGAGCAGGTTAAACCATTGTTCGCTTTTGGTAATTAATCCCGGAaat                                                                                                                                                                                                                                            >  2:19126/1‑151 (MQ=255)
                                                                                                ttgttCCCCATAATTTCCTTATTATTCTTGTTTGCGCAAATAACGCATGGCAAAACATTCCTGCCGAGCGGCTTCGCTCATGGTTGGATAAAATTAATAAATGGAGCAGGTTAATCCATTTTTCGCTTTTGGTATTGAATCTCGAAAATAATAACGATAAACAATTTTTATTGTTGCTTGGTGAAGACCGATCACTTTTTGGTCTGGCCAGTTTGCGTTTTCAGGGTGTCGAACATTTGCTGGATATTGc                                                                                                                                       <  2:856381/246‑1 (MQ=255)
                                                                                                  gATCCCCATTATTACCTTTTTATTCGTGGTTTCGCATATAACGCATGTCCAAACAATCCTGCCGAGCTGCTTCGCTCATTGTTGGATAAAATGAATAAATGGAGCAGGATAAACCATTGGTCGCTTTTGGTAAGTAATCCCGAAAATAATAACGATAAACAATTTTCATTGTTGCTTGAGGAATACCGTTCATTTTTTGGTCTTGCCATTTTGCGTTTTCAGGGTGACCACCATTTGCTGGATATTGCCtt                                                                                                                                    <  2:29090/251‑1 (MQ=255)
                                                                                                  gATCCCCATAATTACCTTTTTATTCTTGTTTGCGCAAATAACGCATGGCAAAACATTCCTGCCGAGCGGCTTCGCTCATGGTTGGATAAAATGAATAAATGGAGCAGGTTAAACCATTGTTCGCTTTTGGTAATTAATCCCGAAAATAATAACGATAAACAATTTTCATTGTTGCTTGAGGAATACCGTTCACTTTTTGGTCTTGCCAGTTTGCGTTTTCAGGGTGACCAACATTTGCTGGATATTGCCtt                                                                                                                                    >  2:733570/1‑251 (MQ=255)
                                                                                                    tCCCCATAATTACCTTTTTATTCTTGTTTGCGCCAATAACGCATGACAAAACATTCCTGCCTAGCGGCTTCGCTCATGGTTGGATAACATGAATAACTGGAGCAGGTTAAACCATTGTTCGCTTTTGGTAATTAATCCCGAAAATAATAACGATCAACACTTTTAATTGTTGCTTGCGGATTACCGTTCACTTTTTGATCTTGCCAGTTTGCGTTTTCAGTGTGAACACGATTTGCAGGATATTGCCTTCt                                                                                                                                  >  2:812828/1‑251 (MQ=255)
                                                                                                     ccccATAATTACCTTATTATTCTTGTTTGCGCAAATAACGCATGGCAAAACATTCCTGCCGAGCGGCTTCGCTCATGGTTGGATAAAATGAATAAATGGAGCAGGTTAAACCATTGTTCGCTTTTGGTAATTAATCCCGAAAATAATACCGCTAAACAATTTTCATTGTTGCTTGAGGAATACCGTTCACTTTTTGGTCTTGCCAGTTTGCGTTTTCAGGGTGACCAACATTTGCTGGATATTGCCTTCTg                                                                                                                                 >  2:518591/1‑251 (MQ=255)
                                                                                                           ccttcccTTTTTATTCTTGTTTTCGCAAATAACGCATGGCAAAACATTCCTGCCGAGCGGCTTCGCTCATGGTTGGATAAAATGAATAAATGGAGCAGGTTAAACCATTGTTCGCTTTTGGTAATTAATCCCGAAAATAATAACGATAAACAATTTTCATTGTTGCTTGAGGAATACCGTTCACTTTTTGGTCTTGCCAGTTTGCGTTTTCAGGGTGACCGACATTTGCTGGATATTGCCTTCTGGTGCa                                                                                                                            >  1:778681/6‑250 (MQ=255)
                                                                                                             ttACCTTTTTATTCTTGTTTGCGCAAATAACGCATGGCAAAACATTCCTGCCGAGCGGCTTCGCTCATGGTTGGATAAAATGAATAAATGGAGCAGGTTAAACCATTGTTCGCTTTTGGTAATTAATCCCGAAAATAATAACGATAAACAATTTTCATTGTTGCTTGAGGAATACCGTTCACTTTTTGGTCTTGCCAGTTTGCGTTTTCAGGGTGACCAACATTTGCTGGATATTGCCTTCTGGTGCAACg                                                                                                                         >  1:759692/1‑251 (MQ=255)
                                                                                                                   ttttATTCTTGTTTGCGCAAATAACGCATGGCAAAACATTCCTGCCGAGCGGCTTCGCTCATGGTTGGATAAAATGAATAAATGGAGCAGGTTAAACCATTGTTCGCTTTTGGTAATTAATCCCGAAAATAATAACGATAAACAATTTTCAttgttg                                                                                                                                                                                                                 >  1:750737/1‑157 (MQ=255)
                                                                                                                   ttttATTCTTGTTTGCGCAAATAACGCATGGCAAAACATTCCTGCCGAGCGGCTTCGCTCATGGTTGGATAAAATGAATAAATGGAGCAGGTTAAACCATTGTTCGCTTTTGGTAATTAATCCCGAAAATAATAACGATAAACAATTTTCAttgttg                                                                                                                                                                                                                 <  2:750737/157‑1 (MQ=255)
                                                                                                                      tATTCTTGTTTGCGCAAATAACGCATGGCAAAACATTCCTGCCGAGCGGCTTCGCTCATGGTTGGATAAAATGAATAAATGGAGCAGGTTAAACCATTGTTCGCTTTTGGTAATTAATCCCGAAAATAATAACGATAAACAATTTTCATTGTTGCTTGAGGAATACCGTTCACTTTTTGGTCTTGCCAGTTTGCGTTTTCAGGGTGACCAACATTTGCTGGATATTGCCTTCTGGTGCAACGAAAAAgggg                                                                                                                >  1:352828/1‑251 (MQ=255)
                                                                                                                       aTTCTTGTTTGCGCAAATAACGCATGGCAAAACATTCCTGCCGAGCGGCTTCGCTCATGGTTGGATAAAATGAATAAATGGAGCAGGTTAAACCATTGTTCGCTTTTGGTAATTAATCCCGAAAATAATAACGATAAACAATTTTCATTGTTGCTTGAGGAATACCGTTCACTTTTTGGTCTTGCCAGTTTGCGTTTTCAGGGTGACCAACATTTGCTGGATATTGCCTTCTGGTGCAACGAAAAAGGGGt                                                                                                               <  1:124176/251‑1 (MQ=255)
                                                                                                                          cTTGTTTGCGCAAATAACGCATGGCAAAACATTCCTGCCGAGCGGCTTCGCTCATGGTTGGATAAAATGAATAAATGGAGCAGGTTAAACCATTGTTCGCTTTTGGTAATTAATCCCGAAAATAATAACGATAAACAATTTTCATTGTTGCTTGAGGAATACCGTTCACTTTTTGGTCTTGCCAGTTTGCGTTTTCAGGGTGACCAACATTTGCTGGATTTTGCCTTCTGGTGCACCGAAAAAGGGGTc                                                                                                              >  2:261972/1‑249 (MQ=255)
                                                                                                                                    cAAATAACGCATGGCAAAACATTCCTGACGAGCGGCTTCGCTCATGGTTGGATAAAATGAATAAATGGAGCAGGTTAAACCATTGTTCGCTTTTGGTAATTAATCCCGAAAATAATAACGATAAACAATTTTCATTGTTGCTTGAGGAATACCGTTCCCTTTTTGGTCTTGCCAGTTTGCGTTTTCAGGGTGACCAACATTTGCTGGATATTGCCTTCTGGTGCAACGAAAAAGGGGTCAGCGCCCGta                                                                                                    >  2:858434/1‑248 (MQ=255)
                                                                                                                                      aaTAACGCATGGCAAAACATTCCTGCCGAGCGGCTTCGCTCATGGTTGGATAAAATGAATAAATGGAGCAGGTTAAACCATTGTTCGCTTTTGGTAATTAATCCCGAAAATAATAACGATAAACAATTTTCATTGTTGCTTGAGGAATACCGTTCACTTTTTGGTCTTGCCAGTTTGCGTTTTCAGGGTGACCAACATTTGCTGGATATTGCCTTCTGGTGCAACGAAAAAGGGGTAAGCGCCCGTcagc                                                                                                 >  2:180187/1‑250 (MQ=255)
                                                                                                                                          aCGCATGGCAAAACATTCCTGCCGAGCGGCTTCGCTCATGGTTGGATAAAATGAATAAATGGAGCAGGTTAAACCATTGTTCGCTTTTGGTAATTAATCCCGAAAATAATAACGATAAACAATTTTCATTGTTGCTTGAGGAATACCGTTCACTTTTTGGTCTTGCCAGTTTGCGTTTTCAGGGTGACCAACATTTGCTGGATATTGCCTTCTGGTGCAACGAAAAAGGGGTCAGCGCCCGTCAGCAGCtt                                                                                            <  2:529130/251‑1 (MQ=255)
                                                                                                                                          aCGCATGGCAAAACATTCCTGCCGAGCGGCTTCGCTCATGGTTGGATAAAATGAATAAATGGAGCAGGTTAAACCATTGTTCGCTTTTGGTAATTAATCCCGAAAATAATAACGATAAACAATTTTCATTGTTGCTTGAGGAATACCGTTCACTTTTTGGTCTTGCCAGTTTGCGTTTTCAGGGTGACCAACATTTGCTGGATATTGCCTTCTGGTGCAACGAAAAAGGGGGCAGCGCCCGTCAGCAGCtt                                                                                            <  2:855246/251‑1 (MQ=255)
                                                                                                                                             cATGGCAAAACATTCCTGCCGAGCGGCTTCGCTCATGGTTGGATAAAATGAATAAATGGAGCAGGTTAAACCATTGTTCGCTTTTGGTAATTAATCCCGAAAATAATAACGATAAACAATTTTCATTGTTGCTTGAGGAATACCGTTCACTTTTTGGTCTTGCCAGTTTGCGTTTTCAGGGTGACCAACATTTGCTGGATATTGCCTTCTGGTGCAACGAAAAAGGGGTCAGCGCCCGTCAGCAGCTTAGc                                                                                         <  1:71657/251‑1 (MQ=255)
                                                                                                                                               tGGCAAAACATTCCTGCCGAGCGGCTTCGCTCATGGTTGGATAAAATGAATAAATGGAGCAGGTTAAACCATTGTTCGCTTTTGGTAATTAATCCCGAAAATAATAACGATAAACAATTTTCATTGTTGCTTGAGGAATACCGTTCACTTTTTGGTCTTGCCAGTTTGCGTTTTCAGGGTGACCAACATTTGCTGGATATTGCCTTCTGGTGCAACGAAAAAGGGGTCAGCGCCCGTCAGCAGCTTAGCGt                                                                                       >  1:723200/1‑251 (MQ=255)
                                                                                                                                                               ccGAGCGGCTTCGCTCATGGTTGGATAAAATGAATAAATGGAGCAGGTTAAACCATTGTTCGCTTTTGGTAATTAATCCCGAAAATAATAACGATAAACAATTTTCATTGTTGCTTGAGGAATACCGTTCACTTTTTGGTCTTGCCAGTTTGCGTTTTCAGGGTGACCAACATTTGCTGGATATTGCCTTCTGGTGCAACGAAAAAGGGGTCAGCGCCCGTCAGCAGCTTAGCGTTCAGCAACAAAATgg                                                                        >  2:554089/1‑250 (MQ=255)
                                                                                                                                                                  agcgtattCGCGCATGGTTGGATAAAATGAATAAATGGAGCAGGTTAAACCATTGTTCTCTTTTGGTAAGTAATCCCGAAAATAATAACGTTAAACAATTTTCATTGTTGCTTGAGGAATACCGTTCTCTTTTTGGTCTTGCCAGTTTGCGTTTTCAGGGTGACCAACATTTGCTGGATAGTGCCTTCTGGTGCAACGAAAAAGGGGTCAGCGCCCGTCAGCAGCTTAGCGTTCAGCAACAAAATGGTagc                                                                    <  2:337762/245‑3 (MQ=255)
                                                                                                                                                                          cGCTCAAGGTTGGATAAAATGAATAAATGGGGCAGGTTAAACCATTGTTCGCTTTTGGTAATTAATCCCAAAAATAATAACGATAAACAATGTTCATTGTTGCTTTAGGAATACCGTTCTCTTTTTGGTCTGGCCAGTTTGCGTTTTCAGGGTGACCAACATTTGCTGGATATTGCCTGCTGGTGCAACGAAAAAGGGGTCAGCGCCCGTCAGCAGCTTAGCGTTCAGTAACAAAATGGTATCTGGACAt                                                             <  2:292969/250‑1 (MQ=255)
                                                                                                                                                                             tCATGGTTGGATAAAATGAATAAATGGAGCAGGTTAAACCATTGTTCGCTTTTGGTAATTAATCCCGAAAATAATAACGATAAACAATTTTCATTGTTGCTTGAGGAATACCGTTCACTTTTTGGTCTTGCCAGTTTGCGTTTTCAGGGAGACCAACATTTGCTGGATATTGCCTTCTGGTGCAACGAAAAAGGGGTCAGCGCCCGTCAGCAGCTTAGCGTTCAGCAACAAAATGGTATCTGGACATTAGt                                                         >  1:809574/1‑251 (MQ=255)
                                                                                                                                                                                 ggTTGGAGAAAATGAATAAATGGAGCAGGTTAAACCATTGTTCGCTTTTGGTAATTAATCCCGAAAATAATAACGATAAACAATTTTCATTGTTGCTTGAGGAATACCGTTCACTTTTTGGTCTTGCCAGTTTGCGTTTTCAGGGTGACCAACATTTGCTGGATATTGCCTTCTGGTGCAACGAAAAAGGGGTCAGCGCCCGTCAGCAGCTTAGCGTTCAGCAACAAAATGGTATCTGGACATTAGTTCa                                                      <  2:222314/250‑1 (MQ=255)
                                                                                                                                                                                          aaaTGAATAAATGGAGCAGGTTAAACCATTGTTCGCTTTTGGTAATTAATCCCGAAAATAATAACGATAAACAATTTTCATTGTTGCTTGAGGAATACCGTTCACTTTTTGGTCTTGCCAGTTTGCGTTTTCAGGGTGACCAACATTTGCTGGATATTGCCTTCTGGTGCAACGAAAAAGGGGTCAGCGCCCGTCAGCAGCTTAGCGTTCAGCCACAAAATGGTATCTGGACATTAGTTCAAAGCGAag                                              >  1:626673/1‑249 (MQ=255)
                                                                                                                                                                                             tGAATAAATGGAGCAGGTTAAACCATTGTTCGCTTTTGGTAATTAATCCCGAAAATAATAACGATAAACAATTTTCATTGTTGCTTGAGGAATACCGTTCACTTTTTGGTCTTGCCAGTTTGCGTTTTCAGGGTGACCAACATTTGCTGGATATTGCCTTCTGGTGCAACGAAAAAGGGGTCAGCGCCCGTCAGCAGCTTAGCGTTCAGCAACAAAATGGTATCTGGACATTAGTTCAAAGCGAAGAGGCg                                         >  1:748220/1‑251 (MQ=255)
                                                                                                                                                                                              gAATAAATGGAGCAGGTTAAACCATTGTTCGCTTTTGGTAATTAATCCCGATAATAATTACTATAAACAATTTTCATTGTTGCTTTAGGAATACCGTTCACTTTTTGTTCTTGCCAGTTTGCGTTTTCAGGGTGACCAACATTTGCTGGATATTGCCTTCTGGTGCACCGAAAAAGGGGTCAGCGCCCGTCAGCAGCTTAGCGTTCAGCAACAAAATGGTATCTGGACATTAGTTCAAAGCGAAGAggagg                                        <  2:688214/251‑4 (MQ=255)
                                                                                                                                                                                              gAATAAATGGAGCAGGTTAAACCATTGTTCGCTTTTGGTAATTAATCCCGAAAATAATAACGATAAACACTTTTCATTGTTGCTTGAGGAATACCGTTCACTTTTTGGACTTGCCAGTTTGCGTTTTCAGGGTGACCAACATTTGCTGGATATTGCATTCTGGTGCAACGAAAAAGGGGTCAGAGCCCGTCAGCAGCTTAGCGTTCAGAAACAAAATGGTATCTGGACATTAGTTCAAAGCGAAGAGGCgg                                        >  2:42177/1‑251 (MQ=255)
                                                                                                                                                                                                aTAAATGGAGCAGGTTAAACCATTGTTCGCTTTTGGTAATTAATCCCGAAAATAATAACGATAAACAATTTTCATTGTTGCTTGAGGAATACCGTTCACTTTTTGGTCTTGCCAGTTTGCGTTTTCAGGGTGACCAACATTTGCTGGATATTGCCTTCTGGTGCAACGAAAAAGGGGTCAGCGCCCGTCAGCAGCTTAGCGTTCAGCAACAAAATGGTATCTGGACATTAGTTCAAAGCGAAGAGGCgg                                        >  1:696102/1‑249 (MQ=255)
                                                                                                                                                                                                  aaaTGGAGCAGGTTAAACCATTGTTCGCTTTTGGTAATTAATCCCGAAAATAATAACGATAAACAATTTTCATTGTTGCTTGAGGAATACCGTTCACTTTTTGGTCTTGCCAGTTTGCGTTTTCAGGGTGACCAACATTTGCTGGATATTGCCTTCTGGTGCAACGAAAAAGGGGTCAGCGCCCGTCAGCAGCTTAGCGTTCAGCAACAAAATGGTATCTGGACATTAGTTCAAAGCGAAGAGGCGGAGAt                                    <  1:89621/251‑1 (MQ=255)
                                                                                                                                                                                                       gAGCAGGTTAAACCATTGTTCGCTTTTGGTAATTAATCCCGAAAATAATAACGATAAACAATTTTCATTGTTGCTTGCGGAATACCGTTCACTTTTTGGTCTTGCCAGTTTGCGTTTTCAGGGTGACCAACATTTGCTGGATATTGCCTTCTGGTGCAACGAAAAAGGGGTCAGCGCCCGTCAGCAGCTTAGCGTTCAGCAACAAAATGGGATCTGGCCATTAGTTCAAAGCGAAGAGGCGGAGATCAAAc                               >  2:598621/1‑251 (MQ=255)
                                                                                                                                                                                                            tgTTAAACCATTGTTCGCTTTTGGTAATTAATCCCGAAAATAATAACGATAAACAATTTTCATTGTTGCTTGAGGAATACCGTTCACTTTTTGGTCTTGCCAGTTTGCGTTTTCAGGGTGACCACCATTTGCTGGATATTGCCTTCTGGTGCAACGAAAAAGGGTTCAGCGCCCGTCAGCAGCTTAGCGTTCAGCAACAAAATGGTATCTGGCCTTTAGTTCAAAGCGAAGAGGCGGAGATCCAACAACGc                          >  2:656900/2‑251 (MQ=255)
                                                                                                                                                                                                            ggTTAAACCATTGTTCGCTTTTGGTAATTAATCCCGAAAATAATAACGATAAACAATTTTCATTGTTGCTTGAGGAATACCGTTCACTTTTTGGACTTGCCAGTTTGCTTTTTCGGGGGGACCACAATTTCCTGGATATTGCCTTCTGGTCCAACGAAAAAGGGGTCAGCGCCGGTCAGCAGCTTACCGTTCAGCAACAAAATGGTATCTGGCCTTTAGTTCAAAGCGAAGAGGCGGAAATCCAACCACGc                          >  2:429418/1‑251 (MQ=255)
                                                                                                                                                                                                                aaaCCATTGTTCGCTTTTGGTAATTAATCCCGAAAATAATAACGATAAACAATTTTCATTGTTGCTTGAGGAATACCGTTCACTTTTTGGTCTTGCCAGTTTGCGTTTTCAGGGTGACCAACATTTGCTGGATATTGCCTTCTGGTGCAACGAAAAAGGGGTCAGCGCCCGTCAGCAGCTTAGCGTTCAGCAACAAAATGGTATCTGGACATTAGTTCAAAGCGAAGAGGCGGAGATCCAACCACGCAGcg                      >  2:77351/1‑251 (MQ=255)
                                                                                                                                                                                                                    cATTGTTCGCTTTTGGTAATTAATCCCGAAAATAATAACGATAAACAATTTTCATTGTTGCTTGAGGAATACCGTTCACTTTTTGGTCTTGCCAGTTTGCGTTTTCAGGGTGACCAACATTTGCTGGATATTGCCTTCTGGTGCAACGAAAAAGGGGTCAGCGCCCGTCAGCAGCTTAGCGTTCAGCAACAAAATGGTATCTGGACATTAGTTCAAAGCGAAGAGGCGGAGATCCAACCACGCAGcgacg                   >  1:609517/1‑250 (MQ=255)
                                                                                                                                                                                                                           cGCTTTTGGTAATTAATCCCGAAAATAATAACGATAAACAATTTTCATTGTTGCTTGAGGAATACCGTTCACTTTTTGGTCTTGCCAGTTTGCGTTTTCAGGGTGACCAACATTTGCTGGATATTGCCTTCTGGTGCAACGAAAAAGGGGTCAGCGCCCGTCAGCAGCTTAGCGTTCAGCAACAAAGTGGTATCTGGCCATTAGTTCAAAGCGAAGAGGCGGAGATACAACCACGACGCGACGAAAAACGc           >  2:201984/1‑251 (MQ=255)
                                                                                                                                                                                                                               tttGGTAATTAATCCCGAAAATAATAACGATAAATAATTTTCATTGTTGCTTGAGGAATACCGTTCACTTTTTGGTCTTGCCAGTTTGCGTTTTCAGGGTGACCAACATTTGCTGGATATTGCCTTCTGGTGCAACGAAAAAGGGGTCAGCGCCCGTCAGCAGCTTAGCGTTCAGCAACAAAATGGTATCTGGACATTAGTTCAAAGCGAAGAGGCGGAGATCCAACCACGCAGCGACGAAAAACGCAt         <  1:238698/249‑1 (MQ=255)
                                                                                                                                                                                                                                   gTAATTAATCCCGAAAATAATAACGATAAACAATTTTCATTGTTGCTTTAGGAATACCGTTCACTTTTTGGTCTTGCCAGTTTGCGTTTTCAGGGTGACCAACATTTGCTGGATATTGCCTTCTGGTGCAACGAAAAAGGGGTCAGCGCCCGTCAGCAGCTTAGCGTTCAGCAACAAAATGGTATCTGGACATTAGTTCAAAGCGAAGAGGCGGAGATCCAACCACGCAGCGACGAAAAACGCATTCTGAg   <  1:589182/251‑1 (MQ=255)
                                                                                                                                                                                                                                     aattaatCCCGAAAATAATAACGATAAACAATTTTCATTGTTGCTTGAGGAATACCGTTCACTTTTTGGTCTTGCCAGTTTGCGTTTTCAGGGTGACCAACATTTGCTGGATATTGCCTTCTGGTGCAACGAAAAAGGGGTCAGCGCCCGTCAGCATCTTAGCGTTCAGCAACAAAATGGTATCTGGACATTAGTTCAAAGCGAAGAGGCGGAGATCCAACCACGCAGCGACGAAAAACGCATTCTGAga  >  1:681841/1‑249 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                              
TAGATGATTCTCAAGGGCCGGAAAAAATAAAATTATTTTCAATGCTAAATCATGAAAAAGGTCTATACTATTTGACCCGTGATTTGCAGTGTTCTATTGATCCCCATAATTACCTTTTTATTCTTGTTTGCGCAAATAACGCATGGCAAAACATTCCTGCCGAGCGGCTTCGCTCATGGTTGGATAAAATGAATAAATGGAGCAGGTTAAACCATTGTTCGCTTTTGGTAATTAATCCCGGAAATAATAACGATAAACAATTTTCATTGTTGCTTGAGGAATACCGTTCACTTTTTGGTCTTGCCAGTTTGCGTTTTCAGGGTGACCAACATTTGCTGGATATTGCCTTCTGGTGCAACGAAAAAGGGGTCAGCGCCCGTCAGCAGCTTAGCGTTCAGCAACAAAATGGTATCTGGACATTAGTTCAAAGCGAAGAGGCGGAGATCCAACCACGCAGCGACGAAAAACGCATTCTGAGT  >  NC_000913/3696669‑3697147

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