Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,623,708 C→T A222V (GCC→GTC)  ppk → polyphosphate kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,623,7080CT80.0% 140.5 / 35.6 70A222V (GCC→GTC) ppkpolyphosphate kinase
Reads supporting (aligned to +/- strand):  ref base C (7/7);  new base T (27/29);  total (34/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.47e-01

TCAGATAAAGTGCCGCGCTTTGTGAATTTACCGCCAGAAGCGCCGCGTCGACGCAAGCCGATGATTCTTCTGGATAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAGCGCGAT  >  NC_000913/2623569‑2623844
                                                                                                                                           |                                                                                                                                        
tCAGATAAAGTGCCGCGCTTTGTGAATTTACCGCCAGAAGCGCCGCGTCGACGCAAGCCGATGATTCTTCTGGATAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGc                                                                                                                                          <  1:738698/140‑1 (MQ=255)
 cAGATAAAGTGCCGCGCTTTGTGAATTTACCGCCAGAAGCGCCGCGTCGACGCAAGCCGATGATTCTTCTGGATAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGc                                                                                                                                          <  2:618129/139‑1 (MQ=255)
 cAGATAAAGTGCCGCGCTTTGTGAATTTACCGCCAGAAGCGCCGCGTCGACGCAAGCCGATGATTCTTCTGGATAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGc                                                                                                                                          >  1:618129/1‑139 (MQ=255)
    aTAAAGTGCCGCGCTTTGTGAATTTACCGCCAGAAGCGCCGCGTCGACGCAAGCCGATGATTCTTCTGGATAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGTCTAt                                                                                                                                      >  2:923234/1‑140 (MQ=255)
     tAAAGTGCCGCGCTTTGTGAATTTACCGCCAGAAGCGCCGCGTCGACGCAAGCCGATGATTCTTCTGGATAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGTCTAtt                                                                                                                                     <  2:129620/140‑1 (MQ=255)
        aGTGCCGCGCTTTGTGAATTTACCGCCAGAAGCGCCGCGTCGACGCAAGCCGATGATTCTTCTGGATAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGTCTATTCaa                                                                                                                                  <  1:696141/140‑1 (MQ=255)
        aGTGCCGCGCTTTGTGAATTTACCGCCAGAAGCGCCGCGTCGACGCAAGCCGATGATTCTTCTGGATAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGTCTATTCaa                                                                                                                                  <  1:923234/140‑1 (MQ=255)
        aGTGCCGCGCTTTGTGAATTTACCGCCAGAAGCGCCGCGTCGACGCAAGCCGATGATTCTTCTGGATAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGTCTATTCaa                                                                                                                                  <  2:341517/140‑1 (MQ=255)
            ccGCGCTTTGTGAATTTACCGCCAGAAGCGCCGCGTCGACGCAAGCCGATGATTCTTCTGGATAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGTCTATTCAATGaa                                                                                                                              >  2:630100/1‑140 (MQ=255)
             cgcgCTTTGTGAATTTACCGCCAGAAGCGCCGCGTCGACGCAAGCCGATGATTCTTCTGGATAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGTCTATTCAATGAAg                                                                                                                             >  2:486399/1‑140 (MQ=255)
             cgcgCTTTGTGAATTTACCGCCAGAAGCGCCGCGTCGACGCAAGCCGATGATTCTTCTGGATAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAg                                                                                                                             >  1:787873/1‑140 (MQ=255)
              gcgcTTTGTGAATTTACCGCCAGAAGCGCCGCGTCGACGCAAGCCGATGATTCTTCTGGATAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGTCTATTCAATGAAGa                                                                                                                            >  1:449747/1‑140 (MQ=255)
              gcgcTTTGTGAATTTACCGCCAGAAGCGCCGCGTCGACGCAAGCCGATGATTCTTCTGGATAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGTCTATTCAATGAAGa                                                                                                                            >  1:279925/1‑140 (MQ=255)
              gcgcTTTGTGAATTTACCGCCAGAAGCGCCGCGTCGACGCAAGCCGATGATTCTTCTGGATAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGTCTATTCAATGAAGa                                                                                                                            >  1:261950/1‑140 (MQ=255)
                     gtgAATTTACCGCCAGAAGCGCCGCGTCGACGCAAGCCGATGATTCTTCTGGATAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGTTGAATGCCTATTCAATGAAGATGACCcg                                                                                                                     <  2:590115/140‑1 (MQ=255)
                       gAATTTACCGCCAGAAGCGCCGCGTCGACGCAAGCCGATGATTCTTCTGGATAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGTCTATTCAATGAAGATGACCcgcg                                                                                                                   <  2:694417/140‑1 (MQ=255)
                          tttACCGCCAGAAGCGCCGCGTCGACGCAAGCCGATGATTCTTCTGGATAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGTCTATTCAATGAAGATGACCCGCGATg                                                                                                                <  1:486399/140‑1 (MQ=255)
                          tttACCGCCAGAAGCGCCGCGTCGACGCAAGCCGATGATTCTTCTGGATAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGTCTATTCAATGAAGATGACCCGCGATg                                                                                                                <  2:279925/140‑1 (MQ=255)
                            tACCGCCAGAAGCGCCGCGTCGACGCAAGCCGATGATTCTTCTGGATAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGTCTATTCAATGAAGATGACCCGCGATGcc                                                                                                              >  2:183469/1‑140 (MQ=255)
                                gCCAGAAGCGCCGCGTCGACGCAAGCCGATGATTCTTCTGGATAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGTCTATTCAATGAAGATGACCCGCGATGCCGAAt                                                                                                          >  2:345145/1‑140 (MQ=255)
                                gCCAGAAGCGCCGCGTCGACGCAAGCCGATGATTCTTCTGGATAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAAt                                                                                                          <  1:59074/140‑1 (MQ=255)
                                     aaGCGCCGCGTCGACGCAAGCCGATGATTCTTCTGGATAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTTATTATGACGCGCTGTATGTCTATTCAATGAAGATGACCCGCGATGCCGAATACGAt                                                                                                     >  2:678873/1‑140 (MQ=255)
                                       gcgcCGCGTCGACGCAAGCCGATGATTCTTCTGGATAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGAttt                                                                                                   <  2:787873/140‑1 (MQ=255)
                                               tCGACGCAAGCCGATGATTCTTCTGGATAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGTCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATg                                                                                           >  2:681309/1‑140 (MQ=255)
                                                        gCCGATGATTCTTCTGGATAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGTCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAag                                                                                  >  2:956425/1‑140 (MQ=255)
                                                         ccGATGATTCTTCTGGATAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGTCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAagc                                                                                 <  1:837468/140‑1 (MQ=255)
                                                             tgatTCTTCTGGATAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGTCTATTCAATGAAGATGACCGCGATGCCGAATACGATTTAGTGCATGAGATGGAagccagcc                                                                            <  2:283865/140‑1 (MQ=255)
                                                                    tctGGATAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGTCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATgg                                                                      >  2:990697/1‑140 (MQ=255)
                                                                    tctGGATAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATgg                                                                      >  1:902675/1‑140 (MQ=255)
                                                                     ctGGATAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGTCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGa                                                                     >  1:421555/1‑140 (MQ=255)
                                                                           aaCATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGTCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGAt                                                               <  2:496986/140‑1 (MQ=255)
                                                                           aaCATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGTCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGAt                                                               <  1:956425/140‑1 (MQ=255)
                                                                                 cTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGTCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTc                                                         <  2:641468/140‑1 (MQ=255)
                                                                                  tGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGTCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTcc                                                        >  2:58827/1‑140 (MQ=255)
                                                                                  tGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTcc                                                        >  1:377693/1‑140 (MQ=255)
                                                                                  tGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTcc                                                        >  2:802248/1‑140 (MQ=255)
                                                                                    cGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGTCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAg                                                      >  1:451343/1‑140 (MQ=255)
                                                                                     gTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGt                                                     >  1:612623/1‑140 (MQ=255)
                                                                                      ttACTGCCTTCATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGTCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGCTGTCtcttatac                                                    <  2:744275/140‑8 (MQ=255)
                                                                                        aCTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGTCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGtctc                                                  >  1:227700/1‑140 (MQ=255)
                                                                                        aCTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGTCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGtctc                                                  >  1:214573/1‑140 (MQ=255)
                                                                                            ccTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAagc                                              >  1:386415/1‑140 (MQ=255)
                                                                                              ttGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGTCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAagcag                                            >  1:252959/1‑140 (MQ=255)
                                                                                              ttGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGTCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAagcag                                            <  2:421555/140‑1 (MQ=255)
                                                                                              ttGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGTCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAagcag                                            <  1:974833/140‑1 (MQ=255)
                                                                                              ttGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGTCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAagcag                                            >  2:94720/1‑140 (MQ=255)
                                                                                              ttGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAagcag                                            <  2:386415/140‑1 (MQ=255)
                                                                                                 atgatATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGTCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGt                                         >  1:331398/1‑140 (MQ=255)
                                                                                                        tttCAAAGGCTTCTTTGATTATGACGCGCTGAATGTCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAActg                                  >  1:495082/1‑140 (MQ=255)
                                                                                                           cAAAGGCTTCTTTGATTATGACGCGCTGAATGTCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAActgctg                               >  1:230162/1‑140 (MQ=255)
                                                                                                             aaGGCTTCTTTGATTATGACGCGCTGAATGTCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAg                             <  2:451343/140‑1 (MQ=255)
                                                                                                             aaGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAg                             <  1:983252/140‑1 (MQ=255)
                                                                                                               ggCTTCTTTGATTATGACGCGCTGAATGTCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGcc                           <  2:495082/140‑1 (MQ=255)
                                                                                                                gCTTCTTTGATTATGACGCGCTGAATGTCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCg                          <  1:94720/140‑1 (MQ=255)
                                                                                                                 cttcttTGATTATGACGCGCTGAATGTCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCgg                         >  2:195597/1‑140 (MQ=255)
                                                                                                                  ttcttTGATTATGACGCGCTGAATGTCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGt                        <  1:315967/140‑1 (MQ=255)
                                                                                                                   tcttTGATTATGACGCGCTGAATGTCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTg                       >  1:308253/1‑140 (MQ=255)
                                                                                                                    cttTGATTATGACGCGCTGAATGTCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGc                      >  1:279939/1‑140 (MQ=255)
                                                                                                                     ttttATTATGACGCGCTGTATGTCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCg                     <  1:678873/140‑1 (MQ=255)
                                                                                                                     tttGATTATTACGCGCTGAATGTCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCg                     <  2:599963/140‑1 (MQ=255)
                                                                                                                     tttGATTATGACGCGCTGAATGTCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCg                     <  1:630100/140‑1 (MQ=255)
                                                                                                                     tttGATTATGACGCGCTGAATGTCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCg                     <  1:675923/140‑1 (MQ=255)
                                                                                                                     tttGATTATGACGCGCTGAATGTCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCg                     <  2:261950/140‑1 (MQ=255)
                                                                                                                     tttGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCg                     <  2:902675/140‑1 (MQ=255)
                                                                                                                       tGATTATGACGCGCTGAATGTCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGtt                   <  2:966301/140‑1 (MQ=255)
                                                                                                                               aCGCGCTGAATGTCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTAt           >  1:376408/1‑140 (MQ=255)
                                                                                                                               aCGCGCTGAATGTCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTAt           <  1:990697/140‑1 (MQ=255)
                                                                                                                                  cgcTGAATGTCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAg        <  1:681309/140‑1 (MQ=255)
                                                                                                                                   gcTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAgc       >  1:944499/1‑140 (MQ=255)
                                                                                                                                    cTGAATGTCTATTCAATGAAGATGATCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAgc       >  1:53635/1‑139 (MQ=255)
                                                                                                                                    cTGAATGTCTATTCAATGAAGATGATCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAgc       <  2:53635/139‑1 (MQ=255)
                                                                                                                                       aaTGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAGCGCGa   <  2:612623/140‑1 (MQ=255)
                                                                                                                                        aTGTCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAGCGCGat  <  2:252959/140‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                        
TCAGATAAAGTGCCGCGCTTTGTGAATTTACCGCCAGAAGCGCCGCGTCGACGCAAGCCGATGATTCTTCTGGATAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAGCGCGAT  >  NC_000913/2623569‑2623844

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