Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_008463 1,128,160 N→A intergenic (+21/‑40) PA14_RS05280 → / → PA14_RS05285 DUF1302 domain‑containing protein/DUF1329 domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0084631,128,1600NA96.5% 1607.9 / NA 433intergenic (+21/‑40)PA14_RS05280/PA14_RS05285DUF1302 domain‑containing protein/DUF1329 domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base N (0/0);  major base A (193/225);  minor base C (6/9);  total (199/234)
Fisher's exact test for biased strand distribution p-value = 7.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 440 total aligned reads displayed.

TGGGGGTGGACGCGGAGTACCAGAACACCTACACCGCGAGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTNNNANCNGGNCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATGGCGGCGGTATCCGCCGACGAGGCGGCCAAGCTCGGCACCA  >  NC_008463/1128012‑1128308
                                                                                                                                                    |                                                                                                                                                    
tGGGGGTGGACGCGGAGTACCAGAACACCTACACCGCGAGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAAc                                                                                                                                                     >  1:9007621/1‑150 (MQ=255)
  ggggTGGACGCGGAGTACCAGAACACCTACACCGCGAGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCg                                                                                                                                                   <  2:2769165/150‑1 (MQ=255)
      tGGACGCGGAGTACCAGAACACCTACACCGCGAGGCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGTAAGTCTGATGTACGCCGGCCCCGTGCGAACCGGGcc                                                                                                                                               >  1:10252529/1‑150 (MQ=255)
             ggAGTACCAGAACACCTACACCGCGAGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCAGGAACTTCGGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGc                                                                                                                                        >  2:7812324/1‑150 (MQ=255)
              gAGTACCAGAACACCTACACCGCGAGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGGACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCa                                                                                                                                       >  2:10468895/1‑150 (MQ=255)
               aGTACCAGAACACCTACCCCGCGAGCCTCTCCTACACCAACTACTTCGCCGGCAAGTACACCACCGTGGACGATCGCGACTGCGTCGCGCTCAGCTTCGGCATGAACTACTCTTGTACGCCGGCCCGGGGCGAACCGGGCCGGCGAGCAg                                                                                                                                      <  1:4400250/150‑1 (MQ=255)
                gTACCAGAACACCTACACCGCGAGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAgg                                                                                                                                     <  2:1552820/150‑1 (MQ=255)
                gTACCAGAACACCTACACCGCGAGCCTCACCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAgg                                                                                                                                     <  2:7858044/150‑1 (MQ=255)
                   ccAGAACACCTACACCGCGAGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGACCCGGGCCGGCGAGCCGGACt                                                                                                                                  >  2:2662651/1‑150 (MQ=255)
                   ccAGAACACCTACACCGCGAGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACt                                                                                                                                  <  1:8058835/150‑1 (MQ=255)
                   ccAGAACACCTACACCGCGAGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACt                                                                                                                                  >  2:8060959/1‑150 (MQ=255)
                      gAACACCTACACCGCGAGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTaa                                                                                                                               <  1:3131420/150‑1 (MQ=255)
                       aaCACCTACACCGCGAGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAAc                                                                                                                              >  2:5042287/1‑150 (MQ=255)
                          acCTACACCGCGAGCCTCTCCTACACCACCTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAg                                                                                                                           <  2:6487006/150‑1 (MQ=255)
                          acCTACACCGCGAGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCCTGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAg                                                                                                                           >  2:9881737/1‑150 (MQ=255)
                            cTACACCGCGAGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGga                                                                                                                         <  1:6333355/150‑1 (MQ=255)
                               cacCGCGAGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGgacga                                                                                                                      >  2:5854023/1‑150 (MQ=255)
                                   gcgAGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACGTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAg                                                                                                                  <  2:7412575/150‑1 (MQ=255)
                                    cgAGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGc                                                                                                                 <  1:11092562/150‑1 (MQ=255)
                                    cgAGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGAATTCGTCGCGCTCAGCTTCGGCATGAACTTATGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGc                                                                                                                 <  1:2968958/150‑1 (MQ=255)
                                      aGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGa                                                                                                               <  2:9191461/150‑1 (MQ=255)
                                      aGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATAGCGACTTCGTCGCGCTCAGCATCGGCATGAACTTCTGATGTAAGCCCGCCCGGTGCGAAACGGGCCGGCGAGGCGGACGTAACGAGGACGgcgagcga                                                                                                               >  1:5611357/1‑150 (MQ=255)
                                       gCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGaa                                                                                                              >  2:10058183/1‑150 (MQ=255)
                                        ccTCTCCTACACCAACTTCTTCGACGGCAAGTACACCCCCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAAt                                                                                                             <  2:11105392/150‑1 (MQ=255)
                                             ccTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGaaaa                                                                                                        <  1:9649780/150‑1 (MQ=255)
                                             ccTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGaaaa                                                                                                        <  2:3982303/150‑1 (MQ=255)
                                               aacacCAACTTCTTCGACGGCAAGTACACCACCGTGGACCATCGCTACTTCGACGCGCTCAGCTTCGGCATGAACTTCTGCTGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAaaca                                                                                                      <  1:4644514/149‑1 (MQ=255)
                                                acacCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAaacaa                                                                                                     <  1:8560911/150‑1 (MQ=255)
                                                 cacCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGACCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAaacaac                                                                                                    <  2:5313139/150‑1 (MQ=255)
                                                 cacCAACTTCTTCGACGGAAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACTAGGACGACGAGCGAATGAAaacaac                                                                                                    <  2:6244046/150‑1 (MQ=255)
                                                      aCTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGa                                                                                               <  1:9901356/150‑1 (MQ=255)
                                                      aCTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGa                                                                                               <  1:10419703/150‑1 (MQ=255)
                                                        ttcttcGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGAtt                                                                                             <  2:6842891/150‑1 (MQ=255)
                                                         tcttcGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGAttc                                                                                            <  1:3425556/150‑1 (MQ=255)
                                                         tcttcGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGAttc                                                                                            >  1:2547434/1‑150 (MQ=255)
                                                          cttcGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGAttct                                                                                           <  2:8695491/150‑1 (MQ=255)
                                                            tcGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTAGCGCTCAGCTTAGGCAAGAACTTAGGATGTACGCCGGCCCCGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGACTGAAAACAACGAAGAttcttc                                                                                         >  1:2103206/1‑150 (MQ=255)
                                                               aCGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGc                                                                                      <  1:3493272/150‑1 (MQ=255)
                                                                cGGCAAGTACACCACCGTGGACGCTCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCa                                                                                     <  1:8775855/150‑1 (MQ=255)
                                                                cGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCCGGGCGAACCGGGCCGGCGAGCAGGACATAACGAGGACGACGAGCGACTGAAAACAACGAAGTTTATTATGCa                                                                                     >  1:1075927/1‑150 (MQ=255)
                                                                cGGCAAGTACACCACCGTGGACGAGCGCGACTTCGTCGCGCTCAGCTTCGGCCTGAACTTCTGATGTACGCCGGCCCGGTGCGAACCTGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGCAAGAAAACAACGAAGAGTCTTCTGCa                                                                                     >  2:3081711/1‑150 (MQ=255)
                                                                  gCAAGTACCCCCCCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACGTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGACTGAAAACAACGAAGATTCTTCTGCATa                                                                                   <  2:2881503/150‑1 (MQ=255)
                                                                  gCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCCTGAACTTCTGATGTACGCCGGCCCGGTGCGCACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGCATGAAAACAACGAAGATTCTTCTGCATa                                                                                   >  2:4876643/1‑150 (MQ=255)
                                                                  gCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATa                                                                                   <  2:9089927/150‑1 (MQ=255)
                                                                  gCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATa                                                                                   <  2:3468439/150‑1 (MQ=255)
                                                                  gCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATa                                                                                   >  2:67653/1‑150 (MQ=255)
                                                                    aaGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATAcc                                                                                 <  1:10011253/150‑1 (MQ=255)
                                                                     aGTACACCACCGTGGACGATCGCGTCTTCGTCGCGCTCAGCTTCGGCCTGAACTTCTGATGTACGCCGGCCCGGTGCGCACCGGGCCGGCGAGCAGGACGTAACGAGGACGACGCGCGAATGAAAACAACGCAGATTCTTCGGCATACCg                                                                                >  2:7135928/1‑150 (MQ=255)
                                                                     aGTACACCACCGTGGACGATCGCGACTTCGGCGCGCTCAGCTTCGGCATGAACTTCTGATGTAAGCCGGCCCGGTGCGTACCGGGCCGGCGAGCAGGACTTAACGAGGAAGACGAGCGAATGAAAACAACGTAGATTATTCTGCATACCg                                                                                >  1:1538605/1‑150 (MQ=255)
                                                                      gTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCgg                                                                               <  1:4577303/150‑1 (MQ=255)
                                                                       tACACCACCGTGGACGATCGCGACTGCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGACCCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGgt                                                                              <  1:3685364/150‑1 (MQ=255)
                                                                           ccaccGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTAc                                                                          <  1:5295488/150‑1 (MQ=255)
                                                                             accGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTg                                                                        <  1:7522273/150‑1 (MQ=255)
                                                                             accGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTg                                                                        <  2:2934372/150‑1 (MQ=255)
                                                                                 tGGACGATCGCGACATCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTAAGCCGGCCCGGTGCGAACCGGGCCCGCGAGCAGGACTTAAAGAGGACGACGAGCGAAAGAAAACAACGAAGATTCTTATGCATACCGGTGTAATGGccc                                                                    >  2:6977761/1‑150 (MQ=255)
                                                                                  ggACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTAATGGCCCt                                                                   >  1:10638591/1‑150 (MQ=255)
                                                                                   gACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTg                                                                  <  1:4011128/150‑1 (MQ=255)
                                                                                      gATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCGTCTGCATACCGGTGTACTGGCCCTGAGc                                                               <  2:4499297/150‑1 (MQ=255)
                                                                                       aTCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCCGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGcc                                                              >  2:4163210/1‑150 (MQ=255)
                                                                                       aTCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGcc                                                              <  1:10129403/150‑1 (MQ=255)
                                                                                          gcgACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCctgc                                                           <  1:6661842/150‑1 (MQ=255)
                                                                                           cgACTTCGTCGCGCTCAGCTGCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCctgct                                                          >  2:470083/1‑150 (MQ=255)
                                                                                            gACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCctgctg                                                         <  2:8617275/150‑1 (MQ=255)
                                                                                              cTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGACCCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGc                                                       <  1:1878660/150‑1 (MQ=255)
                                                                                                   tcgCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCAccc                                                  >  2:10991723/1‑150 (MQ=255)
                                                                                                    cgcgCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCa                                                 >  1:9311640/1‑150 (MQ=255)
                                                                                                      cgcTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAgg                                               <  1:5096883/150‑1 (MQ=255)
                                                                                                       gcTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGt                                              <  1:1204153/150‑1 (MQ=255)
                                                                                                       gcTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGt                                              >  1:7006250/1‑150 (MQ=255)
                                                                                                        cTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTg                                             >  2:772824/1‑150 (MQ=255)
                                                                                                        cTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTg                                             >  1:7748055/1‑150 (MQ=255)
                                                                                                        cTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTg                                             >  1:7673115/1‑150 (MQ=255)
                                                                                                          cAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGAt                                           <  2:2247932/150‑1 (MQ=255)
                                                                                                          cAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGAt                                           >  2:6681038/1‑150 (MQ=255)
                                                                                                          cAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGAt                                           <  1:7334745/150‑1 (MQ=255)
                                                                                                           aGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATg                                          <  2:6124047/150‑1 (MQ=255)
                                                                                                            gCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATgg                                         >  1:2789663/1‑150 (MQ=255)
                                                                                                             cTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATggc                                        <  1:1380533/150‑1 (MQ=255)
                                                                                                             cTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGGGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATggc                                        <  1:4907743/150‑1 (MQ=255)
                                                                                                                  gCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATGGcggcgg                                   <  2:2130641/150‑1 (MQ=255)
                                                                                                                    aTGACCTTCTGATCTACGCCGGCCCGGTGCGCCCCGGGCCGGCGAGCAGGACTTAACGCGGCCGACGAGCGAATGAAAACAACGAAGATTCTTCGGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATGGCGGCGGTa                                 <  2:8836158/150‑1 (MQ=255)
                                                                                                                    aTGAACGTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAATATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATGGCGGCGGTa                                 <  2:4290017/150‑1 (MQ=255)
                                                                                                                      gAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACATAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATGGCGGCGGTATc                               <  2:10746185/150‑1 (MQ=255)
                                                                                                                        aCTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATGGCGGCGGTATccg                             <  1:6272535/150‑1 (MQ=255)
                                                                                                                            cTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATGGCGGCGGTATCCGCcga                         <  2:1951879/150‑1 (MQ=255)
                                                                                                                              gATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATGGCGGCGGTATCCGCcgacg                       <  1:8380119/150‑1 (MQ=255)
                                                                                                                                  tACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTAAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATGGCGGCGGTATCCGCCGACGAggc                   <  2:2290973/150‑1 (MQ=255)
                                                                                                                                   aCGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTGAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCCCCCAGGTGATGGCGGCGGTATCCGCCGACGAggcg                  >  1:2445999/1‑150 (MQ=255)
                                                                                                                                     gCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATGGCGGCGGTATCCGCCGACGAggcggc                >  2:3795631/1‑150 (MQ=255)
                                                                                                                                       cGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGATGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATGGCGGCGGTATCCGCCGACGAGGCGGCCa              >  1:10800112/1‑150 (MQ=255)
                                                                                                                                        ggCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATGGCGGCGGTATCCGCCGACGAGGCGGCCaa             >  2:8163667/1‑150 (MQ=255)
                                                                                                                                         gCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATGGCGGCGGTATCCGCCGACGAGGCGGCCAAg            >  2:9554130/1‑150 (MQ=255)
                                                                                                                                           ccGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATGGCGGCGGTATCCGCCGACGAGGCGGCCAAGCt          <  2:8398167/150‑1 (MQ=255)
                                                                                                                                           ccGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATGGCGGCGGTATCCGCCGACGAGGCGGCCAAGCt          >  1:6467364/1‑150 (MQ=255)
                                                                                                                                             ggTGCGAACCGGGCCGGCGAGCAGGACTTGACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATGGCGGCGGTATCCGCCGACGAGGCGGCCAAGCTCg        >  2:5595001/1‑150 (MQ=255)
                                                                                                                                             ggTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATGGCGGCGGTATCCGCCGACGAGGCGGCCAAGCTCg        >  1:9278963/1‑150 (MQ=255)
                                                                                                                                             ggTGCGAACCGGGCCGGCGAGCAGGACTAAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATGGCGGCGGTATCCGCCGACGAGGCGGCCAAGCTCg        >  1:5003703/1‑150 (MQ=255)
                                                                                                                                                 cgaCCCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAAAAACGAAGATTCTACTGCATACCGTTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATGGCGGCGGTATCCGCCTAAGAGGCGGCCAAGCTCGGAAc    >  2:825074/3‑150 (MQ=255)
                                                                                                                                                   aaCCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATGGCGGCGGTATCCGCCGACGAGGCGGCCAAGCTCGGCACCa  <  2:4781361/150‑1 (MQ=255)
                                                                                                                                                   aCCCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATGGCGGCGGTATCCGCCGACGAGGCGGCCAAGCTCGGCACCa  >  2:9362345/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TGGGGGTGGACGCGGAGTACCAGAACACCTACACCGCGAGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTNNNANCNGGNCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATGGCGGCGGTATCCGCCGACGAGGCGGCCAAGCTCGGCACCA  >  NC_008463/1128012‑1128308

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