Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,367,335 A→G 6.7% F193S (TTT→TCT)  pspF ← psp operon transcriptional activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,367,3350AG 6.7% 263.2 / 7.6 90F193S (TTT→TCT) pspFpsp operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base A (36/48);  new base G (2/4);  total (38/52)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.25e-01

GGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCATACCAAC  >  NC_000913/1367195‑1367484
                                                                                                                                            |                                                                                                                                                     
ggCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTg                                                                                                                                              <  2:1267054/150‑1 (MQ=255)
 gCGATACACTGAACGTTCCACCACGTTTTTCAATTCAAGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGc                                                                                                                                             <  2:1130450/150‑1 (MQ=255)
  cGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGcc                                                                                                                                            <  2:687619/150‑1 (MQ=255)
       cacTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAGAATACTTTGCCATcca                                                                                                                                       <  1:1059575/150‑3 (MQ=255)
        acTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAGAATACTTTGCCATccag                                                                                                                                      <  1:1412305/150‑4 (MQ=255)
        acTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAAc                                                                                                                                      <  1:1643453/150‑1 (MQ=255)
        acTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAAGGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAGAATACTTTGCCATccag                                                                                                                                      <  1:894417/150‑4 (MQ=255)
        acTGAACGTTCCACCACGTTTATCAATTCACGACTATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAGAATACTTTGCCATccag                                                                                                                                      <  1:263510/150‑4 (MQ=255)
         cTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACa                                                                                                                                     >  1:505579/1‑150 (MQ=255)
          tGAACGTTCCACCACGTTTATCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGAAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACAt                                                                                                                                    >  1:975220/1‑150 (MQ=255)
           gAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACAtt                                                                                                                                   <  1:1400048/150‑1 (MQ=255)
           gAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACAtt                                                                                                                                   >  1:1192440/1‑150 (MQ=255)
             aCGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTAt                                                                                                                                 >  2:1419677/1‑150 (MQ=255)
              cGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACACCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATg                                                                                                                                >  2:1711549/1‑150 (MQ=255)
                ttCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTc                                                                                                                              <  2:1499655/150‑1 (MQ=255)
                atacaccacGTTTTTCAATTCACGAATCTTTCCCGGCCAACGATAATTCAGCACTGTTTCTCTGGCGCTCTCCGTAAACCCCGGGAACAGAGGCATCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTc                                                                                                                              <  1:93688/147‑1 (MQ=255)
                  ccaccaCGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGc                                                                                                                            <  1:1239952/150‑1 (MQ=255)
                   caccacGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCt                                                                                                                           <  2:936778/150‑1 (MQ=255)
                   caccacGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCt                                                                                                                           <  1:1038316/150‑1 (MQ=255)
                     ccacGTTTTTCAATTCACGAATATATCCCGGCCTACTATAATTCACCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAAGAGAGGCAGATTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCttt                                                                                                                         <  1:143634/150‑1 (MQ=255)
                      caagTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTtc                                                                                                                        >  2:22160/4‑150 (MQ=255)
                       taatttttCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTtcg                                                                                                                       <  1:985841/147‑1 (MQ=255)
                         gTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCAGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTtcgcg                                                                                                                     <  2:256571/150‑1 (MQ=255)
                          tttttCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACTGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTtcgcgc                                                                                                                    <  1:391196/150‑1 (MQ=255)
                             ttCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTtcgcgctcg                                                                                                                 >  2:242996/1‑150 (MQ=255)
                                   tCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAgt                                                                                                           >  2:443028/1‑150 (MQ=255)
                                    cACGAATATTTCCCGGCCTACGATAATTCAGCAATGTTTCTCTGGCGCGCTGCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAgtg                                                                                                          <  2:340413/150‑1 (MQ=255)
                                       gAATATTTCACGGCCAAGGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCATCTTGATTACCCGACACATCTGGATGGCAAAGTATTATGCCATCAAGATTATGTCGCTTTCGCGCTCGCGCAgtggtg                                                                                                       <  2:186237/150‑1 (MQ=255)
                                         atatTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGc                                                                                                     <  2:970928/150‑1 (MQ=255)
                                         atatTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGc                                                                                                     <  2:1116159/150‑1 (MQ=255)
                                         atatTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGc                                                                                                     >  2:752328/1‑150 (MQ=255)
                                         atatTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGc                                                                                                     <  2:1162709/150‑1 (MQ=255)
                                            tttCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTTATATCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCTGTGGTGGCAga                                                                                                  >  1:80024/1‑149 (MQ=255)
                                                 cGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTac                                                                                             <  1:1609006/150‑1 (MQ=255)
                                                 cGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTac                                                                                             <  1:1424882/150‑1 (MQ=255)
                                                  ggCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTaca                                                                                            <  1:240754/150‑1 (MQ=255)
                                                   gCCAACGATAATTCAGGAATGTTTCTCTGGCGCGCCCGTAAACCACGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGCATTCTGCCATCAACATTATGTCGCTTTCGCTCTCGCGCAGTGGGGGCAGTTGTacaac                                                                                          <  1:470710/150‑1 (MQ=255)
                                                      aaCGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACTACAt                                                                                        >  1:137434/1‑150 (MQ=255)
                                                      aaCGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACAt                                                                                        >  1:1116458/1‑150 (MQ=255)
                                                       aCGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACACCATc                                                                                       >  2:847810/1‑150 (MQ=255)
                                                       aCGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATc                                                                                       >  2:780394/1‑150 (MQ=255)
                                                            aaTTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAg                                                                                  <  1:860157/150‑1 (MQ=255)
                                                              ttCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGcc                                                                                <  2:637801/150‑1 (MQ=255)
                                                              ttCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGcc                                                                                <  2:1604187/150‑1 (MQ=255)
                                                               tCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCa                                                                               >  1:1562510/1‑150 (MQ=255)
                                                               tCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTACCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTTGCAGTTGTCCAACATCAAACGCCa                                                                               >  2:911500/1‑150 (MQ=255)
                                                                   cAATGATACTCTGGCTCCCTACGTAAACACCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGCATGGCAAAGTATTCAGCCAACAACATTATGTCGCTATCGCGCTCGCGCAGAGGTGGCAGTGGTACAACATCAAAAGCCAgtcg                                                                           >  2:387944/1‑150 (MQ=255)
                                                                          tctcTGGCGCGCTCAGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCTAGc                                                                    <  2:593167/150‑1 (MQ=255)
                                                                            tctGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGCCACATCTGGATGGCAAAGTATTCTTCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGGGGAGGCAGTTGTACAACGTCAAAAGCCAGTCGGTCGAGccg                                                                  >  2:252733/1‑148 (MQ=255)
                                                                             ctGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAgg                                                                 <  2:909332/150‑1 (MQ=255)
                                                                              tGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGt                                                                >  1:172355/1‑150 (MQ=255)
                                                                               gtcgcgcTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGCTGGCAAAGTTTTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTc                                                               <  1:883440/148‑1 (MQ=255)
                                                                               ggCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCCTGAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTc                                                               >  1:1275711/1‑150 (MQ=255)
                                                                               ggCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTc                                                               <  1:400630/150‑1 (MQ=255)
                                                                                         gTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGAGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGCTCAGCGCGaaaa                                                     >  1:513212/1‑150 (MQ=255)
                                                                                         gTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTAGCTTTCGCGCTCTCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGaaaa                                                     >  2:1144805/1‑150 (MQ=255)
                                                                                          tAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAg                                                    >  1:330541/1‑150 (MQ=255)
                                                                                          tAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGGCAGCGCGAAAAg                                                    <  2:931012/150‑1 (MQ=255)
                                                                                           aaaCCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGt                                                   >  2:1703466/1‑150 (MQ=255)
                                                                                           aaaCCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGt                                                   >  2:1569867/1‑150 (MQ=255)
                                                                                            aaCCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTg                                                  <  1:391408/150‑1 (MQ=255)
                                                                                              ccccGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGcc                                                >  2:1074502/1‑150 (MQ=255)
                                                                                               cccGGGAACAGAGTCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGcgt                                               <  1:1178217/150‑3 (MQ=255)
                                                                                               cacGGGAACAGAGGGAGCTTGATTTCCCGACACAGCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCt                                               <  1:1178219/148‑1 (MQ=255)
                                                                                                 cGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTc                                             >  2:1532274/1‑150 (MQ=255)
                                                                                                 cGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTc                                             <  2:1021586/150‑1 (MQ=255)
                                                                                                     aaCAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACTTTATGTCGCTTTCGCGCTCGCGCAGTTGTGGCATTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCTGCGCGAAAAGTGCCTTCATTg                                         <  2:1356311/150‑1 (MQ=255)
                                                                                                     aaCAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTg                                         <  1:1260271/150‑1 (MQ=255)
                                                                                                      aCAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGa                                        <  2:1192440/150‑1 (MQ=255)
                                                                                                       cAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGAc                                       <  1:1515081/150‑1 (MQ=255)
                                                                                                           ggCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATc                                   >  2:1681782/1‑150 (MQ=255)
                                                                                                            gCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCg                                  <  1:22160/150‑1 (MQ=255)
                                                                                                                cTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCgg                              <  1:2922/150‑1 (MQ=255)
                                                                                                                cTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGGCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCgg                              <  1:1226862/150‑1 (MQ=255)
                                                                                                                 ttGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCggg                             >  1:1570045/1‑150 (MQ=255)
                                                                                                                      ttCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATc                        <  2:1625676/150‑1 (MQ=255)
                                                                                                                         ccGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGc                     <  1:1419677/150‑1 (MQ=255)
                                                                                                                          cGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCa                    <  1:242996/150‑1 (MQ=255)
                                                                                                                           gACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCTCCGGGAGATCGGCAt                   >  2:1242643/1‑150 (MQ=255)
                                                                                                                           gACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCAt                   >  2:1460112/1‑150 (MQ=255)
                                                                                                                             cacaTCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATtc                 >  1:358128/1‑150 (MQ=255)
                                                                                                                                  gttgATGGCAGAATACTTTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTcgc            >  2:902934/4‑150 (MQ=255)
                                                                                                                                  gttgATGGCAGAATACTTTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTcgc            >  2:987513/4‑150 (MQ=255)
                                                                                                                                   tGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGTGAGATCGGCATTCGTcgcg           >  1:163641/1‑150 (MQ=255)
                                                                                                                                    ggATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCCTTGACCATCGCCGGGAGCTCGGCATTCGTCgcgc          >  2:291570/1‑150 (MQ=255)
                                                                                                                                      gttgCAAAGGATGGTTCCATAAACATTATGTCGCTTTCTCGCTCGCTCAGTGGTGCCAGTTGTAAAACATCAAAAGCCAGTGGGTCGAGCAGTTCAGCGCGAAAAGTGCCTTAATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCAt        <  2:826541/147‑1 (MQ=255)
                                                                                                                                      aTGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGTCCATCGCCGGGAGATCGGCATTCGTCGCGCAt        <  2:116483/150‑1 (MQ=255)
                                                                                                                                         gCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCATAcc     >  2:43432/1‑150 (MQ=255)
                                                                                                                                         gCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCATAcc     <  1:1569867/150‑1 (MQ=255)
                                                                                                                                          cAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCATCCCa    >  1:1284838/1‑150 (MQ=255)
                                                                                                                                           aaaGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCATACCaa   >  1:838386/1‑150 (MQ=255)
                                                                                                                                           aaaGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCAGACCaa   >  1:285795/1‑150 (MQ=255)
                                                                                                                                           aaaGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATAGGCATTCGTCGCGCATACCaa   >  1:1610973/1‑150 (MQ=255)
                                                                                                                                            aaGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCATACCAAc  >  2:515073/1‑150 (MQ=255)
                                                                                                                                            |                                                                                                                                                     
GGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCATACCAAC  >  NC_000913/1367195‑1367484

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