Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,561,368 C→T P147S (CCG→TCG)  kptA → RNA 2'‑phosphotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,561,3680CT96.5% 179.1 / NA 58P147S (CCG→TCG) kptARNA 2'‑phosphotransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (34/22);  minor base G (0/2);  total (34/24)
Fisher's exact test for biased strand distribution p-value = 1.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAAAAGCGTTTTAGTTATTCCAGTGATGGTAGGTGCATCCGCGCGGTTCAGGGGCATTCGACTTCGCAGGTCGCTATTTCCTTTGCCGAAAAAA‑CACCGCCACAGTTTCTCTACCACGGTACGGCAAGCCGTTTCCTTGACGAGATAAAAAAACAGGGGCTGATTGCAGGTGAGCGTCACTACGTACATCTCTCTGCTGATGAAGCAACGGCACGCAAAGTGGGAGCAAGACACGGATCGCCGGTTATTTTAACCGTCAA‑AGCACAGGAAATGGCGAAACGAGGCCTCCCTTTCTGGCAAGCGGAAAACGGGGTCTGGCTGACATCAACGGTTGCGGTTGAGTTTCTTGAGTGGTAAGGGAAACTTCCCTCATGAAGTGACCACTTGTAGCTTAAATAATTGCACTGTGCTGATAAGCGCAGCAATAAACAACGCCGCCGCGCCGCAGAGAGTGGCGCTGGCAAGCGCCGATGTTCCATT  >  NC_000913/4561128‑4561607
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                
aaaaaGCGTTTTGGTTATTCCAGTGATGGTCGGTGCATCCGCGCGGTTCAGGGGCATTCGACTTCGAAGGTCGCTATTACCTTTGCCGAAAACA‑CACCGCCACAGTTTCTCTACCACGGTACGGCAAGCCGTTTCATTGACGAGATAAAAAAACAGGGGCTGATTGCAGGTGAGCGTCACTACGTACATCTCTCTGCTGATGAAGCAACGGCACGAACAGAGGGAGCAAGACCCGGCTCGTCGGTTAtttt                                                                                                                                                                                                                                        >  2:565341/1‑251 (MQ=255)
      cGTTTTAGTTATTCCAGTGATGGTAGGTGCATCCGCGCGGTTCAGGGGCATTCGACTTCGCAGGTCGCTATTTCCTTTGCCGAAAAAA‑CACCGCCACAGTTTCTCTACCACGGTACGGCAAGCCGTTTCCTTGACGAGATAAAAAAACAGGGGCTAATTGCAGGTGAGCGTCACTACGTACATCTCTCTGCTGATGAAGCAACGGCACGCAAAGTGGGAGCAAGACACGGATCGTCGGTTATTTTAACCGt                                                                                                                                                                                                                                  >  1:70373/1‑251 (MQ=255)
       gTTTTAGTTATTCCAGAGATGGTAGGTGCATCCGCGCGGTTCAGGGGCATTCGACTTCGCAGGTCGCTATTTCCTTTGCCGAAAAAA‑CACCGCCACAGTTTCTCTACCACGGTACGGCAAGCCGTTTCCTTGACGAGATAAAAAAACAGGGGCTGATTGCAGGTGAGCGTCACTACGTACATCTCTCTGCTGATGAAGCAACGGCACGCAAAGTGGGAGCAAGACACGGATCGTCGGTTATTTTAACCGTc                                                                                                                                                                                                                                 >  1:416316/1‑251 (MQ=255)
                aTTCCAGTGATGGTAGGTGCATCCGCGCGGTTCAGGGGCCTTCGACTTCGCAGGACGCTATTTCCTTTGCCGAAAAAA‑CACCGCCACAGTTTCTCTACCACGGTACGGCAAGCCGTTTCCTTGACGAGATAAAAAAACAGGGGCTGCTTGCAGGTGAGCGTCCCTACGTACATCTCTCTGATGAGGACGCAACGTCACGCAAAGTGGGGGCAAGACACGGATCGTCGGGTATTTTAAACGTCaa‑atccacg                                                                                                                                                                                                                       >  2:55794/1‑245 (MQ=255)
                 ttCCAGTGATGGGAGGTGCATCCGCGCGGTTCAGGGGCATTCGACTTCGCAGGTCGCTATTTCCTTTGCCGAAAAAA‑CACCGCCACAGTTTCTCTACCACGGTACGGCAAGCCGTTTCCTTGACGAGATAAAAAAACAGGGGCTGATTGCAGGCGAGCGTCACTACGTACATCTCTCTGCGGATGAAGCAACGGCACGCAAAGTGGGAGCAAGCCACGGATCGTCGGTTATTTTACCCGTACA‑AGCACAgg                                                                                                                                                                                                                      >  2:621814/1‑251 (MQ=255)
                  cccggtgttgtgagcttgaTCCGCGCGGTCCAGGGCAATTCGCTTTCGCTGGTCGTTATCCCCTTTGC‑GAAAAAATCCCCGCCACAGTTTCTCTACCACGGTACGGCAAGCCGTTTCCGTGACGAGATAAAAAAACAGGGGCTGATTGCAGGTGAGCGTCACTACGTACATCTCTCTGCTGATGAAGCAACGGCACGCAAAGTGGGAGCAAGACACGTATCGTCGGTTATTTTAACCGTCAA‑AGcaca                                                                                                                                                                                                                        <  1:359748/230‑1 (MQ=255)
                   ccAGTGATGGTAGGTGCATCCGCGCGGTTCCGGGGCATTCGACTTCGCAGGTCGCTATTTCCTTTGCCGAAAAAA‑CACCGCCACAGTTTCTCTACCACGGTACGGCAAGCCGTTTCCTTGACGAGATAAAAAAACAGGGGCTGATTGCAGGTGAGCGTCACTACGTACATCTCTCTGCTGATGAAGCAACGGCACGCAAAGTGGGAGCAAGACACGGATCGTCGGTTATTTTAACCGTCAA‑AGCACAGGa                                                                                                                                                                                                                     >  2:119948/1‑250 (MQ=255)
                        gATGGTAGGGGCATCCGCGCGGTTCAGGGGCATTCGACTTCGCAGGTCGCTATTTCATTTGCCGAAAAAA‑CACCGCCACAGTTTCTCTACCACGGTACGGCAAGCCGTTTCCTTGACGAGATAAAAAAACAGGGGCTGATTGCAGGTGAGCGTCAATACGTACATCTCTCTGCTGTTGAAGCAACGGCACGCAAAGAGGGAGCAAGAAACGGATCGTCGGGTATTTTAACCGTCAA‑AGCACAGGCAATggg                                                                                                                                                                                                               >  2:625116/1‑250 (MQ=255)
                               ggTGCATCCGCGCGGTTCAGGGGCATTAGACTTCGCAGGTCGCTATTTCCTTTGCCGAAAAAA‑CACCGCCACAGTTTCTCTACCACGGTACGGCAAGCCGTTTCCTTGACGAGATAAAAAAACAGGGGCTGATTGCAGGTGAGAGTCACTACGTACATCTCTCTGCTGATGAAGCAACGGCACGCAAAGTGGGAGCAAGACACGGATCGTCGGTTATTTTAACCGTCAA‑AGCACAGGAAATGGAGAAACg                                                                                                                                                                                                         >  2:645096/1‑250 (MQ=255)
                                gTGCATCCGCGCGGTTCAGGGGCATTCGACTTCGCAGGTCGCTATTTCCTTTGCCGAAAAAA‑CACCGCCACAGTTTCTCTACCACGGTACGGCAAGCCGTTTCCTTGACGAGATAAAAAAACAGGGGCTGATTGCAGGTGAGCGTCACTACGTACATCTCTCTGCTGATGAAGCAACGGCACGCAAAGTGGGAGCAAGACACGGATCGTCGGTTATTTTAACCGTCAA‑AGCACAGGAAATGGCGAAACg                                                                                                                                                                                                         >  1:379237/1‑249 (MQ=255)
                                   cATCCGCGCGGTTCAGGGGCATTCGACTTCGCAGGTCGCTATTTCCTTTGCCGAAAAAA‑CACCGCCACAGTTTCTCTACCACGGTACGGCAAGCCGTTTCCTTGACGAGATAAAAAAACAGGGGCTGATTGCAGGTGAGCGTCACTACGTACATCTCTCTGCTGATGAAGCAACGGCACGCAAAGTGGGAGCAAGACACGGATCGTCGGTTATTTTAACCGTCAA‑AGCACAGGAAATGGCGAAACGAGGcc                                                                                                                                                                                                    >  1:408191/1‑251 (MQ=255)
                                    aTCCGCGCGGTTCAGGGGCATTCGACTTCGCAGGTCGCTATTTCCTTTGCCGAAAAAA‑CACCGCCACAGTTTCTCTACCACGGTACGGCAAGCCGTTTCCTTGACGAGATAAAAAAACAGGGGCTGATTGCAGGTGAGCGTCACTACGTACATCTCTCTGCTGATGAAGCAACGGCACGCAAAGTGGGAGCAAGACACGGATCGTCGGTTATTTTAACCGTCAA‑AGCACAGGCAATGGCGAAACGAGGCCt                                                                                                                                                                                                   >  1:22038/1‑251 (MQ=255)
                                          gcgGTTCAGGGGCATTCGACTTCGCAGGTCGCTATTTCCTTTGCCGAAAAAA‑CACCGCCACAGTTTCTCTACCACGGTACGGCAAGCCGTTTCCTTGACGAGATAAAAAAACAGGGGCTGATTGCAGGTGAGCGTCACTACGTACATCTCTCTGCTGATGAAGCAACGGCACGCAAAGTGGGAGCAAGACACGGATCGTCGGTTATTTTAACCGTCAA‑AGCACAGGAAATGGCGAAACGAGGCCTCCCttt                                                                                                                                                                                             >  1:298627/1‑251 (MQ=255)
                                                           gACTTCGCAGGTCGCTATTTCCTTTGCCGAAAAAA‑CACCGCCACAGTTTCTCTACCACGGTACGGCAAGCCGTTTCCTTGACGAGATAACAAAACAGGGGCTGATTGCAGGTGAGCGTCACTACGTACATCTCTCTGCTGATGAAGCAACGGCACGCAAAGTGGGAGCAAGACACGGATCGTCGGTTATTTTAACCGTCAG‑AGCACAGGAAATAGCGCAACGAGGCCTCCCTTTCTGGAAAGAGGAAAACg                                                                                                                                                                            >  2:8248/1‑251 (MQ=255)
                                                            aCTTCGCAGGTCGCTATTTCCTTTGCCGAAAAAA‑CACCGCAACAGTTTCTCTACCACGGTACGGCAAGCCGTTTCCTTGACGAGATAAAAAAACAGGGGCTGATTGCAGGTGAGCGTCACTACGTACATCTCTCTGCTGATGAAGCAACGGCACGCAAAGTGGGAGCAAGACACGGATCGTCGGTTATTTTAACCGTCAA‑AGCACAGGAAATGGCGAAACGAGGCCTCCCTTTCTGGCAAGCGGAAACCgg                                                                                                                                                                           >  2:351560/1‑251 (MQ=255)
                                                                     gTCGCTATTTCCTTTGCCGAAAAAA‑CACCGCCACAGTTTCTCTACCACGGTACGGCAAGCCGTTTCCTTGACGAGATAAAAAAACAGGGGCTGATTGCAGGTGAGCGTCACTACGTACATCTCTCTGCTGATGAAGCAACGGCACGCAAAGTGGGAGCAAGACACGGATCGTCGGtt                                                                                                                                                                                                                                             >  1:375628/1‑177 (MQ=255)
                                                                     gTCGCTATTTCCTTTGCCGAAAAAA‑CACCGCCACAGTTTCTCTACCACGGTACGGCAAGCCGTTTCCTTGACGAGATAAAAAAACAGGGGCTGATTGCAGGTGAGCGTCACTACGTACATCTCTCTGCTGATGAAGCAACGGCACGCAAAGTGGGAGCAAGACACGGATCGTCGGTTATTTTAACCGTCAA‑AGCACAGGAAATGGCGAAACGAGGCCTCCCTTTCTGGCACGCGGAAAACGGGGTCTGGCt                                                                                                                                                                  >  2:199129/1‑251 (MQ=255)
                                                                          tATTTCCTTTGCCGAAAAAA‑CACCGCCACAGTTTCTCTACCACGGTACGGCAAGCCGTTTCCTTGACGAGATAAAAAAACAGGGGCTGATTGCAGGTGAGCGTCACTACGTACATCTCTCTGCTGATGAAGCAACGGCACGCAAAGTGGGAGCAAGACACGGATCGTCGGTTATTTTAACCGTCAA‑AGCACAGGAAATGGCGAAACGAGGCCTCCCTTTCTGGCAAGCGGAAAACGGGGTCTGGCTGACAt                                                                                                                                                             >  1:359866/1‑251 (MQ=255)
                                                                           tttcctgttgacaaaaaaa‑CCCCGCCACAGTTTCTCTACCACGGTACGGCAAGCCGTTTCCTTGACGAGATAAAAAAACAGGGGCTGATTGCAGGTGAGCGTCACTACGTACATCTCTCTGCTGATGAAGCAACGGCACGCAAAGTGGGAGCAAGACACGGATCGTCGGTTATTTTAACCGTCAA‑AGCACAGGAAATGGCGAAACGAGGCCTCACTTTCTGGCAAGCGGAAAACGGGGTCTGGCTg                                                                                                                                                                 <  2:30954/233‑1 (MQ=255)
                                                                                tgtggcaaaaaaaa‑TACCGCCACAGTTTCACTACAACGGTACGGCAAGCCGTTTCCTTGACGAGATAAAAAAACTGGGGCTGATTGCCGGTGAGCGTCACTACGTACATCTCTCTGCTGATGAAGCTACGGCACGCAAAGTGGGAGCAAGACACGGATCGTCGGTTATTTTAACCGTCAA‑AGCACAGGAAATGGCGAAACGAGGCCGCCCTTTCTGGCAAGCGGAAAACGGGGTCTGGCTGACATCaa                                                                                                                                                          <  2:499975/240‑1 (MQ=255)
                                                                                    cccaaaaaac‑CACCGCCACAGTTTCTCTACCACGGTACGGCAAGCCGTTTCCTTGACGAGATAAAAAAACAGGGGCTGATTGCAGGTGAGCGTCACTACGTACATCTCTCTGCTGATGAAGCAACGGCACGCAAAGTGGGAGCAAGACACGGATCGTCGGTTATTTTAACCGTCAA‑AGCACAGGAAATGGCGAAACGAGGCCTCCCTTTCTGGCAAGCGGAAAACGGGGTCTGGCTGACATCAAcggttgc                                                                                                                                                   <  1:316794/247‑1 (MQ=255)
                                                                                     caaaaaaaa‑CACCGCCACAGTTTCTCTACCACGGTACGGCAAGCCGTTTCCTTGACGAGATAAAAAAACAGGGGCTGATTGCAGGTGAGCGTCACTACGTACATCTCTCTGCTGATGAAGCAACGGCACGCAAAGTGGGAGCAAGACACGGATCGGCGGTTATTTTAACCGTCAA‑AGCACAGGAAATGGCGAAACGAGGCCTCCCTTTCTGGCAAGCGGAAAACGGGGTCTGGCTGACATCAAcggttgcg                                                                                                                                                  <  2:70373/248‑1 (MQ=255)
                                                                                             a‑cacCGCCACAGTTTCTCTACCACGGTACGGCAAGCCGTTTCCTTGACGAGATAAAAAAACAGGGGCTGATTGCAGGTGAGCGTCACTACGTACATCTCTCTGCTGATGAAGCAACGGCACGCAAAGTGGGAGCAAGACACGGATCGTCGGTTATTTTAACCGTCAA‑AGCAAAGGAAATGGCGAAACAAGGCCTCCCTTTCTGGCAAGCGGAAAACGGGGGCTGGCTGCCATCAACGGATGCGGTTGAGtt                                                                                                                                          >  2:655916/1‑251 (MQ=255)
                                                                                               ccccGCCACAGTTTCTCTACCACGGTACGGCAAGACGTTTCCTTGACGAGATTAAAAAACAGGGGCTGATGGCAGGTGAGCGTCACTACGTCCATCTCTCTGCTGATGAAGCAACGGCACGCAAAGTGGGAGCAAGACACGGATCGTCGGTTATTTTAACCGTCAA‑AGCACAGGAAATGGCGAAACGAGGCCTCCCTTTCTGGCAAGCTGAAAACGGGGTCTGGCTGACATCAACGGTTGCGGTTGAGttt                                                                                                                                         <  2:597725/249‑1 (MQ=255)
                                                                                                   gCCACAGTTTCTCTACCACGGTACGGCAAGCCGTTTCCTTGACGAGATAAAAAAACAGGGGCTGATTGCAGGTGAGCGTCACTACGTACATCTCTCTGCTGATGAAGCAACGGCACGCAAAGTGGGAGCAAGACACGGATCGTCGGTTATTTTAACCGTCAA‑AGCACAGGAAATGGCGAAACGAGGCCTCCCTTTCTGGCAAGCGGAAAACGGGGTCTGGCTGACATCAACGGTTGCGGTTGGGTTTCTTg                                                                                                                                     >  1:22122/1‑251 (MQ=255)
                                                                                                        aGTTTCTCTACCACGGTACGGCAAGCCGTTTCCTTGACGAGATAAAAAAACAGGGGCTGATTGCAGGTGAGCGTCACTACGTACATCTCTCTGCTGATGAAGCAACGGCACGCAAAGTGGGAGCAAGACACGGATCGTCGGTTATTTTAACCGTCAA‑AGCACAGGAAATGGCGAAACGAGGCCTCCCTTTCTGGCAAGCGGAAAACGGGGTCTGGCTGACATCAACGGTTGCGGTTGAGTTTCTTGAGTgg                                                                                                                                <  2:582451/251‑1 (MQ=255)
                                                                                                              tctACCACGGTACGGCAAGCCGTTTCCTTGACGAGATAAAAAAACAGGGGCTGATTGCAGGTGAGCGTCACTACGTACATCTCTCTGCTGATGAAGCAACGGCACGCAAAGTGGGAGCAAGACACGGATCGTCGGTTATTTTAACCGTCAA‑AGCACAGGAAATGGCGAAACGAGGCCTCCCTTTCTGGCAAGCGGAAAACGGGGTCTGGCTGACATCAACGGTTGCGGTTGGGTTTCTTGAGTGGTAAggg                                                                                                                          >  2:294430/1‑251 (MQ=255)
                                                                                                                       gttcGGCAAGCCGGTTCCTTGACGAGATAAAAAAACAGGGGCTGATTGCAGGTGAGCGTCACTACGTACATCTCTCTGCTGATGAAGCAACGGCACGCAAAGTGGGAGCAAGACACGGATCGTCGGTTATTTTAACCGTCAA‑AGCACAGGAAATGGTGAAACGAGGCCTCCCTTTCTGGCAAGCGGAAAACGGTGTCTGGCTTACATCAACTGTTGCGGTTGAGTTTCTTGAGTGGTAAGGGAAACTTccc                                                                                                                 <  1:310976/248‑1 (MQ=255)
                                                                                                                        tACGGCAAGCCGTTTCCTTGACGAGATAAAAAAACAGGGGCTGATTGCAGGTGAGCGTCACTACGTACATCTCTCTGCTGATGAAGCAACGGCACGCAAAGTGGGAGCAAGACACGGATCGTCGGTTATTTTAACCGTCAA‑AGCACAGGAAATGGCGAAACGAGGCCTCCCTTTCTGGCAAGCGGAAAACGGGGTCTGGCTGACATCAACGGTTGCGGTTGAGTTTCTTGAGGGGTAAGGGAAACTTCCCt                                                                                                                >  2:620936/1‑251 (MQ=255)
                                                                                                                                   gTTTCCTTGACGAGATAAAAAAACAGGGGCTGATTGCAGGTGAGCGTCACTACGTACATCTCTCTGCTGATGAAGCAACGGCACGCAAAGTGGGAGCAAGACACGGATCGTCGGTTATTTTAACCGTCAA‑AGCACAGGAAATGGCGAAACGAGGCCTCCCTTTCTGGCAAGCGGAAAACGGGGt                                                                                                                                                                        <  2:328066/184‑1 (MQ=255)
                                                                                                                                      tccgtGACGAGATAAAAAAACAGGGGCTGATTGCATGTGAGCGTCACTACGTACATCTCTCTGCTGATGAAGCAACGGCACGCAAAGTGGGAGCAAGACACGGATCGTCGGTTATTTTAACCGTCAA‑AGCACAGGAAATGGCGAAACGAGGCCTCCCTTGCTCGCAAGCGGAAAAAGGGGTCTGGCt                                                                                                                                                                  <  2:341364/183‑1 (MQ=255)
                                                                                                                                      tCCTTGACGAGATAAAAAAACAGGGGCTGATTGCAGGTGAGCGTCACTACGTACATCTCTCTGCTGATGAAGCAACGGCACGCACAGTGGGAGCAAGACACGGATCGTCGGTTATTTTAACCGTCAA‑AGCACAGGAAATGGCGAAACGAGGCCTCCCTTTCTGGCAAGCGGAAAACGGGGTCTGGCt                                                                                                                                                                  >  1:341364/1‑187 (MQ=255)
                                                                                                                                       ccTTGACGAGATAAAAAAACAGGGGCTGATTGCAGGTGAGCGTCACTACGTACATCTCTCTGCTGATGAAGCAACGGCACGCAAAGTGGGAGCAAGACACGGATCGTCGGTTATTTTAACCGTCAA‑AGCACAGGAAATGGCGAAACGAGGCCTCCCTTTCTGGCAAGCGGAAAACGGGGTCTGGCTGACATCAACGGTTGCGGTTGAGTTTCTTGAGTGGTAAGGGAAACTTCCCTCATGAAGTGACCACt                                                                                                 <  1:310442/251‑1 (MQ=255)
                                                                                                                                          tGACGAGATAAAAAAACAGGGGCTGATTGCAGGTGAGCGTCACTACGTACATCTCTCTGCTGATGAAGCAACGGCACGCAAAGTGGGAGCAAGACACGGATCGTCGGTTATTTTAACCGTCAA‑AGCACAGGAAATGGCGAAACGAGGCCTCCCTTTCTGGCAAGCGGAAAACGGGGTCTGGCTGACATCAACGGTTGCGGTTGAGTTTCTTGAGTGGTAAGGGAAACTTCCCTCATGAAGTGACCACtt                                                                                                <  1:554777/249‑1 (MQ=255)
                                                                                                                                           gACGAGATAAAAAAACAGGGGCTGATTGCAGGTGAGCGTCACTACGTACATCTCTCTGCTGATGAAGCAACGGCACGCAAAGTGGGAGCAAGACACGGATCGTCGGTTATTTTAACCGTCAA‑AGCACAGGAAATGGCGAAACGAGGCCTCCCTTTCTGGCAAGCGGAAAACGGGGTCTGGCTGACATCAACGGTTGCGGTTGAGTTTCTTGAGTGGTAAGGGAAACTTCCCTCATGAAGTGACCACTTGt                                                                                              <  1:294430/250‑1 (MQ=255)
                                                                                                                                            aCGAGATAAAAAAACAGGGGCTGATTGCAGGTGAGCGTCACTACGTACATCTCTCTGCTGATGAAGCAACGGCACGCAAAGTGGGAGCAAGACACGGATCGTCGGTTATTTTAACCGTCAA‑AGCACAGGAAATGGCGAAACGAGGCCTCCCTTTCTGGCAAGCGGAAAACGGGGTCTGGCTGACATCAACGGTTGCGGTTGAGTTTCTTGAGTGGTAAGGGAAACTTCCCTCATGAAGTGACCACTTGTAg                                                                                            <  1:407510/251‑1 (MQ=255)
                                                                                                                                                    aaaaaaCAGGGGCTGATTGCAGGTGAGCGTCACTACGTACATCTCTCTGCTGATGAAGCAACGGCACGCAAAGTGGGAGCAAGACACGGATCGTCGGTTATTTTAACCGTCAA‑AGCACAGGAAATGGCGAAACGAGGCCTCCCTTTCTGGCAAGCGGAAAACGGGGTCTGGCTGACATCAACGGTTGCGGTTGAGTTTCTTGAGTGGTAAGGGAAACTTCCCTCATGAAGTGACCACTTGTAGCTTAaata                                                                                    <  1:199129/251‑1 (MQ=255)
                                                                                                                                                     caaaaCATGGGCTAATTGCAGGTGAGCGTCACTACGTACATCTCTCTGCTGATGAAGCAACGGCACGCAAAGTGGGAGCAAGACACGGATCGTCGGGTATTTTAACCGTCAA‑AGCACAGGCAATGGCGCAACGAGGCCTCCCTTTCTGGCAAGCGGAAAACGGGGTCTGGCTGAAATCAACGGTTGAGGTTGGGTTTCTTGGGTGGGAAAGGAAGCTTCCCTACTGAAGAGACAACTTgaggcttacata                                                                                    >  2:380198/2‑239 (MQ=255)
                                                                                                                                                        aaCAGGGGCTGATTGCAGGTGAGCGTCACTACGTACATCTCTCTGCTGATGAAGCAACGGCACGCAAAGTGGGAGCAAGACACGGATCGTCGGTTATTTTAACCGTCAA‑AGCACAGGAAATGGCGAAACGAGGCCTCCCTTTCTGGCAAGCGGAAAACGGGGTCTGGCTGACATCAACGGTTGCGGTTGAGTTTCTTGAGTGGTAAGGGAAACTTCCCTCATGAAGTGACCACTTGTAGCTTAAATAATTg                                                                                <  1:645096/251‑1 (MQ=255)
                                                                                                                                                             ggtcTGAGTGCAGGTTAGCGTCACTACCTTCTACTCTCTGCTGATGAAGCACCGGCACGCAAAGTGGGAGCAAGACACGGATCTTCTGTTATTTTAACCGTCAA‑AGCACAGGAAATGTCGAATCGAGGCTTCCCTTTCTTGCAAGCTGAAAACGGGGTCTGGCTGACATCAACTGTTGCGGTTGAGTTTCTTGAGTGGTAAAGGAAACTTCCCTCATGAAGTGACCACTTTTAGCTGAAAGAATTGCACtg                                                                           <  2:593332/248‑1 (MQ=255)
                                                                                                                                                                cTGATTGCAGGTGAGCGTCACTACGTACATCTCTCTGCTGATGAAGCAACGGCACGCAAAGTGGGAGCAAGACACGGATCGTCGGATATTTTAACCGTCAA‑AGCACAGGAAATGGCGAAACGAGGCCTCCCTTTCTGGCAAGCGGAAAACGGGGTCTGGCTGACATCAACGGTTGCGGTTGAGTTTCTTGAGTGGTAAGGGAAACTTCCCTCATGAAGTGACCACTTGTAGCTTAAATAATTGCACTGTGc                                                                        >  1:575777/1‑251 (MQ=255)
                                                                                                                                                                  gATTGCAGGTGAGCGTCACTTAGTACAACCCTTTGCTGATGAAGCAACTGCACGCAAAGTGGGAGCAAGACACGGATCGTCGGTTATTTTAACCGGCAA‑AGCACAGGAAATGGCGAAACGAGGCCTCCCTGTCTGGCAAGCGGAAAACGGGGTCTGGCTGACAGCAACGGTTGCGGTTGAGTTTCTTGAGTGGTAAGGGAATCTTCCCTCATGAAGTGACCACTTGTAGCTTAAAAAATTGCACTGTGCt                                                                       <  2:457731/250‑1 (MQ=255)
                                                                                                                                                                   aTTGCAGGTGAGCGTCACTACGTACATCTCTCTGCTGATGAAGCAACGGCACGCAAAGTGGGAGCAAGACACGGATCGTCGGTTATTTTAACCGTCAA‑AGCACAGGAAATGGCGAAACGAGGCCTCCCTTTCTGGCAAGCGGAAAACGGGGTCTGGCTGACATCCACGGTTGCGGGTGAGTTTCTTGAGTGGTAAGGGACACTTCCCTCATGAAGTGACCACTTGTAGCTTACATAATTGCACTGTGCTg                                                                      >  2:82802/1‑250 (MQ=255)
                                                                                                                                                                                        gTACATCTCTCTGCTGATGAAGCAACGGCACGCAAAGTGGGAGCAAGACACGGATCGTCGGTTATTTTAACCGTCAA‑AGCACAGGAAATGGCGAAACGAGGCCTCCCTTTCTGGCAAGCGGAAAACGGGGTCTGGCTGACATCAACGGTTGCGGTTGAGTTTCTTGAGTGGTAAGGGAAACTTCCCTCATGAAGTGACCACTTGTAGCTTAAATAATTGCACTGTGCTGATAAGCGCAGCAATAAACAAcg                                                >  1:88014/1‑251 (MQ=255)
                                                                                                                                                                                                cactgctgAGGAAGCAACGGCACGCAAAGTGGGAGCAAGACACTGATCATCTGTTATGTTAACCGGCAAGAG‑ACAGGAAATGGCGAAACGAGGCCTCCCTTTCGGGCAAGCGGAAAACGTGGTCTGGCTTACATCATCGGTTGCGGTTGAGTTTCTTGAGTGGTAAGGGAAACTTCCCTCATGAAGTGACCACTTGTAGCTTAAATAAGTGCGCTGTGCAGATAAGTGCAGCAATAAACTAcgccgccg                                          <  2:410203/247‑1 (MQ=255)
                                                                                                                                                                                                          gaAGCAACGGCACGCAAAGTGGGGGCAAGACACGGATCGTCGGTTATTTTAACCGTCAA‑AGCACAGGAAATGGCGAAACGAGGCCTCCCTTTCTGGCAAGCGGAAAACGGGGTCTGGCTGACATCAACGGTTGCGGTTGAGTTTCTTGAGTGGTAAGGGAAACTTCCCTCATGAAGTGACCACTTGTAGCTTAAATAATTGCACTGTGCTGATAAGCGCAGCAATAAACAACGCCTCCGCGCCGCagag                                >  2:10509/1‑249 (MQ=255)
                                                                                                                                                                                                          gaAGCAACGGCACGCAAAGTGGGAGCAAGACACGGATCGTCGGTTATTTTAACCGTCAA‑AGCACAGGAAATGGCGAAACGAGGCCTCCCTTTCTGGCAAGCGGAAAACGGGGTCTGGCTGACATCAACGGTTGCGGTTGAGTTTCTTGAGTGGTAAGGGAAACTTCCCTCATGAAGTGACCACTTGTAGCTTAAATAAtt                                                                                 <  1:430825/200‑1 (MQ=255)
                                                                                                                                                                                                          gaAGCAACGGCACGCAAAGTGGGAGCAAGACACGGATCGTCGGTTATTTTAACCGTCAA‑AGCACAGGAAATGGCGAAACGAGGCATCCCTTTCTGGCAAGCGGAAAACGGGGTCTGGCTGACATCAACGGTTGCGGTTGAGTTTCTTGAGTGGTAAGGGAAACTTCCCTCATGAAGTGACCACTTGTATCTTCAATAAtt                                                                                 >  2:430825/1‑200 (MQ=255)
                                                                                                                                                                                                              caccGGCACGCAAAGTGGGAGCAAGACACGGATCGTCGGTTATTTTAACCGTCAA‑AGCACAGGAAATGGCGAAACGAGGCCTCCCTTTCTGGCAAGCGGAAAACGGGGTCTGGCTGACATCAACGGTTGCGGTTGAGTTTCTTGAGTGGTAAGGGAAACTTCCCTCATGAAGTGACCACTTGTAGCTTAAATAATTGCACTGTGCTGATAAGCGCAGCAATAAACAACGCCGCCGCGCCGCagagaggagc                          >  2:136906/4‑247 (MQ=255)
                                                                                                                                                                                                                aCGGCACGCAAAGTGGGAGCAAGACACGGATCGTCGGTTATTTTAACCGTCAA‑AGCACAGGAAATGGCGAAACGAGGCCTCCCTTTCTGGCAAGCGGAAAACGGGGTCTGGCTGACATCAACGGTTGCGGTTGAGCTTATTGAGTGGTAAGGGAAACTTCCCTCATGAAGTGACCACTTGGAGCTTAAATAATTGCACTGTGCTGATAAGCGCAGCAATAAACAACGCCGCCGCGCCGCCGAGGGTGgcgc                        >  2:530795/1‑251 (MQ=255)
                                                                                                                                                                                                                  ggCACGCAAAGTGGGAGCAAGACACGGATCGTCGGTTATTTTAACCGGCAA‑AGCACAGGAAATGGCGAAACGAGGCCTCCCTTTCTGGCAAGCGGAAAACGGGGTCTGGCTGACATCAACGGTTGCGGTTGAGTTTCTTGAGTGGTAAGGGAAACTTCCCTCATGAAGTGACCACTTGTAGCTTAAATAATTGCACTGTGCTGATAAGCGCAGCAATAAACAACGCCGCCGCGCCGCAGAGAGTGGCGCTg                      <  2:408191/251‑1 (MQ=255)
                                                                                                                                                                                                                   gCACGCAAAGTGGGAGCAAGACACGGATCGTCGGTTATTTTAACCGTCAA‑AGCACAGTAAATGGCGAAACGAGGCCTCCCTTTCTGGCAAGCGGAAAACGGGGTCTGGCTGACATCAACTGTTGCGGTTGAGTTTCTTGAGTGGTAAGGGAAACTTCCCTCATGAAGTTACCACTTGTAGCTTAaata                                                                                    <  2:313143/188‑1 (MQ=255)
                                                                                                                                                                                                                   gCACGCAAAGTGGGAGCAAGACACGGATCGTCGGTTATTTTAACCGTCAA‑AGCACAGGAAATGGCGAAACGAGGCCTCCCTTTCTGGCAAGCGGAAAACGGGGTCTGGCTGACATCAACGGTTGCGGTTGAGTTTCTTGAGTGGTAAGGGAAACTTCCCTCATGAAGTGACCACTTGTAGCTTAaata                                                                                    >  1:313143/1‑188 (MQ=255)
                                                                                                                                                                                                                         aaaGTGGGAGCAAGACACGGATCGTCGGTTATTTTAACCGTCAA‑AGCACAGGAAATGGCGAAACGAGGCCTCCCTTTCTGGCAAGCGGAAAACGGGGTCTGGCTGACATCAACGGTTGCGGTTGGGTTTCTTGAGTGGTAAGGGAAACTTCCCTCATGAAGTGACCACTTGTAGCTTAAATAATTGAACTGTGCTGATAAGCGCAGCAATAAACAACGCCGCCGCGCCGCAGAGAGTGGCGCTGGCAAgcg               >  2:259507/1‑251 (MQ=255)
                                                                                                                                                                                                                               ggAGCAAGACACGGATCGTCGGTTATTTTAACCGTCAA‑AGCACAGGAAATGGCGAAACGAGGCCTCCCTTTCTGGCAAGCGGAAAACGGGGTCTGGCTGACATCAACGGTTGCGGTTGAGTTTCTTGAGTGGTAAGGGAAACTTCCCTCATGAAGTGACCACTTGTAGCTTAAATAATTGCACTGTGCTGATAAGCGCAGCAATAAACAACGCCGCCGCGCCGCAGAGAGTGGCGCTGGCAAGCGCCGATg         >  1:645685/1‑251 (MQ=255)
                                                                                                                                                                                                                               ggAGCAAGACACGGATCGTCGGTTATTTTAACCGTCAA‑AGCACAGGAAATGGCGAAACGAGGCCTCCCTTTCTGGCAAGCGGAAAACGGGGTCTGGCTGACATCAACGGTTGCGGTTGAGTTTCTTGAGTGGTAAGGGAAACTTCCCTCATGAAGTGACCACTTGTAGCTTAAATAATTGCACTGTGCTGATAAGCGCAGCAATAAACAACGCCGCCGCGCCGCAGAGAGTGGCGCTGGCAAGCGCCGATg         <  1:620936/251‑1 (MQ=255)
                                                                                                                                                                                                                                 aGCAAGACACGGATCGTCGGTTATTTTAACCGTCAA‑AGCACAGGAAATGGCGAAACGAGGCCTCCCTTTCTGGCAAGCGGAAAACGGGGTCTGGCTGACATCAACGGTTGCGGTTGAGTTTCTTGAGTGGTAAGGGAAACTTCCCTCATGAAGTGACCACTTGTAGCTTAAATAATTGCACTGTGCTGATAAGCGCAGCAATAAACAACGCCGCCGCGCCGCAGAGAGTGGCGCTGGCAAGCGCCGATGtt       >  1:244299/1‑251 (MQ=255)
                                                                                                                                                                                                                                      gACACTGATCGGTGGTTATTTTAACCGGCAA‑AGCACAGGAAATGGCGAAACGAGGCCTCCCTGTCTGGCAAGCGGAAAACGGGGTCTGGCTGACAGCAACGGTTGCGGTTGAGTTTCTTGAGTGGTAAGGGAAACTTCCCTCAAGAAGTGACCACTTGTAGCTTAAATAAGTGCACTGTGCTGATAAGCTCAGCAATAAACAACGCCGCCGCGCCGCAGAGAGTGGCGCTGGCAAGCGCCGATGTTCCAtt  <  2:124640/251‑1 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                
AAAAAGCGTTTTAGTTATTCCAGTGATGGTAGGTGCATCCGCGCGGTTCAGGGGCATTCGACTTCGCAGGTCGCTATTTCCTTTGCCGAAAAAA‑CACCGCCACAGTTTCTCTACCACGGTACGGCAAGCCGTTTCCTTGACGAGATAAAAAAACAGGGGCTGATTGCAGGTGAGCGTCACTACGTACATCTCTCTGCTGATGAAGCAACGGCACGCAAAGTGGGAGCAAGACACGGATCGCCGGTTATTTTAACCGTCAA‑AGCACAGGAAATGGCGAAACGAGGCCTCCCTTTCTGGCAAGCGGAAAACGGGGTCTGGCTGACATCAACGGTTGCGGTTGAGTTTCTTGAGTGGTAAGGGAAACTTCCCTCATGAAGTGACCACTTGTAGCTTAAATAATTGCACTGTGCTGATAAGCGCAGCAATAAACAACGCCGCCGCGCCGCAGAGAGTGGCGCTGGCAAGCGCCGATGTTCCATT  >  NC_000913/4561128‑4561607

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