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

AAAAAGCGTTTTAGTTATTCCAGTGATGGTAGGTGCATCCGCGCGGTTCAGGGGCATTCGACTTCGCAGGTCGCTATTTCCTTTGCCGAAAAAACACCGCCACAGTTTCTCTACCACGGTACGGCAAGCCGTTTCCTTGACGAGATAAAAAAACAGGGGCTGATTGCAGGTGAGCGTCACTACGTACATCTCTCTGCTGATGAAGCAACGGCACGCAAAGTGGGAGCAAGACACGGATCGCCGGTTATTTTAACCGTCAAAGCACAGGAAATGGCGAAACGAGGCCTCCCTTTCTGGCAAGCGGAAAACGGGGTCTGGCTGACATCAACGGTTGCGGTTGAGTTTCTTGAGTGGTAAGGGAAACTTCCCTCATGAAGTGACCACTTGTAGCTTAAATAATTGCACTGTGCTGATAAGCGCAGCAATAAACAACGCCGCCGCGCCGCAGAGAGTGGCGCTGGCAAGCGCCGATGTTCCATT  >  NC_000913/4561128‑4561607
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                
aaaaaGCGTTTTGGTTATTCCAGTGATGGTCGGTGCATCCGCGCGGTTCAGGGGCATTCGACTTCGAAGGTCGCTATTACCTTTGCCGAAAACACACCGCCACAGTTTCTCTACCACGGTACGGCAAGCCGTTTCATTGACGAGATAAAAAAACAGGGGCTGATTGCAGGTGAGCGTCACTACGTACATCTCTCTGCTGATGAAGCAACGGCACGAACAGAGGGAGCAAGACCCGGCTCGTCGGTTAtttt                                                                                                                                                                                                                                        >  2:565341/1‑251 (MQ=255)
      cGTTTTAGTTATTCCAGTGATGGTAGGTGCATCCGCGCGGTTCAGGGGCATTCGACTTCGCAGGTCGCTATTTCCTTTGCCGAAAAAACACCGCCACAGTTTCTCTACCACGGTACGGCAAGCCGTTTCCTTGACGAGATAAAAAAACAGGGGCTAATTGCAGGTGAGCGTCACTACGTACATCTCTCTGCTGATGAAGCAACGGCACGCAAAGTGGGAGCAAGACACGGATCGTCGGTTATTTTAACCGt                                                                                                                                                                                                                                  >  1:70373/1‑251 (MQ=255)
       gTTTTAGTTATTCCAGAGATGGTAGGTGCATCCGCGCGGTTCAGGGGCATTCGACTTCGCAGGTCGCTATTTCCTTTGCCGAAAAAACACCGCCACAGTTTCTCTACCACGGTACGGCAAGCCGTTTCCTTGACGAGATAAAAAAACAGGGGCTGATTGCAGGTGAGCGTCACTACGTACATCTCTCTGCTGATGAAGCAACGGCACGCAAAGTGGGAGCAAGACACGGATCGTCGGTTATTTTAACCGTc                                                                                                                                                                                                                                 >  1:416316/1‑251 (MQ=255)
                aTTCCAGTGATGGTAGGTGCATCCGCGCGGTTCAGGGGCCTTCGACTTCGCAGGACGCTATTTCCTTTGCCGAAAAAACACCGCCACAGTTTCTCTACCACGGTACGGCAAGCCGTTTCCTTGACGAGATAAAAAAACAGGGGCTGCTTGCAGGTGAGCGTCCCTACGTACATCTCTCTGATGAGGACGCAACGTCACGCAAAGTGGGGGCAAGACACGGATCGTCGGGTATTTTAAACGTCaaatccacg                                                                                                                                                                                                                       >  2:55794/1‑245 (MQ=255)
                 ttCCAGTGATGGGAGGTGCATCCGCGCGGTTCAGGGGCATTCGACTTCGCAGGTCGCTATTTCCTTTGCCGAAAAAACACCGCCACAGTTTCTCTACCACGGTACGGCAAGCCGTTTCCTTGACGAGATAAAAAAACAGGGGCTGATTGCAGGCGAGCGTCACTACGTACATCTCTCTGCGGATGAAGCAACGGCACGCAAAGTGGGAGCAAGCCACGGATCGTCGGTTATTTTACCCGTACAAGCACAgg                                                                                                                                                                                                                      >  2:621814/1‑251 (MQ=255)
                  cccggtgttgtgagcttgaTCCGCGCGGTCCAGGGCAATTCGCTTTCGCTGGTCGTTATCCCCTTTGCGAAAAAATCCCCGCCACAGTTTCTCTACCACGGTACGGCAAGCCGTTTCCGTGACGAGATAAAAAAACAGGGGCTGATTGCAGGTGAGCGTCACTACGTACATCTCTCTGCTGATGAAGCAACGGCACGCAAAGTGGGAGCAAGACACGTATCGTCGGTTATTTTAACCGTCAAAGcaca                                                                                                                                                                                                                        <  1:359748/230‑1 (MQ=255)
                   ccAGTGATGGTAGGTGCATCCGCGCGGTTCCGGGGCATTCGACTTCGCAGGTCGCTATTTCCTTTGCCGAAAAAACACCGCCACAGTTTCTCTACCACGGTACGGCAAGCCGTTTCCTTGACGAGATAAAAAAACAGGGGCTGATTGCAGGTGAGCGTCACTACGTACATCTCTCTGCTGATGAAGCAACGGCACGCAAAGTGGGAGCAAGACACGGATCGTCGGTTATTTTAACCGTCAAAGCACAGGa                                                                                                                                                                                                                     >  2:119948/1‑250 (MQ=255)
                        gATGGTAGGGGCATCCGCGCGGTTCAGGGGCATTCGACTTCGCAGGTCGCTATTTCATTTGCCGAAAAAACACCGCCACAGTTTCTCTACCACGGTACGGCAAGCCGTTTCCTTGACGAGATAAAAAAACAGGGGCTGATTGCAGGTGAGCGTCAATACGTACATCTCTCTGCTGTTGAAGCAACGGCACGCAAAGAGGGAGCAAGAAACGGATCGTCGGGTATTTTAACCGTCAAAGCACAGGCAATggg                                                                                                                                                                                                               >  2:625116/1‑250 (MQ=255)
                               ggTGCATCCGCGCGGTTCAGGGGCATTAGACTTCGCAGGTCGCTATTTCCTTTGCCGAAAAAACACCGCCACAGTTTCTCTACCACGGTACGGCAAGCCGTTTCCTTGACGAGATAAAAAAACAGGGGCTGATTGCAGGTGAGAGTCACTACGTACATCTCTCTGCTGATGAAGCAACGGCACGCAAAGTGGGAGCAAGACACGGATCGTCGGTTATTTTAACCGTCAAAGCACAGGAAATGGAGAAACg                                                                                                                                                                                                         >  2:645096/1‑250 (MQ=255)
                                gTGCATCCGCGCGGTTCAGGGGCATTCGACTTCGCAGGTCGCTATTTCCTTTGCCGAAAAAACACCGCCACAGTTTCTCTACCACGGTACGGCAAGCCGTTTCCTTGACGAGATAAAAAAACAGGGGCTGATTGCAGGTGAGCGTCACTACGTACATCTCTCTGCTGATGAAGCAACGGCACGCAAAGTGGGAGCAAGACACGGATCGTCGGTTATTTTAACCGTCAAAGCACAGGAAATGGCGAAACg                                                                                                                                                                                                         >  1:379237/1‑249 (MQ=255)
                                   cATCCGCGCGGTTCAGGGGCATTCGACTTCGCAGGTCGCTATTTCCTTTGCCGAAAAAACACCGCCACAGTTTCTCTACCACGGTACGGCAAGCCGTTTCCTTGACGAGATAAAAAAACAGGGGCTGATTGCAGGTGAGCGTCACTACGTACATCTCTCTGCTGATGAAGCAACGGCACGCAAAGTGGGAGCAAGACACGGATCGTCGGTTATTTTAACCGTCAAAGCACAGGAAATGGCGAAACGAGGcc                                                                                                                                                                                                    >  1:408191/1‑251 (MQ=255)
                                    aTCCGCGCGGTTCAGGGGCATTCGACTTCGCAGGTCGCTATTTCCTTTGCCGAAAAAACACCGCCACAGTTTCTCTACCACGGTACGGCAAGCCGTTTCCTTGACGAGATAAAAAAACAGGGGCTGATTGCAGGTGAGCGTCACTACGTACATCTCTCTGCTGATGAAGCAACGGCACGCAAAGTGGGAGCAAGACACGGATCGTCGGTTATTTTAACCGTCAAAGCACAGGCAATGGCGAAACGAGGCCt                                                                                                                                                                                                   >  1:22038/1‑251 (MQ=255)
                                          gcgGTTCAGGGGCATTCGACTTCGCAGGTCGCTATTTCCTTTGCCGAAAAAACACCGCCACAGTTTCTCTACCACGGTACGGCAAGCCGTTTCCTTGACGAGATAAAAAAACAGGGGCTGATTGCAGGTGAGCGTCACTACGTACATCTCTCTGCTGATGAAGCAACGGCACGCAAAGTGGGAGCAAGACACGGATCGTCGGTTATTTTAACCGTCAAAGCACAGGAAATGGCGAAACGAGGCCTCCCttt                                                                                                                                                                                             >  1:298627/1‑251 (MQ=255)
                                                           gACTTCGCAGGTCGCTATTTCCTTTGCCGAAAAAACACCGCCACAGTTTCTCTACCACGGTACGGCAAGCCGTTTCCTTGACGAGATAACAAAACAGGGGCTGATTGCAGGTGAGCGTCACTACGTACATCTCTCTGCTGATGAAGCAACGGCACGCAAAGTGGGAGCAAGACACGGATCGTCGGTTATTTTAACCGTCAGAGCACAGGAAATAGCGCAACGAGGCCTCCCTTTCTGGAAAGAGGAAAACg                                                                                                                                                                            >  2:8248/1‑251 (MQ=255)
                                                            aCTTCGCAGGTCGCTATTTCCTTTGCCGAAAAAACACCGCAACAGTTTCTCTACCACGGTACGGCAAGCCGTTTCCTTGACGAGATAAAAAAACAGGGGCTGATTGCAGGTGAGCGTCACTACGTACATCTCTCTGCTGATGAAGCAACGGCACGCAAAGTGGGAGCAAGACACGGATCGTCGGTTATTTTAACCGTCAAAGCACAGGAAATGGCGAAACGAGGCCTCCCTTTCTGGCAAGCGGAAACCgg                                                                                                                                                                           >  2:351560/1‑251 (MQ=255)
                                                                     gTCGCTATTTCCTTTGCCGAAAAAACACCGCCACAGTTTCTCTACCACGGTACGGCAAGCCGTTTCCTTGACGAGATAAAAAAACAGGGGCTGATTGCAGGTGAGCGTCACTACGTACATCTCTCTGCTGATGAAGCAACGGCACGCAAAGTGGGAGCAAGACACGGATCGTCGGtt                                                                                                                                                                                                                                             >  1:375628/1‑177 (MQ=255)
                                                                     gTCGCTATTTCCTTTGCCGAAAAAACACCGCCACAGTTTCTCTACCACGGTACGGCAAGCCGTTTCCTTGACGAGATAAAAAAACAGGGGCTGATTGCAGGTGAGCGTCACTACGTACATCTCTCTGCTGATGAAGCAACGGCACGCAAAGTGGGAGCAAGACACGGATCGTCGGTTATTTTAACCGTCAAAGCACAGGAAATGGCGAAACGAGGCCTCCCTTTCTGGCACGCGGAAAACGGGGTCTGGCt                                                                                                                                                                  >  2:199129/1‑251 (MQ=255)
                                                                          tATTTCCTTTGCCGAAAAAACACCGCCACAGTTTCTCTACCACGGTACGGCAAGCCGTTTCCTTGACGAGATAAAAAAACAGGGGCTGATTGCAGGTGAGCGTCACTACGTACATCTCTCTGCTGATGAAGCAACGGCACGCAAAGTGGGAGCAAGACACGGATCGTCGGTTATTTTAACCGTCAAAGCACAGGAAATGGCGAAACGAGGCCTCCCTTTCTGGCAAGCGGAAAACGGGGTCTGGCTGACAt                                                                                                                                                             >  1:359866/1‑251 (MQ=255)
                                                                           tttcctgttgacaaaaaaaCCCCGCCACAGTTTCTCTACCACGGTACGGCAAGCCGTTTCCTTGACGAGATAAAAAAACAGGGGCTGATTGCAGGTGAGCGTCACTACGTACATCTCTCTGCTGATGAAGCAACGGCACGCAAAGTGGGAGCAAGACACGGATCGTCGGTTATTTTAACCGTCAAAGCACAGGAAATGGCGAAACGAGGCCTCACTTTCTGGCAAGCGGAAAACGGGGTCTGGCTg                                                                                                                                                                 <  2:30954/233‑1 (MQ=255)
                                                                                tgtggcaaaaaaaaTACCGCCACAGTTTCACTACAACGGTACGGCAAGCCGTTTCCTTGACGAGATAAAAAAACTGGGGCTGATTGCCGGTGAGCGTCACTACGTACATCTCTCTGCTGATGAAGCTACGGCACGCAAAGTGGGAGCAAGACACGGATCGTCGGTTATTTTAACCGTCAAAGCACAGGAAATGGCGAAACGAGGCCGCCCTTTCTGGCAAGCGGAAAACGGGGTCTGGCTGACATCaa                                                                                                                                                          <  2:499975/240‑1 (MQ=255)
                                                                                    cccaaaaaacCACCGCCACAGTTTCTCTACCACGGTACGGCAAGCCGTTTCCTTGACGAGATAAAAAAACAGGGGCTGATTGCAGGTGAGCGTCACTACGTACATCTCTCTGCTGATGAAGCAACGGCACGCAAAGTGGGAGCAAGACACGGATCGTCGGTTATTTTAACCGTCAAAGCACAGGAAATGGCGAAACGAGGCCTCCCTTTCTGGCAAGCGGAAAACGGGGTCTGGCTGACATCAAcggttgc                                                                                                                                                   <  1:316794/247‑1 (MQ=255)
                                                                                     caaaaaaaaCACCGCCACAGTTTCTCTACCACGGTACGGCAAGCCGTTTCCTTGACGAGATAAAAAAACAGGGGCTGATTGCAGGTGAGCGTCACTACGTACATCTCTCTGCTGATGAAGCAACGGCACGCAAAGTGGGAGCAAGACACGGATCGGCGGTTATTTTAACCGTCAAAGCACAGGAAATGGCGAAACGAGGCCTCCCTTTCTGGCAAGCGGAAAACGGGGTCTGGCTGACATCAAcggttgcg                                                                                                                                                  <  2:70373/248‑1 (MQ=255)
                                                                                             acacCGCCACAGTTTCTCTACCACGGTACGGCAAGCCGTTTCCTTGACGAGATAAAAAAACAGGGGCTGATTGCAGGTGAGCGTCACTACGTACATCTCTCTGCTGATGAAGCAACGGCACGCAAAGTGGGAGCAAGACACGGATCGTCGGTTATTTTAACCGTCAAAGCAAAGGAAATGGCGAAACAAGGCCTCCCTTTCTGGCAAGCGGAAAACGGGGGCTGGCTGCCATCAACGGATGCGGTTGAGtt                                                                                                                                          >  2:655916/1‑251 (MQ=255)
                                                                                               ccccGCCACAGTTTCTCTACCACGGTACGGCAAGACGTTTCCTTGACGAGATTAAAAAACAGGGGCTGATGGCAGGTGAGCGTCACTACGTCCATCTCTCTGCTGATGAAGCAACGGCACGCAAAGTGGGAGCAAGACACGGATCGTCGGTTATTTTAACCGTCAAAGCACAGGAAATGGCGAAACGAGGCCTCCCTTTCTGGCAAGCTGAAAACGGGGTCTGGCTGACATCAACGGTTGCGGTTGAGttt                                                                                                                                         <  2:597725/249‑1 (MQ=255)
                                                                                                   gCCACAGTTTCTCTACCACGGTACGGCAAGCCGTTTCCTTGACGAGATAAAAAAACAGGGGCTGATTGCAGGTGAGCGTCACTACGTACATCTCTCTGCTGATGAAGCAACGGCACGCAAAGTGGGAGCAAGACACGGATCGTCGGTTATTTTAACCGTCAAAGCACAGGAAATGGCGAAACGAGGCCTCCCTTTCTGGCAAGCGGAAAACGGGGTCTGGCTGACATCAACGGTTGCGGTTGGGTTTCTTg                                                                                                                                     >  1:22122/1‑251 (MQ=255)
                                                                                                        aGTTTCTCTACCACGGTACGGCAAGCCGTTTCCTTGACGAGATAAAAAAACAGGGGCTGATTGCAGGTGAGCGTCACTACGTACATCTCTCTGCTGATGAAGCAACGGCACGCAAAGTGGGAGCAAGACACGGATCGTCGGTTATTTTAACCGTCAAAGCACAGGAAATGGCGAAACGAGGCCTCCCTTTCTGGCAAGCGGAAAACGGGGTCTGGCTGACATCAACGGTTGCGGTTGAGTTTCTTGAGTgg                                                                                                                                <  2:582451/251‑1 (MQ=255)
                                                                                                              tctACCACGGTACGGCAAGCCGTTTCCTTGACGAGATAAAAAAACAGGGGCTGATTGCAGGTGAGCGTCACTACGTACATCTCTCTGCTGATGAAGCAACGGCACGCAAAGTGGGAGCAAGACACGGATCGTCGGTTATTTTAACCGTCAAAGCACAGGAAATGGCGAAACGAGGCCTCCCTTTCTGGCAAGCGGAAAACGGGGTCTGGCTGACATCAACGGTTGCGGTTGGGTTTCTTGAGTGGTAAggg                                                                                                                          >  2:294430/1‑251 (MQ=255)
                                                                                                                       gttcGGCAAGCCGGTTCCTTGACGAGATAAAAAAACAGGGGCTGATTGCAGGTGAGCGTCACTACGTACATCTCTCTGCTGATGAAGCAACGGCACGCAAAGTGGGAGCAAGACACGGATCGTCGGTTATTTTAACCGTCAAAGCACAGGAAATGGTGAAACGAGGCCTCCCTTTCTGGCAAGCGGAAAACGGTGTCTGGCTTACATCAACTGTTGCGGTTGAGTTTCTTGAGTGGTAAGGGAAACTTccc                                                                                                                 <  1:310976/248‑1 (MQ=255)
                                                                                                                        tACGGCAAGCCGTTTCCTTGACGAGATAAAAAAACAGGGGCTGATTGCAGGTGAGCGTCACTACGTACATCTCTCTGCTGATGAAGCAACGGCACGCAAAGTGGGAGCAAGACACGGATCGTCGGTTATTTTAACCGTCAAAGCACAGGAAATGGCGAAACGAGGCCTCCCTTTCTGGCAAGCGGAAAACGGGGTCTGGCTGACATCAACGGTTGCGGTTGAGTTTCTTGAGGGGTAAGGGAAACTTCCCt                                                                                                                >  2:620936/1‑251 (MQ=255)
                                                                                                                                   gTTTCCTTGACGAGATAAAAAAACAGGGGCTGATTGCAGGTGAGCGTCACTACGTACATCTCTCTGCTGATGAAGCAACGGCACGCAAAGTGGGAGCAAGACACGGATCGTCGGTTATTTTAACCGTCAAAGCACAGGAAATGGCGAAACGAGGCCTCCCTTTCTGGCAAGCGGAAAACGGGGt                                                                                                                                                                        <  2:328066/184‑1 (MQ=255)
                                                                                                                                      tccgtGACGAGATAAAAAAACAGGGGCTGATTGCATGTGAGCGTCACTACGTACATCTCTCTGCTGATGAAGCAACGGCACGCAAAGTGGGAGCAAGACACGGATCGTCGGTTATTTTAACCGTCAAAGCACAGGAAATGGCGAAACGAGGCCTCCCTTGCTCGCAAGCGGAAAAAGGGGTCTGGCt                                                                                                                                                                  <  2:341364/183‑1 (MQ=255)
                                                                                                                                      tCCTTGACGAGATAAAAAAACAGGGGCTGATTGCAGGTGAGCGTCACTACGTACATCTCTCTGCTGATGAAGCAACGGCACGCACAGTGGGAGCAAGACACGGATCGTCGGTTATTTTAACCGTCAAAGCACAGGAAATGGCGAAACGAGGCCTCCCTTTCTGGCAAGCGGAAAACGGGGTCTGGCt                                                                                                                                                                  >  1:341364/1‑187 (MQ=255)
                                                                                                                                       ccTTGACGAGATAAAAAAACAGGGGCTGATTGCAGGTGAGCGTCACTACGTACATCTCTCTGCTGATGAAGCAACGGCACGCAAAGTGGGAGCAAGACACGGATCGTCGGTTATTTTAACCGTCAAAGCACAGGAAATGGCGAAACGAGGCCTCCCTTTCTGGCAAGCGGAAAACGGGGTCTGGCTGACATCAACGGTTGCGGTTGAGTTTCTTGAGTGGTAAGGGAAACTTCCCTCATGAAGTGACCACt                                                                                                 <  1:310442/251‑1 (MQ=255)
                                                                                                                                          tGACGAGATAAAAAAACAGGGGCTGATTGCAGGTGAGCGTCACTACGTACATCTCTCTGCTGATGAAGCAACGGCACGCAAAGTGGGAGCAAGACACGGATCGTCGGTTATTTTAACCGTCAAAGCACAGGAAATGGCGAAACGAGGCCTCCCTTTCTGGCAAGCGGAAAACGGGGTCTGGCTGACATCAACGGTTGCGGTTGAGTTTCTTGAGTGGTAAGGGAAACTTCCCTCATGAAGTGACCACtt                                                                                                <  1:554777/249‑1 (MQ=255)
                                                                                                                                           gACGAGATAAAAAAACAGGGGCTGATTGCAGGTGAGCGTCACTACGTACATCTCTCTGCTGATGAAGCAACGGCACGCAAAGTGGGAGCAAGACACGGATCGTCGGTTATTTTAACCGTCAAAGCACAGGAAATGGCGAAACGAGGCCTCCCTTTCTGGCAAGCGGAAAACGGGGTCTGGCTGACATCAACGGTTGCGGTTGAGTTTCTTGAGTGGTAAGGGAAACTTCCCTCATGAAGTGACCACTTGt                                                                                              <  1:294430/250‑1 (MQ=255)
                                                                                                                                            aCGAGATAAAAAAACAGGGGCTGATTGCAGGTGAGCGTCACTACGTACATCTCTCTGCTGATGAAGCAACGGCACGCAAAGTGGGAGCAAGACACGGATCGTCGGTTATTTTAACCGTCAAAGCACAGGAAATGGCGAAACGAGGCCTCCCTTTCTGGCAAGCGGAAAACGGGGTCTGGCTGACATCAACGGTTGCGGTTGAGTTTCTTGAGTGGTAAGGGAAACTTCCCTCATGAAGTGACCACTTGTAg                                                                                            <  1:407510/251‑1 (MQ=255)
                                                                                                                                                    aaaaaaCAGGGGCTGATTGCAGGTGAGCGTCACTACGTACATCTCTCTGCTGATGAAGCAACGGCACGCAAAGTGGGAGCAAGACACGGATCGTCGGTTATTTTAACCGTCAAAGCACAGGAAATGGCGAAACGAGGCCTCCCTTTCTGGCAAGCGGAAAACGGGGTCTGGCTGACATCAACGGTTGCGGTTGAGTTTCTTGAGTGGTAAGGGAAACTTCCCTCATGAAGTGACCACTTGTAGCTTAaata                                                                                    <  1:199129/251‑1 (MQ=255)
                                                                                                                                                     caaaaCATGGGCTAATTGCAGGTGAGCGTCACTACGTACATCTCTCTGCTGATGAAGCAACGGCACGCAAAGTGGGAGCAAGACACGGATCGTCGGGTATTTTAACCGTCAAAGCACAGGCAATGGCGCAACGAGGCCTCCCTTTCTGGCAAGCGGAAAACGGGGTCTGGCTGAAATCAACGGTTGAGGTTGGGTTTCTTGGGTGGGAAAGGAAGCTTCCCTACTGAAGAGACAACTTgaggcttacata                                                                                    >  2:380198/2‑239 (MQ=255)
                                                                                                                                                        aaCAGGGGCTGATTGCAGGTGAGCGTCACTACGTACATCTCTCTGCTGATGAAGCAACGGCACGCAAAGTGGGAGCAAGACACGGATCGTCGGTTATTTTAACCGTCAAAGCACAGGAAATGGCGAAACGAGGCCTCCCTTTCTGGCAAGCGGAAAACGGGGTCTGGCTGACATCAACGGTTGCGGTTGAGTTTCTTGAGTGGTAAGGGAAACTTCCCTCATGAAGTGACCACTTGTAGCTTAAATAATTg                                                                                <  1:645096/251‑1 (MQ=255)
                                                                                                                                                             ggtcTGAGTGCAGGTTAGCGTCACTACCTTCTACTCTCTGCTGATGAAGCACCGGCACGCAAAGTGGGAGCAAGACACGGATCTTCTGTTATTTTAACCGTCAAAGCACAGGAAATGTCGAATCGAGGCTTCCCTTTCTTGCAAGCTGAAAACGGGGTCTGGCTGACATCAACTGTTGCGGTTGAGTTTCTTGAGTGGTAAAGGAAACTTCCCTCATGAAGTGACCACTTTTAGCTGAAAGAATTGCACtg                                                                           <  2:593332/248‑1 (MQ=255)
                                                                                                                                                                cTGATTGCAGGTGAGCGTCACTACGTACATCTCTCTGCTGATGAAGCAACGGCACGCAAAGTGGGAGCAAGACACGGATCGTCGGATATTTTAACCGTCAAAGCACAGGAAATGGCGAAACGAGGCCTCCCTTTCTGGCAAGCGGAAAACGGGGTCTGGCTGACATCAACGGTTGCGGTTGAGTTTCTTGAGTGGTAAGGGAAACTTCCCTCATGAAGTGACCACTTGTAGCTTAAATAATTGCACTGTGc                                                                        >  1:575777/1‑251 (MQ=255)
                                                                                                                                                                  gATTGCAGGTGAGCGTCACTTAGTACAACCCTTTGCTGATGAAGCAACTGCACGCAAAGTGGGAGCAAGACACGGATCGTCGGTTATTTTAACCGGCAAAGCACAGGAAATGGCGAAACGAGGCCTCCCTGTCTGGCAAGCGGAAAACGGGGTCTGGCTGACAGCAACGGTTGCGGTTGAGTTTCTTGAGTGGTAAGGGAATCTTCCCTCATGAAGTGACCACTTGTAGCTTAAAAAATTGCACTGTGCt                                                                       <  2:457731/250‑1 (MQ=255)
                                                                                                                                                                   aTTGCAGGTGAGCGTCACTACGTACATCTCTCTGCTGATGAAGCAACGGCACGCAAAGTGGGAGCAAGACACGGATCGTCGGTTATTTTAACCGTCAAAGCACAGGAAATGGCGAAACGAGGCCTCCCTTTCTGGCAAGCGGAAAACGGGGTCTGGCTGACATCCACGGTTGCGGGTGAGTTTCTTGAGTGGTAAGGGACACTTCCCTCATGAAGTGACCACTTGTAGCTTACATAATTGCACTGTGCTg                                                                      >  2:82802/1‑250 (MQ=255)
                                                                                                                                                                                        gTACATCTCTCTGCTGATGAAGCAACGGCACGCAAAGTGGGAGCAAGACACGGATCGTCGGTTATTTTAACCGTCAAAGCACAGGAAATGGCGAAACGAGGCCTCCCTTTCTGGCAAGCGGAAAACGGGGTCTGGCTGACATCAACGGTTGCGGTTGAGTTTCTTGAGTGGTAAGGGAAACTTCCCTCATGAAGTGACCACTTGTAGCTTAAATAATTGCACTGTGCTGATAAGCGCAGCAATAAACAAcg                                                >  1:88014/1‑251 (MQ=255)
                                                                                                                                                                                                cactgctgAGGAAGCAACGGCACGCAAAGTGGGAGCAAGACACTGATCATCTGTTATGTTAACCGGCAAGAGACAGGAAATGGCGAAACGAGGCCTCCCTTTCGGGCAAGCGGAAAACGTGGTCTGGCTTACATCATCGGTTGCGGTTGAGTTTCTTGAGTGGTAAGGGAAACTTCCCTCATGAAGTGACCACTTGTAGCTTAAATAAGTGCGCTGTGCAGATAAGTGCAGCAATAAACTAcgccgccg                                          <  2:410203/247‑1 (MQ=255)
                                                                                                                                                                                                          gaAGCAACGGCACGCAAAGTGGGGGCAAGACACGGATCGTCGGTTATTTTAACCGTCAAAGCACAGGAAATGGCGAAACGAGGCCTCCCTTTCTGGCAAGCGGAAAACGGGGTCTGGCTGACATCAACGGTTGCGGTTGAGTTTCTTGAGTGGTAAGGGAAACTTCCCTCATGAAGTGACCACTTGTAGCTTAAATAATTGCACTGTGCTGATAAGCGCAGCAATAAACAACGCCTCCGCGCCGCagag                                >  2:10509/1‑249 (MQ=255)
                                                                                                                                                                                                          gaAGCAACGGCACGCAAAGTGGGAGCAAGACACGGATCGTCGGTTATTTTAACCGTCAAAGCACAGGAAATGGCGAAACGAGGCCTCCCTTTCTGGCAAGCGGAAAACGGGGTCTGGCTGACATCAACGGTTGCGGTTGAGTTTCTTGAGTGGTAAGGGAAACTTCCCTCATGAAGTGACCACTTGTAGCTTAAATAAtt                                                                                 <  1:430825/200‑1 (MQ=255)
                                                                                                                                                                                                          gaAGCAACGGCACGCAAAGTGGGAGCAAGACACGGATCGTCGGTTATTTTAACCGTCAAAGCACAGGAAATGGCGAAACGAGGCATCCCTTTCTGGCAAGCGGAAAACGGGGTCTGGCTGACATCAACGGTTGCGGTTGAGTTTCTTGAGTGGTAAGGGAAACTTCCCTCATGAAGTGACCACTTGTATCTTCAATAAtt                                                                                 >  2:430825/1‑200 (MQ=255)
                                                                                                                                                                                                              caccGGCACGCAAAGTGGGAGCAAGACACGGATCGTCGGTTATTTTAACCGTCAAAGCACAGGAAATGGCGAAACGAGGCCTCCCTTTCTGGCAAGCGGAAAACGGGGTCTGGCTGACATCAACGGTTGCGGTTGAGTTTCTTGAGTGGTAAGGGAAACTTCCCTCATGAAGTGACCACTTGTAGCTTAAATAATTGCACTGTGCTGATAAGCGCAGCAATAAACAACGCCGCCGCGCCGCagagaggagc                          >  2:136906/4‑247 (MQ=255)
                                                                                                                                                                                                                aCGGCACGCAAAGTGGGAGCAAGACACGGATCGTCGGTTATTTTAACCGTCAAAGCACAGGAAATGGCGAAACGAGGCCTCCCTTTCTGGCAAGCGGAAAACGGGGTCTGGCTGACATCAACGGTTGCGGTTGAGCTTATTGAGTGGTAAGGGAAACTTCCCTCATGAAGTGACCACTTGGAGCTTAAATAATTGCACTGTGCTGATAAGCGCAGCAATAAACAACGCCGCCGCGCCGCCGAGGGTGgcgc                        >  2:530795/1‑251 (MQ=255)
                                                                                                                                                                                                                  ggCACGCAAAGTGGGAGCAAGACACGGATCGTCGGTTATTTTAACCGGCAAAGCACAGGAAATGGCGAAACGAGGCCTCCCTTTCTGGCAAGCGGAAAACGGGGTCTGGCTGACATCAACGGTTGCGGTTGAGTTTCTTGAGTGGTAAGGGAAACTTCCCTCATGAAGTGACCACTTGTAGCTTAAATAATTGCACTGTGCTGATAAGCGCAGCAATAAACAACGCCGCCGCGCCGCAGAGAGTGGCGCTg                      <  2:408191/251‑1 (MQ=255)
                                                                                                                                                                                                                   gCACGCAAAGTGGGAGCAAGACACGGATCGTCGGTTATTTTAACCGTCAAAGCACAGTAAATGGCGAAACGAGGCCTCCCTTTCTGGCAAGCGGAAAACGGGGTCTGGCTGACATCAACTGTTGCGGTTGAGTTTCTTGAGTGGTAAGGGAAACTTCCCTCATGAAGTTACCACTTGTAGCTTAaata                                                                                    <  2:313143/188‑1 (MQ=255)
                                                                                                                                                                                                                   gCACGCAAAGTGGGAGCAAGACACGGATCGTCGGTTATTTTAACCGTCAAAGCACAGGAAATGGCGAAACGAGGCCTCCCTTTCTGGCAAGCGGAAAACGGGGTCTGGCTGACATCAACGGTTGCGGTTGAGTTTCTTGAGTGGTAAGGGAAACTTCCCTCATGAAGTGACCACTTGTAGCTTAaata                                                                                    >  1:313143/1‑188 (MQ=255)
                                                                                                                                                                                                                         aaaGTGGGAGCAAGACACGGATCGTCGGTTATTTTAACCGTCAAAGCACAGGAAATGGCGAAACGAGGCCTCCCTTTCTGGCAAGCGGAAAACGGGGTCTGGCTGACATCAACGGTTGCGGTTGGGTTTCTTGAGTGGTAAGGGAAACTTCCCTCATGAAGTGACCACTTGTAGCTTAAATAATTGAACTGTGCTGATAAGCGCAGCAATAAACAACGCCGCCGCGCCGCAGAGAGTGGCGCTGGCAAgcg               >  2:259507/1‑251 (MQ=255)
                                                                                                                                                                                                                               ggAGCAAGACACGGATCGTCGGTTATTTTAACCGTCAAAGCACAGGAAATGGCGAAACGAGGCCTCCCTTTCTGGCAAGCGGAAAACGGGGTCTGGCTGACATCAACGGTTGCGGTTGAGTTTCTTGAGTGGTAAGGGAAACTTCCCTCATGAAGTGACCACTTGTAGCTTAAATAATTGCACTGTGCTGATAAGCGCAGCAATAAACAACGCCGCCGCGCCGCAGAGAGTGGCGCTGGCAAGCGCCGATg         >  1:645685/1‑251 (MQ=255)
                                                                                                                                                                                                                               ggAGCAAGACACGGATCGTCGGTTATTTTAACCGTCAAAGCACAGGAAATGGCGAAACGAGGCCTCCCTTTCTGGCAAGCGGAAAACGGGGTCTGGCTGACATCAACGGTTGCGGTTGAGTTTCTTGAGTGGTAAGGGAAACTTCCCTCATGAAGTGACCACTTGTAGCTTAAATAATTGCACTGTGCTGATAAGCGCAGCAATAAACAACGCCGCCGCGCCGCAGAGAGTGGCGCTGGCAAGCGCCGATg         <  1:620936/251‑1 (MQ=255)
                                                                                                                                                                                                                                 aGCAAGACACGGATCGTCGGTTATTTTAACCGTCAAAGCACAGGAAATGGCGAAACGAGGCCTCCCTTTCTGGCAAGCGGAAAACGGGGTCTGGCTGACATCAACGGTTGCGGTTGAGTTTCTTGAGTGGTAAGGGAAACTTCCCTCATGAAGTGACCACTTGTAGCTTAAATAATTGCACTGTGCTGATAAGCGCAGCAATAAACAACGCCGCCGCGCCGCAGAGAGTGGCGCTGGCAAGCGCCGATGtt       >  1:244299/1‑251 (MQ=255)
                                                                                                                                                                                                                                      gACACTGATCGGTGGTTATTTTAACCGGCAAAGCACAGGAAATGGCGAAACGAGGCCTCCCTGTCTGGCAAGCGGAAAACGGGGTCTGGCTGACAGCAACGGTTGCGGTTGAGTTTCTTGAGTGGTAAGGGAAACTTCCCTCAAGAAGTGACCACTTGTAGCTTAAATAAGTGCACTGTGCTGATAAGCTCAGCAATAAACAACGCCGCCGCGCCGCAGAGAGTGGCGCTGGCAAGCGCCGATGTTCCAtt  <  2:124640/251‑1 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                
AAAAAGCGTTTTAGTTATTCCAGTGATGGTAGGTGCATCCGCGCGGTTCAGGGGCATTCGACTTCGCAGGTCGCTATTTCCTTTGCCGAAAAAACACCGCCACAGTTTCTCTACCACGGTACGGCAAGCCGTTTCCTTGACGAGATAAAAAAACAGGGGCTGATTGCAGGTGAGCGTCACTACGTACATCTCTCTGCTGATGAAGCAACGGCACGCAAAGTGGGAGCAAGACACGGATCGCCGGTTATTTTAACCGTCAAAGCACAGGAAATGGCGAAACGAGGCCTCCCTTTCTGGCAAGCGGAAAACGGGGTCTGGCTGACATCAACGGTTGCGGTTGAGTTTCTTGAGTGGTAAGGGAAACTTCCCTCATGAAGTGACCACTTGTAGCTTAAATAATTGCACTGTGCTGATAAGCGCAGCAATAAACAACGCCGCCGCGCCGCAGAGAGTGGCGCTGGCAAGCGCCGATGTTCCATT  >  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: