Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,367,337 G→A 6.3% Y192Y (TAC→TAT pspF ← psp operon transcriptional activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,367,3370GA 6.3% 298.4 / 6.4 95Y192Y (TAC→TATpspFpsp operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base G (41/48);  new base A (2/4);  total (43/52)
Fisher's exact test for biased strand distribution p-value = 6.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.68e-01

GGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCATACCAACCG  >  NC_000913/1367195‑1367486
                                                                                                                                              |                                                                                                                                                     
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)
                                                                                                                                              |                                                                                                                                                     
GGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCATACCAACCG  >  NC_000913/1367195‑1367486

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: