Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 848,362 A→G 5.2% intergenic (‑358/+46) glnH ← / ← dps L‑glutamine ABC transporter periplasmic binding protein/stationary phase nucleoid protein that sequesters iron and protects DNA from damage

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913848,3620AG 5.2% 237.7 / 2.7 77intergenic (‑358/+46)glnH/dpsL‑glutamine ABC transporter periplasmic binding protein/stationary phase nucleoid protein that sequesters iron and protects DNA from damage
Reads supporting (aligned to +/- strand):  ref base A (37/36);  new base G (2/2);  total (39/38)
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.15e-01

TAAAACCTGATGAAAAAACAATCTTTATGTAACGATTGTGTGATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAGGATATCTGCGGTGTCGTCATCTTTCGCTTCGCCAATCG  >  NC_000913/848214‑848508
                                                                                                                                                    |                                                                                                                                                  
tAAAACCTGATGAAAAAACAATCTTTATGTAACGATTGTGTGATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAg                                                                                                                                                   <  1:1267388/150‑1 (MQ=255)
       tgatgaAAAAACAATCTTTATGTAACGATTGTGTGATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCttt                                                                                                                                            <  2:371893/150‑1 (MQ=255)
       tgatgaAAAAACAAGCTTTATGTAACGATTGGGTGATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCttt                                                                                                                                            <  2:487544/150‑1 (MQ=255)
         atgaAAAAACAATCTTTATGTAACGATTGTGTGATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATTCACTAAATAAGTTCGTTGAGGTGGTCTGCAAAACAAAACGGCGTCCTGTCAGGAAGCCGCTTTTa                                                                                                                                          <  2:256705/150‑1 (MQ=255)
          tgaAAAAACAATCTTTATGTAACGATTGTGTGATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTAt                                                                                                                                         >  2:457650/1‑150 (MQ=255)
          tgaAAAAACAATCTTTATGTAACGATTGTGTGATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTAt                                                                                                                                         >  2:114747/1‑150 (MQ=255)
            aaaaaaCAATCTTTATGTAACGATTGTGTGATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCg                                                                                                                                       >  1:730053/1‑150 (MQ=255)
             aaaaaCAATCTTTATGTAACGATTGTGTGATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCTTCCTGACAGGAGGCCGttttgttttg                                                                                                                                      <  2:611150/150‑7 (MQ=255)
              aaaaCAATCTTTATGTAACGATTGTGTGATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCTTCCTGACAGGAGGCCGttttgttttgc                                                                                                                                     <  1:1321839/150‑8 (MQ=255)
                aactttCTTTATTTAACGATTGTGTGATGATGTCGATACAAAAAATTTACTTTATCTGGTGTAGATAGCAGATGGATGCACTAAATAAGTGGGTGGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTa                                                                                                                                   <  2:635434/145‑1 (MQ=255)
                      cTTTATGTAACGATTGTGTGATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAg                                                                                                                             >  2:320222/1‑150 (MQ=255)
                           tGTAACGATTGTTTGATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTg                                                                                                                        <  1:441126/150‑1 (MQ=255)
                             tAACGATTGTGTGATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCa                                                                                                                      <  1:1210832/150‑1 (MQ=255)
                              aaCGATTGTGTGATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCAc                                                                                                                     <  1:842963/150‑1 (MQ=255)
                                    tgtgtgATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAg                                                                                                               <  1:705199/150‑1 (MQ=255)
                                      tgtgATGATGTGGATACATAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCg                                                                                                             >  1:634889/1‑150 (MQ=255)
                                      tgtgATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCg                                                                                                             >  1:1249009/1‑150 (MQ=255)
                                          atgatgTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATgg                                                                                                         <  1:630050/150‑1 (MQ=255)
                                             atgTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGAtt                                                                                                      <  1:114747/150‑1 (MQ=255)
                                              tgtgGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGAttt                                                                                                     >  1:895318/1‑150 (MQ=255)
                                                tgGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTAt                                                                                                   >  2:450272/1‑150 (MQ=255)
                                                    tACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGa                                                                                               <  2:608439/150‑1 (MQ=255)
                                                     aCAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATTGATTTATTCGAt                                                                                              <  2:615937/150‑1 (MQ=255)
                                                     aCAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGAt                                                                                              <  1:306888/150‑1 (MQ=255)
                                                      cAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCTGCGATGGATTTATTCGATg                                                                                             >  2:1163431/1‑150 (MQ=255)
                                                      cAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTc                                                                                                 <  1:1355277/146‑1 (MQ=255)
                                                      cAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTc                                                                                                 >  2:1355277/1‑146 (MQ=255)
                                                      cAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATg                                                                                             >  1:81243/1‑150 (MQ=255)
                                                      cAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATg                                                                                             >  1:672572/1‑150 (MQ=255)
                                                      cAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATg                                                                                             >  2:578900/1‑150 (MQ=255)
                                                      cAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATg                                                                                             <  1:946842/150‑1 (MQ=255)
                                                      cAAAAAATTTACTTAATCAGCTGGAGATAGAAGATGGATGCACTAAATAAGTGAGTTGAGGTGGGCTGCAAAACAAAACGGCCTCATGTCAGGAAGCCGCTTTTATCGGGTACTAAAGATCTGCACCATCAGCGATGGATATATTCGAtt                                                                                             >  1:230213/1‑149 (MQ=255)
                                                       aaaaaaTTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGt                                                                                            >  2:516720/1‑150 (MQ=255)
                                                         aaaaTTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTa                                                                                          <  2:1389850/150‑1 (MQ=255)
                                                            aTTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGAc                                                                                       <  1:11102/150‑1 (MQ=255)
                                                              ttACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTc                                                                                     <  2:730053/150‑1 (MQ=255)
                                                                 cTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGAt                                                                                  <  1:902769/150‑1 (MQ=255)
                                                                 attAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGAt                                                                                  <  2:230213/149‑1 (MQ=255)
                                                                   tAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATaa                                                                                >  1:1391903/1‑150 (MQ=255)
                                                                   tAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATaa                                                                                >  1:155071/1‑150 (MQ=255)
                                                                     aTCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAAc                                                                              >  1:73671/1‑150 (MQ=255)
                                                                            ggAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGa                                                                       >  1:774508/1‑150 (MQ=255)
                                                                                 tAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTa                                                                  >  2:768579/1‑150 (MQ=255)
                                                                                 tAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTa                                                                  >  1:821504/1‑150 (MQ=255)
                                                                                  aGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTAt                                                                 <  1:340498/150‑1 (MQ=255)
                                                                                    cAGATGGATGCACTAAATAAGTGCGTTGAGGTTGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATcc                                                               <  2:155071/150‑1 (MQ=255)
                                                                                      gatggatgCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAg                                                             >  1:118375/1‑150 (MQ=255)
                                                                                      cctcaacgCACTTATTTAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAg                                                             >  1:928595/8‑150 (MQ=255)
                                                                                        tggatgCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAATGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGt                                                           >  2:157327/1‑150 (MQ=255)
                                                                                              cACTAAATTAGTGCGTTGATGTGGGCTGCAAAACAAAAGGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTGGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCgag                                                     <  2:727524/150‑1 (MQ=255)
                                                                                              cACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCgag                                                     >  2:784471/1‑150 (MQ=255)
                                                                                               aCTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCgaga                                                    >  1:442827/1‑150 (MQ=255)
                                                                                               aCTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCgaga                                                    >  2:731582/1‑150 (MQ=255)
                                                                                                    aTAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACgcgg                                               <  1:1365839/150‑1 (MQ=255)
                                                                                                     tAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCgag                                                     <  2:1021535/143‑1 (MQ=255)
                                                                                                     tAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCgag                                                     >  1:1021535/1‑143 (MQ=255)
                                                                                                       aGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGcggcgg                                            >  2:923234/1‑150 (MQ=255)
                                                                                                       aGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGcggcgg                                            <  2:1249009/150‑1 (MQ=255)
                                                                                                        gTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGt                                           >  1:943139/1‑150 (MQ=255)
                                                                                                         tGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTc                                          <  2:626324/150‑1 (MQ=255)
                                                                                                          gCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCa                                         <  1:450272/150‑1 (MQ=255)
                                                                                                              tGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAGGat                                     <  1:335173/150‑1 (MQ=255)
                                                                                                                 ggTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAGGATATc                                  <  2:634889/150‑1 (MQ=255)
                                                                                                                    gggCTGCAAAACAAAACGGCCTCCTGTGAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGGCAGGATATCTGc                               <  1:757907/150‑1 (MQ=255)
                                                                                                                     ggCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAGGATATCTGCg                              >  1:1310158/1‑150 (MQ=255)
                                                                                                                       cTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAGGATATCTGCGgt                            >  2:860856/1‑150 (MQ=255)
                                                                                                                           aaaacaaaacGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCAGCGGTCAGGATATCTGCGGTgtcg                        >  2:1197749/1‑150 (MQ=255)
                                                                                                                            aaacaaaacGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAGGATATCTGCGGTgtcgt                       >  1:733949/1‑150 (MQ=255)
                                                                                                                            aaacaaaacGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAGGATATCTGCGGTgtcgt                       >  1:272658/1‑150 (MQ=255)
                                                                                                                                aaaacGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAGGATATCTGCGGTGTCGTCATc                   <  2:672572/150‑1 (MQ=255)
                                                                                                                                aaaacGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAGGATATCTGCGGTGTCGTCATc                   >  1:1035357/1‑150 (MQ=255)
                                                                                                                                 aagcGGCTTCCTGACAGGAGGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAGGATATCTGCGGTGTCGTCATCt                  >  2:109196/4‑150 (MQ=255)
                                                                                                                                   gcGGCTTCCTGACAGGAGGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAGGATATCTGCGGTGTCGTCATCttt                >  1:1431533/2‑150 (MQ=255)
                                                                                                                                      gCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAGGATATCTGCGGTGTCGTCATCTttcgc             <  1:11967/150‑1 (MQ=255)
                                                                                                                                            tGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAGGATATCTGCGGTGTCGTCATCTTTCGCTTCGcc       <  2:81243/150‑1 (MQ=255)
                                                                                                                                              tCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAGGATATCTGCGGTGTCGTCATCTTTCGCTTCGCCaa     <  1:516720/150‑1 (MQ=255)
                                                                                                                                                 ggAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAGGATATCTGCGGTGTCGTCATCTTTCGCTTCGCCAATCg  >  2:1112733/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TAAAACCTGATGAAAAAACAATCTTTATGTAACGATTGTGTGATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAGGATATCTGCGGTGTCGTCATCTTTCGCTTCGCCAATCG  >  NC_000913/848214‑848508

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