Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,367,342 C→T 6.2% E191K (GAA→AAA) ‡ pspF ← psp operon transcriptional activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,367,3420CT 6.2% 309.7 / 6.0 99E191K (GAA→AAA) ‡pspFpsp operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base C (44/47);  new base T (2/4);  total (46/53)
Fisher's exact test for biased strand distribution p-value = 6.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.71e-01

GGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCATACCAACCGCACA  >  NC_000913/1367195‑1367490
                                                                                                                                                   |                                                                                                                                                    
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:752328/1‑150 (MQ=255)
                                         atatTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGc                                                                                                           <  2:1162709/150‑1 (MQ=255)
                                         atatTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGc                                                                                                           <  2:1116159/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:1604187/150‑1 (MQ=255)
                                                              ttCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGcc                                                                                      <  2:637801/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:1569867/1‑150 (MQ=255)
                                                                                           aaaCCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGt                                                         >  2:1703466/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)
                                                                                                                                             aGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCATACCCAcc       >  1:541109/1‑150 (MQ=255)
                                                                                                                                              gTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTGCAACATCAAAAGCCAGTCGGTCGCGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCTTCGCCGGGAGATCGGCATTCGTCGCGCATACCAACCg      >  1:275539/1‑150 (MQ=255)
                                                                                                                                              gTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCATACCAACCg      <  2:510387/150‑1 (MQ=255)
                                                                                                                                              gTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCATACCAACCg      >  2:890529/1‑150 (MQ=255)
                                                                                                                                                  tCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACGACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCATACCACCCGcaca  >  2:548683/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCATACCAACCGCACA  >  NC_000913/1367195‑1367490

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: