Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,855,349 G→A 5.1% E338K (GAA→AAA) ‡ fhlA → DNA‑binding transcriptional activator FhlA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,855,3490GA 5.1% 316.0 / 3.4 97E338K (GAA→AAA) ‡fhlADNA‑binding transcriptional activator FhlA
Reads supporting (aligned to +/- strand):  ref base G (48/44);  new base A (3/2);  total (51/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.20e-01

TGTTGTTCGACACCTGGGGCAACCAGATTCAAACCTTGTGCCTGTTACCGCTGATGTCTGGCGACACCATGCTGGGCGTGCTGAAACTGGCGCAATGCGAAGAGAAAGTGTTTACCACTACCAATCTGAATTTACTGCGCCAGATTGCCGAACGTGTGGCAATCGCTGTCGATAACGCCCTCGCCTATCAGGAAATCCATCGTCTGAAAGAACGGCTGGTTGATGAAAACCTCGCCCTGACCGAGCAGCTCAACAATGTTGATAGTGAATTTGGCGAGATTATTGGCCGCAGCGAAG  >  NC_000913/2855200‑2855496
                                                                                                                                                     |                                                                                                                                                   
tgttgtTCGACACCTGGGGCAACCAGATTCAAACCTTGTGCCTGTTACCGCTGATGTCTGGCGACACCATGCTGGGCGTGCTGAAACTGGCGCAATGCGAAGAGAAAGTGTTTACCACTACCAATCTGAATTTACTGCGCCAGATTGCCg                                                                                                                                                     <  2:985915/150‑1 (MQ=255)
  ttgttCGACACCTGGGGCAACCAGATTCAAACCTTGTGCCTGTTACCGCTGATGTCTGGCGACACCATGCTGGGCGTGCTGAAACTGGCGCAATGCGAAGAGAAAGTGTTTACCACTACCAATCTGAATTTACTGCGCCAGATTGCCGaa                                                                                                                                                   >  1:1137694/1‑150 (MQ=255)
   tgttCGACACCTGGGGCAACCAGATTCAAACCTTGTGCCTGTTACCGCTGATGTCTGGCGACACCATGCTGGGCGTGCTGAAACTGGCGCAATGCGAAGAGAAAGTGTTTACCACTACCAATCTGAATTTACTGCGCCAGATTGCCGAAc                                                                                                                                                  <  2:1137694/150‑1 (MQ=255)
    gttCGACACCTGGGGCAACCAGATTCAAACCTTGTGCCTGTTACCGCTGATGTCTGGCGACACAATGCTGGGCGTGCTGAAACTGGCGCAATGCGAAGACAAAGTGTTTACCACTACCAATCTGAATTTACTGCGCCAGATTGCCGAACg                                                                                                                                                 >  2:1171119/1‑150 (MQ=255)
       cGACACCTGGGGCAACCAGATTCAAACCTTGTGCCTGTTACCGCTGATGTCTGGCGACACCATGCTGGGCGTGCTGAAACTGGCGCAATGCGAAGAGAAAGTGTTTACCACTACCAATCTGAATTTACTGCGCCAGATTGCCGAACgtgt                                                                                                                                              >  2:985888/1‑150 (MQ=255)
       cGACACCTGGGGCAACCAGATTCAAACCTTGTGCCTGTTACCGCTGATGTCTGGCGACACCATGCTGGGCGTGCTGAAACTGGCGCAATGCGAAGAGAAAGTGTTTACCACTACCAATCTGAATTTACTGCGCCAGATTGCCGAACgtgt                                                                                                                                              >  2:1036165/1‑150 (MQ=255)
       cGACACCTGGGGCAACCAGATTCAAACCTTGTGCCTGTTACCGCTGATGTCTGGCGACACCATGCTGGGCGTGCTGAAACTGGCGCAATGCGAAGAGAAAGTGTTTACCACTACCAATCTGAATTTACTGCGCCAGATTGCCGAACgtgt                                                                                                                                              >  1:493268/1‑150 (MQ=255)
        gACACCTGGGGCAACCAGATTCAAACCTTGTGCCTGTTACCGCTGATGTCTGGCGACACCATGCTGGGCGTGCTGAAACTGGCGCAATGCGAAGAGAAAGTGTTTACCACTACCAATCTGAATTTACTGCGCCAGATTGCCGAACGtgtg                                                                                                                                             >  1:185742/1‑150 (MQ=255)
          cagcTGGGGCAACCAGATTCAAACCTTGTGCCTGTTACCGCTGATGTCTGGCGACACCATGCTGGGCGTGCTGAAACTGGCGCAATGCGAAGAGAAAGTGTTTACCACTACCAATCTGAATTTACTGCGCCAGATTGCCGAACGTGTGGc                                                                                                                                           <  1:138812/147‑1 (MQ=255)
           acCTGGGGCAACCAGATTCAAACCTTGTGCCTGTTACCGCTGATGTCTGGCGACACCATGCTGGGCGTGCTGAAACTGGCGCAATGCGAAGAGAAAGTGTTTACCACTACCAATCTGAATTTACTGCGCCAGATTGCCGAACGTGTGGCa                                                                                                                                          <  1:1261175/150‑1 (MQ=255)
              tGGGGCAACCAGATTCAAACCTTGTGCCTGTTACCGCTGATGTCTGGCGACACCATGCTGGGCGTGCTGAAACTGGCGCAATGCGAAGAGAAAGTGTTTACCACTACCAATCTGAATTTACTGCGCCAGATTGCCGAACGTGTGGCAATc                                                                                                                                       >  1:1596166/1‑150 (MQ=255)
                gggCAACCAGATTCAAACCTTGTGCCTGTTACCGCTGATGTCTGGCGACACCATGCTGGGCGTGCTGAAACTGGCGCAATGCGAAGAGAAAGTGTTTACCACTACCAATCTGAATTTACTGCGCCAGATTGCCGAACGTGTGGCAATCGc                                                                                                                                     <  1:985888/150‑1 (MQ=255)
                 ggCAACCAGATTCAAACCTTGTGCCTGTTACCGCTGATGTCTGGCGACACCATGCTGGGCGTGCTGAAACTGGCGCAATGCGAAGAGAAAGTGTTTACCACTACCAATCTGAATTTACTGCGCCAGATTGCCGAACGTGTGGCAATCGCt                                                                                                                                    >  1:1190673/1‑150 (MQ=255)
                  gCAACCAGATTCAAACCTTGTGCCTGTTACCGCTGATGTCTGGCGACACCATGCTGGGCGTGCTGAAACTGGCGCAATGCGAAGAGAAAGTGTTTACCACTACCAATCTGAATTTACTGCGCCAGATTGCCACACGTTCGGCAATctggc                                                                                                                                   <  1:1295904/150‑5 (MQ=255)
                  gCAACCAGATTCAAACCTTGTGCCTGTTACCGCTGATGTCTGGCGACACCATGCTGGGCGTGCTGAAACTGGCGCAATGCGAAGAGAAAGTGTTTACCACTACCAATCTGAATTTACTGCGCCAGATTGCCACACGTTCGGCAATctggc                                                                                                                                   <  2:401341/150‑5 (MQ=255)
                    aaCCAGATTCAAACCTTGTGCCTGTTACCGCTGATGTCTGGCGACACCATGCTGGGCGTGCTGAAACTGGCGCAATGCGAAGAGAAAGTGTTTACCACTACCAATCTGAATTTACTGCGCCAGATTGCCGAACGTGTGGCAATCGCTGTc                                                                                                                                 >  1:631867/1‑150 (MQ=255)
                    aaCCAGATTCAAACCTTGTGCCTGTTACCGCTGATGTCTGGCGACACCATGCTGGGCGTGCTGAAACTGGCGCAATGCGAAGAGAAAGTGTTTACCACTACCAATCTGAATTTACTGCGCCAGATTGCCGAACGTGTGGCAATCGCTGTc                                                                                                                                 <  1:1110680/150‑1 (MQ=255)
                    aaCCAGATTCAAACCTTGTGCCTGTTACCGCTGATGTCTGGCGACACCATGCTGGGCGTGCTGAAACTGGCGCAATGCGAAGAGAAAGTGTTTACCACTACCAATCTGAATTTACTGCGCCAGATTGCCGAACGTGTGGCAATCGCTGTc                                                                                                                                 >  2:1065083/1‑150 (MQ=255)
                     aCCAGATTCAAACCTTGTGCCTGTTACCGCTGATGTCTGGCGACACCATGCTGGGCGTGCTGAAACTGGCGCAATGCGAAGAGAAAGTGTTTACCACTACCAATCTGAATTTACTGCGCCAGATTGCCGAACGTGTGGCAATCGCTGTCg                                                                                                                                >  2:822938/1‑150 (MQ=255)
                       cAGATTCAAACCTTGTGCCTGTTACCGCTGATGTCTGGCGACACCATGCTGGGCGTGCTGAAACTGGCGCAATGCGAAGAGAAAGTGTTTACCACTACCAATCTGAATTTACTGCGCCAGATTGCCGAACGTGTGGCAATCGCTGTCGAt                                                                                                                              >  2:232034/1‑150 (MQ=255)
                       cAGATTCAAACCTTGTGCCTGTTACCGCTGATGTCTGGCGACACCATGCTGGGCGTGCTGAAACTGGCGCAATGCGAAGAGAAAGTGTTTACCACTACCAATCTGAATTTACTGCGCCAGATTGCCGAACGTGTGGCAATCGCTGTCGAt                                                                                                                              <  1:1036165/150‑1 (MQ=255)
                           ttCAAACCTTGTGCCTGTTACCGCTGATGTCTGGCGACACCATGCTGGGCGTGCTGAAACTGGCGCAATGCGAAGAGAAAGTGTTTACCACTACCAATCTGAATTTACTGCGCCAGATTGCCGAACGTGTGGCAATCGCTGTCGATAACg                                                                                                                          >  1:565302/1‑150 (MQ=255)
                           ttCAAACCTTGTGCCTGTTACCGCTGATGTCTGGCGACACCATGCTGGGCGTGCTGAAACTGGCGCAATGCGAAGAGAAAGTGTTTACCACTACCAATCTGAATTTACTGCGCCAGATTGCCGAACGTGTGGCAATCGCTGTCGATAACg                                                                                                                          >  1:471893/1‑150 (MQ=255)
                              aaaCCTTGTGCCTGTTACCGCTGATGTCTGGCGACACCATGCTGGGCGTGCTGAAACTGGCGCAATGCGAAGAGAAAGTGTTTACCACTACCAATCTGAATTTACTGCGCCAGATTGCCGAACGTGTGGCAATCGCTGTCGATAACGccc                                                                                                                       >  2:450292/1‑150 (MQ=255)
                                 ccTTGTGCCTGTTACCGCTGATGTCTGGCGACACCATGCTGGGCGTGCTGAAACTGGCGCAATGCGAAGAGAAAGTGTTTACCACTACCAATCTGAATTTACTGCGCCAGATTGCCGAACGTGTGGCAATCGCTGTCGATAACGCCCTCg                                                                                                                    >  2:803243/1‑150 (MQ=255)
                                 ccTTGTGCCTGTTACCGCTGATGTCTGGCGACACCATGCTGGGCGTGCTGAAACTGGCGCAATGCGAAGAGAAAGTGTTTACCACTACCAATCTGAATTTACTGCGCCAGATTGCCGAACGTGTGGCAATCGCTGTCGATAACGCCCTCg                                                                                                                    >  1:74047/1‑150 (MQ=255)
                                   ttGTGCCTGTTACCGCTGATGTCTGGCGACACCATGCTGGGCGTGCTGAAACTGGCGCAATGCGAAGAGAAAGTGTTTACCACTACCAATCTGAATTTACTGCGCCAGATTGCCGAACGTGTGGCAATCGCTGTCGATAACGCCCTCGcc                                                                                                                  <  2:1534610/150‑1 (MQ=255)
                                       gCCTGTTACCGCTGATGTCTGGCGACACCATGCTGGGCGTGCTTAAACTGGCGCAATGCGAAGAGAAAGTGTTTACCACTACCAATCTGAATTTACTGCGCCAGATTGCCGAACGTGTGGCAATCGCTGTCGATAACGCCCTCGCCTATc                                                                                                              >  1:80370/1‑150 (MQ=255)
                                        ccTGTTACCGCTGATGTCTGGCGACACCATGCTGGGCGTGCTGAAACTGGCGCAATGCGAAGAGAAAGTGTTTACCACTACCAATCTGAATTTACTGCGCAAGATTGCCGAACGTGTGGCAATCGCTGTCGATAACGCCCTCGcc                                                                                                                  >  2:1107849/1‑145 (MQ=255)
                                        ccTGTTACCGCTGATGTCTGGCGACACCATGCTGGGCGTGCTGAAACTGGCGCAATGCGAAGAGAAAGTGTTTACCACTACCAATCTGAATTTACTGCGCAAGATTGCCGAACGTGTGGCAATCGCTGTCGATAACGCCCTCGcc                                                                                                                  <  1:1107849/145‑1 (MQ=255)
                                          tGTTACCGCTGATGTCTGGCGACACCATGCTGGGCGTGCTGAAACTGGCGCAATGCGAAGAGAAAGTGTTTACCACTACCAATCTGAATTTACTGCGCCAGATTGCCGAACGTGTGGCAATCGCTGTCGATACCGCCCGCGCCGATCAgg                                                                                                           >  1:778248/1‑150 (MQ=255)
                                          tGTTACCGCTGATGTCTGGCGACACCATGCTGGGCGTGCTGAAACTGGCGCAATGCGAAGAGAAAGTGTTTACCACTACCAATCTGAATTTACTGCGCCAGATTGCCGAACGTGTGGCAATCGCTGTCGATAACGCCCTCGCCTATCAgg                                                                                                           >  2:881121/1‑150 (MQ=255)
                                          tGTTACCGCTGATGTCTGGCGACACCATGCTGGGCGTGCTGAAACTGGCGCAATGCGAAGAGAAAGTGTTTACCACTACCAATCTGAATTTACTGCGCCAGATTGCCGAACGTGTGGCAATCGCTGTCGATAACGCCCTCGCCTATCAgg                                                                                                           >  1:784054/1‑150 (MQ=255)
                                                 gCTGATGTCTGGCGACACCATGCTGGGCGTGCTGAAACTGGCGCAATGCGAAGAGAAAGTGTTTACCACTACCAATCTGAATTTACTGCGCCAGATTGCCGAACGTGTGGCAATCGCTGTCGATAACGCCCTCGCCTATCAGGAAATCCa                                                                                                    >  1:130392/1‑150 (MQ=255)
                                                  cTGATGTCTGGCGACACCATGCTGGGCGTGCTGAAACTGGCGCAATGCGAAGAGAAAGTGTTTACCACTACCAATCTGAATTTACTGCGCCAGATTGCCGAACGTGTGGCAATCGCTGTCGATAACGCCCTCGCCTATCAGGAAATCCAt                                                                                                   >  2:1298858/1‑150 (MQ=255)
                                                     aTGTCTGGCGACACCATGCTGGGCGTGCTGAAACTGGCGCAATGCGAAGAGAAAGTGTTTACCACTACCACTCTGAATTTACTGCGCCAGATTGCCGAACGTGTGGCAATCGCTGTCGATAACGCCCTCGCCTATCAGGAAATCCATctt                                                                                                <  2:87631/150‑3 (MQ=255)
                                                      tGTCTGGCGACACCATGCTGGGCGTGCTGAAACTGGCGCAATGCGAAGAGAAAGTGTTTACCACTACCAATCTGAATTTACTGCGCCAGATTGCCGAACGTGTGGCAATCGCTGTCGATAACGCCCTCGCCTATCAGGAAATCCATCGTc                                                                                               <  1:822938/150‑1 (MQ=255)
                                                       gTCTGGCGACACCATGCTGGGCGTGCTGAAACTGGCGCAATGCGAAGAGAAAGTGTTTACCACTACCAATCTGAATTTACTGCGCCAGATTGCCGAACGTGTGGCAATCGCTGTCGATAACGCCCTCGCCTATCAGGAAATCCATCGTCt                                                                                              >  1:650698/1‑150 (MQ=255)
                                                        tCTGGCGACACCATGCTGGGCGTGCTGAAACTGGCGCAATGCGAAGAGAAAGTGTTTACCACTACCAATCTGAATTTACTGCGCCAGATTGCCGAACGTGTGGCAATCGCTGTCGATAACGCCCTCGCCTATCAGGAAATCCATCGTCTg                                                                                             <  2:185742/150‑1 (MQ=255)
                                                        tCTGGCGACACCATGCTGGGCGTGCTGAAACTGGCGCAATGCGAAGAGAAAGTGTTTACCACTACCAATCTGAATTTACTGCGCCAGATTGCCGAACGTGTGGCAATCGCTGTCGATAACGCCCTCGCCTATCAGGAAATCCATCGTCTg                                                                                             <  2:1540473/150‑1 (MQ=255)
                                                         cTGGCGACACCATGCTGGGCGTGCTGAAACTGGCGCAATGCGAAGAGAAAGTGTTTACCACTACCAATCTGAATTTACTGCGCCAGATTGCCGAACGTGTGGCAATCGCTGTCGATAACGCCCTCGCCTATCAGGAAATCCATCGTCTGa                                                                                            <  1:232034/150‑1 (MQ=255)
                                                         cTGGCGACACCATGCTGGGCGTGCTGAAACTGGCGCAATGCGAAGAGAAAGTGTTTACCACTACCAATCTGAATTTACTGCGCCAGATTGCCGAACGTGTGGCAATCGCTGTCGATAACGCCCTCGCCTATCAGGAAATCCATCGTCTGa                                                                                            >  1:528873/1‑150 (MQ=255)
                                                         cTGGCGACACCATGCTGGGCGTGCTGAAACTGGCGCAATGCGAAGAGAAAGTGTTTACCACTACCAATCTGAATTTACTGCGCCAGATTGCCGAACGTGTGGCAATCGCTGTCGATAACGCCCTCGCCTATCAGGAAATCCATCGTCTGa                                                                                            >  1:1345518/1‑150 (MQ=255)
                                                          tGGCGACACCATGCTGGGCGTGCTGAAACTGGCGCAATGCGAAGAGAAAGTGTTTACCACTACCAATCTGAATTTACTGCGCCAGATTGCCGAACGTGTGGCAATCGCTGTCGATAACGCCCTCGCCTATCAGGAAATCCATCGTCTGaa                                                                                           >  2:35931/1‑150 (MQ=255)
                                                           ggCGACACCATGCTGGGCGTGCTGAAACTGGCGCAATGCGAAGAGAAAGTGTTTACCACTACCAATCTGAATTTACTGCGCCAGATTGCCGAACGTGTGGCAATCGCTGTCGATAACGCCCTCGCCTATCAGGAAATCCATCGTCTGaaa                                                                                          >  2:1628356/1‑150 (MQ=255)
                                                                  ccATGCTGGGCGTGCTGAAACTGGCGCAATGCGAAGAGAAAGTGTTTACCACTACCAATCTGAATTTACTGCGCCAGATTGCCGAACGTGTGGCAATCGCTGTCGATAACGCCCTCGCCTATCAGGAAATCCATCGTCTGAAAGAACGGc                                                                                   >  2:1430181/1‑150 (MQ=255)
                                                                         gggCGTGCTGAAACTGGCGCAATGCGAAGAGAAAGTGTTTACCACTACCAATCTGAATTTACTGCGCCAGATTGCCGAACGTGTGGCAATCGCTGTCGATAACGCCCTCGCCTATCAGGAAATCCATCGTCTGAAAGAACGGCTGGTtga                                                                            >  1:924033/1‑150 (MQ=255)
                                                                          ggCGTGCTGAAACTGGCGCAATGCGAAGAGAAAGTGTTTACCACTACCAATCTGAATTTACTGCGCCAGATTGCCGAACGTGTGGCAATCGCTGTCGATAACGCCCTCGCCTATCAGGAAATCCATCGTCTGAAAGAACGGCTGGTtgat                                                                           <  2:631867/150‑1 (MQ=255)
                                                                          ggCGTGCTGAAACTGGCGCAATGCGAAGAGAAAGTGTTTACCACTACCAATCTGAATTTACTGCGCCAGATTGCCGAACGTGTGGCAATCGCTGTCGATAACGCCCTCGCCTATCAGGAAATCCATCGTCTGAAAGAACGGCTGGTtgat                                                                           <  1:450292/150‑1 (MQ=255)
                                                                           gCGTGCTGAAACTGGCGCAATGCGAAGAGAAAGTGTTTACCACTACCAATCTGAATTTACTGCGCCAGATTGCCGAACGTGTGGCAATCGCTGTCGATAACGCCCTCGCCTATCAGGAAATCCATCGTCTGAAAGAACGGCTGGTtgatg                                                                          <  2:1190673/150‑1 (MQ=255)
                                                                           gCGTGCTGAAACTGGCGCAATGCGAAGAGAAAGTGTTTACCACTACCAATCTGAATTTACTGCGCCAGATTGCCGAACGTGTGGCAATCGCTGTCGATAACGCCCTCGCCTATCAGGAAATCCATCGTCTGAAAGAACGGCTGGTtgatg                                                                          <  2:60182/150‑1 (MQ=255)
                                                                              tGCGGAAACTGGCGCAATGCGAAGAGAAAGTGTTTACCACTACCAATCTGAATTTACTGCGCCAGATTGCCGAACGTGTGGCAATCGCTGTCGATAACGCCCTCGCCTAGCAGGAAATCCATCGTCTGAAAGAACGGCTGGTtgatga                                                                         >  1:1661816/1‑148 (MQ=255)
                                                                              tGCCGAAACTGGCGCAATGCGAAGAGAAAGTGTTTACCACTACCAATCTGAATTTACTGCGCCAGATTGCCGAACGTGTGGCAATCGCTGTCGATAACGCCCTCGCCTATCAGGAAATCCATCGTCTGAAAGAACGGCTGGTtgatga                                                                         <  2:1661816/148‑1 (MQ=255)
                                                                                cTGAAACTGGCGCAATGCGAAGAGAAAGTGTTTACCACTACCAATCTGAATTTACTGCGCCAGATTGCCGAACGTGTGGCAATCGCTGTCGATAACGCCCTCGCCTATCAGGAAATCCATCGTCTGAAAGAACGGCTGGTTGATGAAAAc                                                                     <  1:1036577/150‑1 (MQ=255)
                                                                                cTGAAACTGGCGCAATGCGAAGAGAAAGTGTTTACCACTACCAATCTGAATTTACTGCGCCAGATTGCCGAACGTGTGGCAATCGCTGTCGATAACGCCCTCGCCTATCAGGAAAACCATCGTCTGAAAGAACGGCTGGTTGATGAAAAc                                                                     >  2:393217/1‑150 (MQ=255)
                                                                                       tGGCGCAATGCGAAGAGAAAGTGTTTACCACTACCAATCTGAATTTACTGCGCCAGATTGCCGAACGTGTGGCAATCGCTGTCGATAACGCCCTCGCCTATCAGGAAATCCATCGTCTGAAAGAACGGCTGGTTGATGAAAACCTCGccc                                                              <  2:109311/150‑1 (MQ=255)
                                                                                           gcAATGCGAAGAGAAAGTGTTTACCACTACCAATCTGAATTTACTGCGCCAGATTGCCGAACGTGTGGCAATCGCTGTCGATAACGCCCTCGCCTATCAGGAAATCCATCGTCTGAAAGAACGGCTGGTTGATGAAAACCTCGCCCTGAc                                                          >  2:289948/1‑150 (MQ=255)
                                                                                              aTGCGAAGAGAAAGTGTTTACCACTACCAATCTGAATTTACTGCGCCAGATTGCCGAACGTGTGGCAATCGCTGTCGATAACGCCCTCGCCTATCAGGAAATCCATCGTCTGAAAGAACGGCTGGTTGATGAAAACCTCGCCCTGACCGa                                                       <  1:1430181/150‑1 (MQ=255)
                                                                                               tGCGAAGAGAAAGTGTTTACCACTACCAATCTGAATTTACTGCGCCAGATTGCCGAACGTGTGGCAATCGCTGTCGATAACGCCCTCGCCTATCAGGAAATCCATCGTCTGAAAGAACGGCTGGTTGATGAAAACCTCGCCCTGACCGag                                                      <  2:345780/150‑1 (MQ=255)
                                                                                                gCGAAGAGAAAGTGTTTACCACTACCAATCTGAATTTACTGCGCCAGATTGCCGAACGTGTGGCAATCGCTGTCGATAACGCCCTCGCCTATCAGGAAATCCATCGTCTGAAAGAACGGCTGGTTGATGAAAACCTCGCCCTGACCGagc                                                     >  1:576349/1‑150 (MQ=255)
                                                                                                gCGAAGAGAAAGTGTTTACCACCACCAATCTGAATTTACTGCGCCAGATTGCCGAACGTGTGGCAATCGCTGTCAATAACGCCCTCGCCTATCAGGAAATCCATCGTCTGAAAGAACGGCTGGTTGATGAAAACCTCGCCCTGACCGagc                                                     <  1:661400/150‑1 (MQ=255)
                                                                                                 cGAAGAGAAAGTGTTTACCACTACCAATCTGAATTTACTGCGCCAGATTGCCGAACGTGTGGCAATCGCTGTCGATAACGCCCTCGCCTATCAGGAAATCCATCGTCTGAAAGAACGGCTGGTTGATGAAAACCTCGCCCTGACCGagca                                                    <  1:881121/150‑1 (MQ=255)
                                                                                                     gagaAAGTGTTTACCACTACCAATCTGAATTTACTGCGCCAGATTGCCGAACGTGTGGCAATCGCTGTCGATAACGCCCTCGCCTATCAGGAAATCCATCGTCTGAAAGAACGGCTGGTTGATGAAAACCTCGCCCTGACCGAGCAGCTc                                                <  2:1276267/150‑1 (MQ=255)
                                                                                                      agaAAGTGTTTACCACTACCAATCTGAATTTACTGCGCCAGATTGCCGAACGTGTGGCAATCGCTGTCGATAACGCCCTCGCCTATCAGGAAATCCATCGTCTg                                                                                             <  2:1236379/104‑1 (MQ=255)
                                                                                                      agaAAGTGTTTACCACTACCAATCTGAATTTACTGCGCCAGATTGCCGAACGTGTGGCAATCGCTGTCGATAACGCCCTCGCCTATCAGGAAATCCATCGTCTg                                                                                             >  1:1236379/1‑104 (MQ=255)
                                                                                                       gaAAGTGTTTACCACTACCAATCTGAATTTACTGCGCCAGATTGCCGAACGTGTGGCAATCGCTGTCGATAACGCCCTCGCCTATCAGGAAATCCATCGTCTGAAAGAACGGCTGGTTGATGAAAACCTCGCCCTGACCGAGCAGCTcaa                                              >  2:961473/1‑150 (MQ=255)
                                                                                                        aaaGTGTTTACCACTACCAATCTGAATTTACTGCGCCAGATTGCCGAACGTGTGGCAATCGCTGTCGATAACGCCCTCGCCTATCAGGAAATCCATCGTCTGAAAGAACGGCTGGTTGATGAAAACCTCGCCCTGACCGAGCAGCTcaac                                             >  2:1415031/1‑150 (MQ=255)
                                                                                                        aaaGTGTTTACCACTACCAATCTGAATTTACTGCGCCAGATTGCCGAACGTGTGGCAATCGCTGTCGATAACGCCCTCGCCTATCAGGAAATCCATCGTCTGAAAGAACGGCTGGTTGATGAAAACCTCGCCCTGACCGAGCAGCTcaac                                             >  1:1350664/1‑150 (MQ=255)
                                                                                                          aGTGTTTACCACTACCAATCTGAATTTACTGCGCCAGATTGCCGAACGTGTGGCAATCGCTGTCGATAACGCCCTCGCCTATCAGGAAATCCATCGTCTGAAAGAACGGCTGGTTGATGAAAACCTCGCCCTGACCGAGCAGCTcaacaa                                           >  1:1515413/1‑150 (MQ=255)
                                                                                                             gtTTACCACTACCAATCTGAATTTACTGCGCCAGATTGCCGAACGTGTGGCAATCGCTGTCGATAACGCCCTCGCCTATCAGGAAATCCATCGTCTGAAAGAACGGCTGGTTGATGAAAACCTCGCCCTGACCGAGCAGCTCAACAATGt                                        >  1:112519/1‑150 (MQ=255)
                                                                                                              tttACCACTACCAATCTGAATTTACTGCGCCAGATTGCCGAACGTGTGGCAATCGCTGTCGATAACGCCCTCGCCTATCAGGAAATCCATCGTCTGAAAGAACGGCTGGTTGATGAAAACCTCGCCCTGACCGAGCAGCTCAACAATGtt                                       <  2:471893/150‑1 (MQ=255)
                                                                                                               ttACCACTACCAATCTGAATTTACTGCGCCAGATTGCCGAACGTGTGGCAATCGCTGTCGATAACGCCCTCGCCTATCAGGAAATCCATCGTCTGAAAGAACGGCTGGTTGATGAAAACCTCGCCCTGACCGAGCAGCTCAACAATGTTg                                      <  2:130392/150‑1 (MQ=255)
                                                                                                                tACCACTACCAATCTGAATTTACTGCGCCAGATTGCCGAACGTGTGGCAATCGCTGTCGATAACGCCCTCGCCTATCAGGAAATCCATCGTCTGAAAGAACGGCTGGTTGATGAAAACCTCGCCCTGACCGAGCAGCTCAACAATGTTGa                                     >  2:1535012/1‑150 (MQ=255)
                                                                                                                 aCCACTAACAATCTGAATTTACTTCGCCAGTTTGCCGAACTTGTGGCAATCGCTGTCGATAACGCCCTCGCCTTTCAGGAAATCCATAGTCTGAAAGAACGGCTGGTTGATGAAAACCTCGCCCTGACCGAGCAGCTCAACAATGTTGAt                                    <  1:1171119/150‑1 (MQ=255)
                                                                                                                  ccACTACCAATCTGAATTTACTGCGCCAGATTGCCGAACGTGTGGCAATCGCTGTCGATAACGCCCTCGCCTATCAGGAAATCCATCGTCTGAAAGAACGGCTGGTTGATGAAAACCTCGCCCTGACCGAGCAGCTCAACAATGTTGATa                                   >  1:911132/1‑150 (MQ=255)
                                                                                                                   cACTACCAATCTGAATTTACTGCGCCAGATTGCCGAACGTGTGGCAATCGCTGTCGATAACGCCCTCGCCTATCAGGAAATCCATCGTCTGAAAGAACGGCTGGTTGATGAAAACCTCGCCCTGACCGAGCAGCTCAACAATGTTGATAg                                  <  1:1343873/150‑1 (MQ=255)
                                                                                                                       aCCAATCTGAATTTACTGCGCCAGATTGCCGAACGTGTGGCAATCGCTGTCGATAACGCCCTCGCCTATCAGGAAATCCATCGTCTGAAAGAACGGCTGGTTGATGAAAACCTCGCCCTGACCGAGCAGCTCAACAATGTTGATAGTGaa                              >  2:942386/1‑150 (MQ=255)
                                                                                                                            tCTGAATTTACTGCGCCAGATTGCCGAACGTGTGGCAATCGCTGTCGATAACGCCCTCGCCTATCAGGAAATCCATCGTCTGAAAGAACGGCTGGTTGATGAAAACCTCGCCCTGACCGAGCAGCTCAACAATGTTGATAGTGAATTTgg                         <  2:802437/150‑1 (MQ=255)
                                                                                                                             cTGAATTTACTGCGCCAGATTGCCGAACGTGTGGCAATCGCTGTCGATAACGCCCTCGCCTATCAGGAAATCCATCGTCTGAAAGAACGGCTGGTTGATGAAAACCTCGCCCTGACCGAGCAGCTCAACAATGTTGATAGTGAATTTGGc                        <  2:1345518/150‑1 (MQ=255)
                                                                                                                              tGAATTTACTGCGCCAGATTGCCGAACGTGTGGCAATCGCTGTCGATAACGCCCTCGCCTATCAGGAAATCCATCGTCTGAAAGAACGGCTGGTTGATGAAAACCTCGCCCTGACCGAGCAGCTCAACAATGTTGATAGTGAATTTGGCg                       >  1:672928/1‑150 (MQ=255)
                                                                                                                               gAATTTACTGCGCCAGATTGCCGAACGTGTGGCAATCGCTGTCGATAACGCCCTCGCCTATCAGGAAATCCATCGTCTGAAAGAACGGCTGGTTGATGAAAACCTCGCCCTGACCGAGCATCTCAACAATGTTGATAGTGAATTTGGCga                      <  2:924033/150‑1 (MQ=255)
                                                                                                                               gAATTTACTGCGCCAGATTGCCGAACGTGTGGCAATCGCTGTCGATAACGCCCTCGCCTATCAGGAAATCCATCGTCTGAAAGAACGGCTGGTTGATGAAAACCTCGCCCTGACCGAGCAGCTCAACAATGTTGATAGTGAATTTGGCga                      <  1:961473/150‑1 (MQ=255)
                                                                                                                                  tttACTGCGCCAGATTGCCGAACGTGTGGCAATCGCTGTCGATAACGCCCTCGCCTATCAGGAAATCCATCGTCTGAAAGAACGGCTGGTTGATGAAAACCTCGCCCTGACCGAGCAGCTCAACAATGTTGATAGTGAATTTGGCGAGat                   <  1:803243/150‑1 (MQ=255)
                                                                                                                                   ttACTGCGCCAGATTGCCGAACGTGTGGCAATCGCTGTCGATAACGCCCTCGCCTATCAGGAAATCCATCGTCTGAAAGAACGGCTGGTTGATGAAAACCTCGCCCTGACCGAGCAGCTCAACAATGTTGATAGTGAATTTGGCGAGatt                  <  2:493268/150‑1 (MQ=255)
                                                                                                                                    tACTGCGCCAGATTGCCGAACGTGTGGCAATCGCTGTCGATAACGCCCTCGCCTATCAGGAAATCCATCGTCTGAAAGAACGGCTGGTTGATGAAAACCTCGCCCTGACCGAGCAGCTCAACAATGTTGATAGTGAATTTGGCGAGatta                 >  1:650223/1‑150 (MQ=255)
                                                                                                                                      cTGCGCCAGATTGCCGAACGTGTGGCAATCGCTGTCGATAACGCCCTCGCCTATCAGGAAATCCATCGTCTGAAAGAACGGCTGGTTGATGAAAACCTCGCCCTGACCGAGCAGCTCAACAATGTTGATAGTGAATTTGGCGAGattatt               <  2:1596166/150‑1 (MQ=255)
                                                                                                                                       tGCGCCAGATTGCCGAACGTGTGGCAATCGCTGTCGATAACGCCCTCGCCTATCAGGAAATCCATCGTCTGAAAGAACGGCTGGTTGATGAAAACCTCGCCCTGACCGAGCAGCTCAACAATGTTGATAGTGAATTTGGCGAGATTATTg              <  2:576349/150‑1 (MQ=255)
                                                                                                                                       cgcgcCAGATTGCCGAACGTGTGGCAATCGCTGTCGCTAACGCCCTCGCCTATCCGGAAATCCATCGTCTGAAAGAACGGCTGGTTGATGAAAACCTCGCCCTGCCCGAGCACCTCAACAATGTTGATAGTGAATTTGGCGAGATTATTg              <  2:1570108/149‑1 (MQ=255)
                                                                                                                                        gcgcCAGATTGCCGAACGTGTGGCAATCGCTGTCGATAACGCCCTCGCCTATCAGGAAATCCATCGTCTGAAAGAACGGCTGGTTGATGAAAACCTCGCCCTGACCGAGCAGCTCAACAATGTTGATAGTGAATTTGGCGAGATTATTgg             >  1:653512/1‑150 (MQ=255)
                                                                                                                                         cgcCAGATTGCCGAACGTGTGGCAATCGCTGTCGATAACGCCCTCGCCTATCAGGAAATCCATCGTCTGAAAGAACGGCTGGTTGATGAAAACCTCGCCCTGACCGAGCAGCTCAACAATGTTGATAGTGAATTTGGCGAGATTATTGGc            <  1:791877/150‑1 (MQ=255)
                                                                                                                                         cgcCAGATTGCCGAACGTGGGGCAATCGCTGTCGATAACGCCCTCGCCTATCAGGAAATCCATCGTCTGAAAGAACGGCTGGTTGATGAAAACCTCGCCCTGACCGAGCAGCTCAACAATGTTGATAGTGAATTTGGCGAGATTATTGGc            <  2:796488/150‑1 (MQ=255)
                                                                                                                                         acagcgATTGCCACACGTTCGGCAATCGCTGTCGATAACGCCCTCGCCTATCAGGAAATCCATCGTCTGAAAGAACGGCTGGTTGATGAAAACCTCGCCCTGACCGAGCAGCTCAACAATGTTGATAGTGAATTTGGCGAGATTATTGGc            >  1:1566600/6‑150 (MQ=255)
                                                                                                                                          gcCAGATTGCCGAACGTGTGGCAATCGCTGTCGATAACGCCCTCGCCTATCAGGAAATCCATCGTCTGAAAGAACGGCTGGTTGATGAAAACCTCGCCCTGACCGAGCAGCTCAACAATGTTGATAGTGAATTTGGCGAGATTATTGGcc           <  1:1065083/150‑1 (MQ=255)
                                                                                                                                            gcgATTGCCACACGTTCGGCAATCGCTGTCGATAACGCCCTCGCCTATCAGGAAATCCATCGTCTGAAAGAACGGCTGGTTGATGAAAACCTCGCCCTGACCGAGCAGCTCAACAATGTTGATAGTGAATTTGGCGAGATTATTGGCCGc         >  1:1592452/3‑150 (MQ=255)
                                                                                                                                            gcgATTGCCACACGTTCGGCAATCGCTGTCGATAACGCCCTCGCCTATCAGGAAATCCATCGTCTGAAAGAACGGCTGGTTGATGAAAACCTCGCCCTGACCGAGCAGCTCAACAATGTTGATAGTGAATTTGGCGAGATTATTGGCCGc         >  1:681135/3‑150 (MQ=255)
                                                                                                                                               aTTGCCGAACGGGTGGCAATCGCTGTCGATAACGCCCTCGCCTATCAGGAAATCCATCGTCTGAAAGAACGGCTGGTTGATGAAAACCTCGCCCTGACCGAGCAGCTCAACAATGTTGATAGTGAATTTGGCGAGATTATTGGCCGCAGc      <  2:947151/150‑1 (MQ=255)
                                                                                                                                                  gCCGAACGTGTGGCAATCGCTGTCGATAACGCCCTCGCCTATCAGGAAATCCATCGTCTGAAAGAACGGCTGGTTGATGAAAACCTCGCCCTGACCGAGCAGCTCAACAATGTTGATAGTGAATTTGGCGAGATTATTGGCCGCAGCGaa   >  2:1359637/1‑150 (MQ=255)
                                                                                                                                                   ccGAACGTGTGGCAATCGCTGTCGATAACGCCCTCGCCTATCAGGAAATCCATCGTCTGAAAGAACGGCTGGTTGATGAAAACCTCACCCTGACCGAGCAGCTCAACAATGTTGATAGTGAATTTGGCGAGATTATTGGCCGCAGCGAAg  <  1:1298858/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
TGTTGTTCGACACCTGGGGCAACCAGATTCAAACCTTGTGCCTGTTACCGCTGATGTCTGGCGACACCATGCTGGGCGTGCTGAAACTGGCGCAATGCGAAGAGAAAGTGTTTACCACTACCAATCTGAATTTACTGCGCCAGATTGCCGAACGTGTGGCAATCGCTGTCGATAACGCCCTCGCCTATCAGGAAATCCATCGTCTGAAAGAACGGCTGGTTGATGAAAACCTCGCCCTGACCGAGCAGCTCAACAATGTTGATAGTGAATTTGGCGAGATTATTGGCCGCAGCGAAG  >  NC_000913/2855200‑2855496

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