Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,687,214 T→G 8.9% V65V (GTA→GTC)  glnB ← nitrogen regulatory protein PII‑1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,687,2140TG 8.9% 195.3 / 9.4 66V65V (GTA→GTC) glnBnitrogen regulatory protein PII‑1
Reads supporting (aligned to +/- strand):  ref base T (26/34);  new base G (4/2);  total (30/36)
Fisher's exact test for biased strand distribution p-value = 3.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGGTTTAAATTGCCGCGTCGTCCTCTTCACCGGTACGGATGCGAATGACCCGTGCCACGTCAAAGACGAAGATTTTACCGTCACCGATTTTGCCGGTTTGCGCCGTGCGAATAATGGTATCGACACAGGTATCGACAATGTCGTCCGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATCCACCATATACTCCGCGCCGCGGTACAGCTCGGTATGGCCTTTCTGGCGACCAAAGCCTTTCACTTCGGTCACCGTCATGCCGGTAATACCGACTTCGGCCAGTGCTTCGCGGAC  >  NC_000913/2687065‑2687363
                                                                                                                                                     |                                                                                                                                                     
gcgGTTTAAATTGCCGCGTCGTCCTCTTCACCGGTACGGATGCGAATGACCCGTGCCACGTCAAAGACGAAGATTTTACCGTCACCGATTTTGCCGGTTTGCGCCGTGCGAATAATGGTATCGACACAGGTATCGACAATGTCGTCCGGt                                                                                                                                                       >  2:1262965/1‑150 (MQ=255)
   gTTTAAATTGCCGCGTCGTCCTCTTCACCGGTACGGATGCGAATGACCCGTGCCACGTCAAAGACGAAGATTTTACCGTCACCGATTTTGCCGGTTTGCGCCGTGCGAATAATGGTATCGACACAGGTATCTACAATGTCGTCCGGTACg                                                                                                                                                    <  2:951934/150‑1 (MQ=255)
      tAAATTGCCGCGTCGTCCTCTTCACCGGTACGGATGCGAATGACCCGTGCCACGTCAAAGACGAAGATTTTACCGTCACCGATTTTGCCGGTTTGCGCCGTGCGAATAATGGTATCGACACAGGTATCGACAATGTCGTCCGGTACGACa                                                                                                                                                 >  2:1408790/1‑150 (MQ=255)
       aaaTTGCCGCGTCGTCCTCTTCACCGGTACGGATGCGAATGACCCGTGCCACGTCAAAGACGAAGATTTTACCGTCACCGATTTTGCCGGTTTGCGCCGTGCGAATAATGGTATCGACACAGGTATCGACAATGTCGTCCGGTACGACaa                                                                                                                                                >  1:22925/1‑150 (MQ=255)
       aaaTTGCCGCGTCGTCCTCTTCACCGGTACGGATGCGAATGACCCGTGCCACGTCAAAGACGAAGATTTTACCGTCACCGATTTTGCCGGTTTGCGCCGTGCGAATAATGGTATCGACACAGGTATCGACAATGTCGTCCGGTACGACaa                                                                                                                                                >  1:189790/1‑150 (MQ=255)
       aaaTTGCCGCGTCGTCCTCTTCACCGGTACGGATGCGAATGACCCGTGCCACGTCAAAGACGAAGATTTTACCGTCACCGATTTTGCCGGTTTGCGCCGTGCGAATAATGGTATCGACACAGGTATCGACAATGTCGTCCGGTACGACaa                                                                                                                                                >  1:17524/1‑150 (MQ=255)
         aTTGCCGCGTCGTCCTCTTCACCGGTACGGATGCGAATGACCCGTGCCACGTCAAAGACGAAGATTTTACCGTCACCGATTTTGCCGGTTTGCGCCGTGCGAATAATGGTATCGACACAGGTATCGACAATGTCGTCCGGTACGACAAtc                                                                                                                                              <  1:1383215/150‑1 (MQ=255)
          ttGCCGCGTCGTCCTCTTCACCGGTACGGATGCGAATGACCCGTGCCACGTCAAAGACGAAGATTTTACCGTCACCGATTTTGCCGGTTTGCGCCGTGCGAATAATGGTATCGACACAGGTATCGACAATGTCGTACCGGACGACattgt                                                                                                                                             <  2:913814/150‑5 (MQ=255)
            gCCGCGTCGTCCTCTTCACCGGTACGGATGCGAATGACCCGTGCCACGTCAAAGACGAAGATTTTACCGTCACCGATTTTGCCGGTTTGCGCCGTGCGAATAATGGTATCGACACAGGTATCGACAATGTCGTCCGGTACGACAATCTCa                                                                                                                                           >  2:723098/1‑150 (MQ=255)
                  tcgtcCTCTTCACCGGTACGGATGCGAATGACCCGTGCCACGTCAAAGACGAAGATTTTACCGTCACCGATTTTGCCGGTTTGCGCCGTGCGAATAATGGTATCGACACAGGTATCGACAATGTCGTCCGGTACGACAATCTCAATTTTc                                                                                                                                     <  2:1394046/150‑1 (MQ=255)
                     tcCTCTTCACCGGTACGGATGCGAATGACCCGTGCCACGTCAAAGACGAAGATTTTACCGTCACCGATTTTGCCGGTTTGCGCCGTGCGAATAATGGTATCGACACAGGTATCGACAATGTCGTCCGGTACGACAATCTCAATTTTCACt                                                                                                                                  <  1:784439/150‑1 (MQ=255)
                     tcCTCTTCACCGGTACGGATGCGAATGACCCGTGCCACGTCAAAGACGAAGATTTTACCGTCACCGATTTTGCCGGTTTGCGCCGTGCGAATAATGGTATCGACACAGGTATCGACAATGTCGTCCGGTACGACAATCTCAATTTTCACt                                                                                                                                  <  2:504794/150‑1 (MQ=255)
                        tctTCACCGGTACGGATGCGAATGACCCGTGCCACGTCAAAGACGAAGATTTTACCGTCACCGATTTTGCCGGTTTGCGCCGTGCGAATAATGGTATCGACACAGGTATCGACAATGTCGTCCGGTACGACAATCTCAATTTTCACTTTc                                                                                                                               <  1:315432/150‑1 (MQ=255)
                          ttCACCGGTACGGATGCGAATGACCCGTGCCACGTCAAAGACGAAGATTTTACCGTCACCGATTTTGCCGGTTTGCGCCGTGCGAATAATGGTATCGACACAGGTATCGACAATGTCGTCCGGTACGACAATCTCAATTTTCACTTTCgg                                                                                                                             >  1:351643/1‑150 (MQ=255)
                           tCACCGGTACGGATGCGAATGACCCGTGCCACGTCAAAGACGAAGATTTTACCGTCACCGATTTTGCCGGTTTGCGCCGTGCGAATAATGGTATCGACACAGGTATCGACAATGTCGTCCGGTACGACAATCTCAATTTTCACTTTCGGc                                                                                                                            <  2:351643/150‑1 (MQ=255)
                            cACCGGTACGGATGCGAATGACCCGTGCCACGTCAAAGACGAAGATTTTACCGTCACCGATTTTGCCGGTTTGCGCCGTGCGAATAATGGTATCGACACAGGTATCGACAATGTCGTACCGGACGACattgtcg                                                                                                                                           >  2:1267888/1‑128 (MQ=255)
                            cACCGGTACGGATGCGAATGACCCGTGCCACGTCAAAGACGAAGATTTTACCGTCACCGATTTTGCCGGTTTGCGCCGTGCGAATAATGGTATCGACACAGGTATCGACAATGTCGTACCGGACGACattgtcg                                                                                                                                           <  1:1267888/134‑7 (MQ=255)
                                     ggATGCGAATGACCCGTGCCACGTCAAAGACGAAGATTTTACCGTCACCGATTTTGCCGGTTTGCGCCGTGCGAATAATGGTATCGACACAGGTATCGACAATGTCGTCCGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATcca                                                                                                                  <  1:388950/150‑1 (MQ=255)
                                      gATGCGAATGACCCGTGCCACGTCAAAGACGAAGATTTTACCGTCACCGATTTTGCCGGTTTGCGCCGTGCGAATAATGGTATCGACACAGGTATCGACAATGTCGTCCGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATccac                                                                                                                 <  1:1377895/150‑1 (MQ=255)
                                       aTGCGAATGACCCGTGCCACGTCAAAGACGAAGATTTTACCGTCACCGATTTTGCCGGTTTGCGCCGTGCGAATAATGGTATCGACACAGGTATCGACAATGTCGTCCGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATccacc                                                                                                                <  1:845038/150‑1 (MQ=255)
                                       aTGCGAATGACCCGTGCCACGTCAAAGACGAAGATTTTACCGTCACCGATTTTGCCGGTTTGCGCCGTGCGAATAATGGTATCGACACAGGTATCGACAATGTCGTCCGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATccacc                                                                                                                <  1:1408790/150‑1 (MQ=255)
                                             aTGACCCGTGCCACGTCAAAGACGAAGATTTTACCGTCACCGATTTTGCCGGTTTGCGCCGTGCGAATAATGGTATCGACACAGGTATCGACAATGTCGTCCGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATCCACCATATAc                                                                                                          <  1:936013/150‑1 (MQ=255)
                                               gACCCGTGCCACGTCAAAGACGAAGATTTTACCGTCACCGATTTTGCCGGTTTGCGCCGTGCGAATAATGGTATCGACACAGGTATCGACAATGTCGTCCGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATCCACCATATACTc                                                                                                        >  2:129249/1‑150 (MQ=255)
                                                aCCCGTGCCACGTCAAAGACGAAGATTTTACCGTCACCGATTTTGCCGGTTTGCGCCGTGCGAATAATGGTATCGACACAGGTATCGACAATGTCGTCCGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATCCACCATATACTcc                                                                                                       >  2:1172856/1‑150 (MQ=255)
                                                          cGTCAAAGACGAAGATTTTACCGTCACCGATTTTGCCGGTTTGCGCCGTGCGAATAATGGTATCGACACAGGTATCGACAATGTCGTCCGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATCCACCATATACTCCGCGCCGCGGt                                                                                             <  1:723098/150‑1 (MQ=255)
                                                           gTCAAAGACGAAGATTTTACCGTCACCGATTTTGCCGGTTTGCGCCGTGCGAATAATGGTATCGACACAGGTATCGACAATGTCGTCCGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATCCACCATATACTCCGCGCCGCGGTa                                                                                            >  2:1189187/1‑150 (MQ=255)
                                                             cAAAGACGAAGATTTTACCGTCACCGATTTTGCCGGTTTGCGCCGTGCGAATAATGGTATCGACACAGGTATCGACAATGTCGTCCGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATCCACCATATACTCCGCGCCGCGGTACa                                                                                          <  2:97850/150‑1 (MQ=255)
                                                                           ttACCGTCACCGATTTTGCCGGTTTGCGCCGTGCGAATAATGGTATCGACACAGGTATCGACAATGTCGTCCGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATCCACCATATACTCCGCGCCGCGGTACAGCTCGGTATGGCCt                                                                            <  1:802060/150‑1 (MQ=255)
                                                                               cGTCACCGATTTTGCCGGTTTGCGCCGTGCGAATAATGGTATCGACACAGGTATCGACAATGTCGTCCGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATCCACCATATACTCCGCGCCGCGGTACAGCTCGGTATGGCCTTTCt                                                                        <  1:775357/150‑1 (MQ=255)
                                                                                  cACCGATTTTGCCGGTTTGCGCCGTGCGAATAATGGTATCGACACAGGTATCGACAATGTCGTCCGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATCCACCATATACTCCGCGCCGCGGTACAGCTCGGTATGGCCTTTCTGGc                                                                     >  2:1363543/1‑150 (MQ=255)
                                                                                   aCCGATTTTGCCGGTTTGCGCCGTGCGAATAATGGTATCGACACAGGTATCGACAATGTCGTCCGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATCCACCATATACTCCGCGCCGCGGTACAGCTCGGTATGGCCTTTCTGGCg                                                                    <  2:17524/150‑1 (MQ=255)
                                                                                   aCCGATTTTGCCGGTTTGCGCCGTGCGAATAATGGTATCGACACAGGTATCGACAATGTCGTCCGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATCCACCATATACTCCGCGCCGCGGTACAGCTCGGTATGGCCTTTCTGGCg                                                                    <  2:22925/150‑1 (MQ=255)
                                                                                      gATTTTGCCGGTTTGCGCCGTGCGAATAATGGTATCGACACAGGTATCGACAATGTCGTCCGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATCCACCATATACTCCGCGCCGCGGTACAGCTCGGTATGGCCTTTCTGGCGAcc                                                                 >  2:172984/1‑150 (MQ=255)
                                                                                         tttGCCGGTTTGCGCCGTGCGAATAATGGTATCGACACAGGTATCGACAATGTCGTCCGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATCCACCATATACTCCGCGCCGCGGTACAGCTCGGTATGGCCTTTCTGGCGACCaaa                                                              <  2:335126/150‑1 (MQ=255)
                                                                                          ttcccGGTTTGCGCCGTGCGAATAATGGTATCGACACAGGTATCGACAATGTCGTCCGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATCCACCATATACTCCGCGCCGCGGTACAGCTCGGTATGGCCTTTCTGGCGACCAAAg                                                             >  2:810358/4‑150 (MQ=255)
                                                                                           tGCCGGTTTGCGCCGTGCGAATAATGGTATCGACACAGGTATCGACAATGTCGTCCGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATCCACCATATACTCCGCGCCGCGGTACAGCTCGGTATGGCCTTTCTGGCGACCAAAGc                                                            >  1:648886/1‑150 (MQ=255)
                                                                                                gTTTGCGCCGTGCGAATAATGGTATCGACACAGGTATCGACAATGTCGTCCGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATCCACCATATACTCCGCGCCGCGGTACAGCTCGGTATGGCCTTTCTGGCGACCAAAGCCTTTc                                                       <  2:148889/150‑1 (MQ=255)
                                                                                                 tttGCGCCGTGCGAATAATGGTATCGACACAGGTATCGACAATGTCGTCCGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATCCACCATATACTCCGCGCCGCGGTACAGCTCGGTATGGCCTTTCTGGCGACCAAAGCCTTTCa                                                      <  2:1345017/150‑1 (MQ=255)
                                                                                                 tttGCGCCGTGCGAATAATGGTATCGACACAGGTATCGACAATGTCGTCCGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATCCACCATATACTCCGCGCCGCGGTACAGCTCGGTATGGCCTTTCTGGCGACCAAAGCCTTTCa                                                      <  1:1262965/150‑1 (MQ=255)
                                                                                                     cgcCGTGCGAATAATGGTATCGACACAGGTATCGACAATGTCGTCCGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATCCACCATATACTCCGCGCCGCGGTACAGCTCGGTATGGCCTTTCTGGCGACCAAAGCCTTTCACTTc                                                  <  1:1326106/150‑1 (MQ=255)
                                                                                                        cGTGCGAATAATGGTATCGACACAGGTATCGACAATGTCGTCCGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATCCACCATATACTCCGCGCCGCGGTACAGCTCGGTATGGCCTTTCTGGCGACCAAAGCCTTTCACTTCGGt                                               <  1:129249/150‑1 (MQ=255)
                                                                                                         gTGCGAATAATGGTATCGACACAGGTATCGACAATGTCGTCCGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATCCACCATATACTCCGCGCCGCGGTACAGCTCGGTATGGCCTTTCTGGCGACCAAAGCCTTTCACTTCGGTc                                              >  1:382377/1‑150 (MQ=255)
                                                                                                            cGAATAATGGTATCGACACAGGTATCGACAATGTCGTCCGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATCCACCATATACTCCGCGCCGCGGTACAGCTCGGTATGGCCTTTCTGGCGACCAAAGCCTTTCACTTCGGTCAcc                                           >  1:234112/1‑150 (MQ=255)
                                                                                                            cGAATAATGGTATCGACACAGGTATCGACAATGTCGTCCGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATCCACCATATACTCCGCGCCGCGGTACAGCTCGGTATGGCCTTTCTGGCGACCAAAGCCTTTCACTTCGGTCAcc                                           >  1:1245890/1‑150 (MQ=255)
                                                                                                              aataatGGTATCGACACAGGTATCGACAATGTCGTCCGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATCCACCATATACTCCGCGCCGCGGTACAGCTCGGTATGGCCTTTCTGGCGACCAAAGCCTTTCACTTCGGTCACCGt                                         >  2:30855/1‑150 (MQ=255)
                                                                                                              aataatGGTATCGACACAGGTATCGACAATGTCGTCCGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATCCACCATATACTCCGCGCCGCGGTACAGCTCGGTATGGCCTTTCTGGCGACCAAAGCCTTTCACTTCGGTCACCGt                                         >  1:1133375/1‑150 (MQ=255)
                                                                                                                taatGGTATCGACACAGGTATCGACAATGTCGTCCGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATCCACCATATACTCCGCGCCGCGGTACAGCTCGGTATGGCCTTTCTGGCGACCAAAGCCTTTCACTTCGGTCACCGTCa                                       <  2:648886/150‑1 (MQ=255)
                                                                                                                    ggTATCGACACAGGTATCGACAATGTCGTCCGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATCCACCATATACTCCGCGCCGCGGTACAGCTCGGTATGGCCTTTCTGGCGACCAAAGCCTTTCACTTCGGTCACCGTCATGcc                                   >  2:272223/1‑150 (MQ=255)
                                                                                                                      tATCGACACAGGTATCGACAATGTCGTCCGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATCCACCATATACTCCGCGCCGCGGTACAGCTCGGTATGGCCTTTCTGGCGACCAAAGCCTTTCACTTCGGTCACCGTCATGCCgg                                 <  2:594411/150‑1 (MQ=255)
                                                                                                                       aTCGACACAGGTATCGACAATGTCGTCCGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATCCACCATATACTCCGCGCCGCGGTACAGCTCGGTATGGCCTTTCTGGCGACCAAAGCCTTTCACTTCGGTCACCGTCATGCCGGt                                >  1:804441/1‑150 (MQ=255)
                                                                                                                       aTCGACACAGGTATCGACAATGTCGTCCGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATCCACCATATACTCCGCGCCGCGGTACAGCTCGGTATGGCCTTTCTGGCGACCAAAGCCTTTCACTTCGGTCACCGTCATGCCGGt                                <  2:382377/150‑1 (MQ=255)
                                                                                                                           acacAGGTATCGACAATGTCGTCCGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATCCACCATATACTCCGCGCCGCGGTACAGCTCGGTATGGCCTTTCTGGCGACCAAAGCCTTTCACTTCGGTCACCGTCATGCCGGTAATa                            >  1:1219185/1‑150 (MQ=255)
                                                                                                                                    ttgagattGTCGTACCGGACGACAATCTCAATTTTCACTTTCGGCAGAAAATCCACCATATACTCCGCGCCGCGGTACAGCTCGGTATGGCCTTTCTGGCGACCAAAGCCTTTCACTTCGGTCACCGTCATGCCGGTAATACCGACTTCg                   >  1:802345/8‑150 (MQ=255)
                                                                                                                                     tgagattGTCGTACCGGACGACAATCTCAATTTTCACTTTCGGCAGAAAATCCACCATATACTCCGCGCCGCGGTACAGCTCGGTATGGCCTTTCTGGCGACCAAAGCCTTTCACTTCGGTCACCGTCATGCCGGTAATACCGACTTCgg                  >  2:1218421/7‑150 (MQ=255)
                                                                                                                                     cGACAATGTCGTCCGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATCCACCATATACTCCGCGCCGCGGTACAGCTCGGTATGGCCTTTCTGGCGACCAAAGCCTTTCACTTCGGTCACCGTCATGCCGGTAATACCGACTTCgg                  >  2:934345/1‑150 (MQ=255)
                                                                                                                                      gagattGTCGTACCGGACGACAATCTCAATTTTCACTTTCGGCAGAAAATCCACCATATACTCCGCGCCGCGGTACAGCTCGGTATGGCCTTTCTGGCGACCAAAGCCTTTCACTTCGGTCACCGTCATGCCGGTAATACCGACTTCGGc                 >  2:224257/6‑150 (MQ=255)
                                                                                                                                      gACAATGTCGTCCGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATCCACCATATACTCCGCGCCGCGGTACAGCTCGGTATGGCCTTTCTGGCGACCAAAGCCTTTCACTTCGGTCACCGTCATGCCGGTAATACCGACTTCGGc                 >  2:983064/1‑150 (MQ=255)
                                                                                                                                       aCAATGTCGTCCGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATCCACCATATACTCCGCGCCGCTGTACAGCTCGGTATGGCCTTTCTGGCGACCAAAGCCTTTCACTTCGGTCACCGTCATGCCGGTAATACCGACTTCGGcc                <  1:526936/150‑1 (MQ=255)
                                                                                                                                       aCAATGTCGTCCGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATCCACCATATACTCCGCGCCGCGGTACAGCTCGGTATGGCCTTTCTGGCGACCAAAGCCTTTCACTTCGGTCACCGTCATGCCGGTAATACCGACTTCGGcc                <  1:696336/150‑1 (MQ=255)
                                                                                                                                          aTGTCGTCCGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATCCACCATATACTCCGCGCCGCGGTACAGCTCGGTATGGCCTTTCTGGCGACCAAAGCCTTTCACTTCGGTCACCGTCATGCCGGTAATACCGACTTCGGCCAGt             <  1:1172856/150‑1 (MQ=255)
                                                                                                                                            gtcgtcCGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATCCACCATATACTCCGCGCCGCGGTACAGCTCGGTATGGCCTTTCTGGCGACCAAAGCCTTTCACTTCGGTCACCGTCATGCCGGTAATACCGACTTCGGCCAGTGc           >  2:742111/1‑150 (MQ=255)
                                                                                                                                              cgtcCGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATCCACCATATACTCCGCGCCGCGGTACAGCTCGGTATGGCCTTTCTGGCGACCAAAGCCTTTCa                                                      <  2:331962/105‑1 (MQ=255)
                                                                                                                                              cgtcCGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATCCACCATATACTCCGCGCCGCGGTACAGCTCGGTATGGCCTTTCTGGCGACCAAAGCCTTTCa                                                      >  1:331962/1‑105 (MQ=255)
                                                                                                                                                 ccGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATCCACCATATACTCCGCGCCGCGGTACAGCTCGGTATGGCCTTTCTGGCGACCAAAGCCTTTCACTTCGGTCACCGTCATGCCGGTAATACCGACTTCGGCCAGTGCTTcgc      >  2:1125432/1‑150 (MQ=255)
                                                                                                                                                  cGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATCCACCATATACTCCGCGCCGCGGTACAGCTCGGTATGGCCTTTCTGGCGACCAAAGCCTTTCACTTCGGTCACCGTCATGCCGGTAATACCGACTTCGGCCAGTGCTTcgcg     <  1:742111/150‑1 (MQ=255)
                                                                                                                                                  cGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATCCACCATATACTCCGCGCCGCGGTACAGCTCGGTATGGCCTTTCTGGCGACCAAAGCCTTTCACTTCGGTCACCGTCATGCCGGTAATACCGACTTCGGCCAGTGCTTcgcg     <  1:1334342/150‑1 (MQ=255)
                                                                                                                                                    gTACGACAATCTCAATTTTCACTTTCGGCAGAAAATCCACCATATACTCCGCGCCGCGGTACAGCTCGGTATGGCCTTTCTGGCGACCAAAGCCTTTCACTTCGGTCACCGTCATGCCGGTAATACCGACTTCGGCCAGTGCTTCGCGGa   <  2:1245890/150‑1 (MQ=255)
                                                                                                                                                     tACGACAATCTCAATTTTCACTTTCGGCAGAAAATCCACCATATACTCCGCGCCGCGGTACAGCTCGGTATGGCCTTTCTGGCGACCAAAGCCTTTCACTTCGGTCACCGTCATGCCGGTAATACCGACTTCGGCCAGTGCTTCGCGGAc  <  2:896057/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCGGTTTAAATTGCCGCGTCGTCCTCTTCACCGGTACGGATGCGAATGACCCGTGCCACGTCAAAGACGAAGATTTTACCGTCACCGATTTTGCCGGTTTGCGCCGTGCGAATAATGGTATCGACACAGGTATCGACAATGTCGTCCGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATCCACCATATACTCCGCGCCGCGGTACAGCTCGGTATGGCCTTTCTGGCGACCAAAGCCTTTCACTTCGGTCACCGTCATGCCGGTAATACCGACTTCGGCCAGTGCTTCGCGGAC  >  NC_000913/2687065‑2687363

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: ‑