Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,446,786 T→A 7.4% E166D (GAA→GAT) ‡ feaR ← transcriptional activator for tynA and feaB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,446,7860TA 7.4% 280.1 / 10.0 104E166D (GAA→GAT) ‡feaRtranscriptional activator for tynA and feaB
Reads supporting (aligned to +/- strand):  ref base T (53/43);  new base A (5/3);  total (58/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.93e-01

Only 100 of 106 total aligned reads displayed.

TAATATCTCTTCGCGAATATTATCGTCTATCAACGTAACCACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGTTCCAGCAGAGTGCGTGGCAAAAGTAATGAAA  >  NC_000913/1446642‑1446934
                                                                                                                                                |                                                                                                                                                    
tAATATCTCTTCGCGAATATTATCGTCTATCAACGTAACCAATTTTTGAAACTGACGTTCACGAAGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGCCACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGa                                                                                                                                                 >  1:460864/1‑150 (MQ=255)
 aaTATCTCTTCGCGAATATTATCGTCTATCAACGTAACCACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGaa                                                                                                                                                <  2:1314096/150‑1 (MQ=255)
    atCTCTTCTCGAATATTATAGTCTATCAACGTAACCACTTTTTGAAACTGACGTTCACTACGAGGTTGAACAGATTACCGCTGATGAAGTACCGGTCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGt                                                                                                                                             >  2:392847/1‑150 (MQ=255)
    atCTCTTCGCGAATATTATCGTCTATCAACGTAACCACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGt                                                                                                                                             >  2:663104/1‑150 (MQ=255)
       tctTCGCGAATATTATCGTCTATCAACGTAACCACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGAAACAGAAAGTgaa                                                                                                                                          <  1:1484498/150‑3 (MQ=255)
       tctTCGCGAATATTATCGTCTATCAACGTAACCACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGAAACAGAAAGTgaa                                                                                                                                          <  1:1346916/150‑3 (MQ=255)
          tCGCGAATATTATCGTCTATCAACGTAACCACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGa                                                                                                                                       >  1:1290297/1‑150 (MQ=255)
          tCGCGAATATTATCGTCTATCAACGTAACCACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCAGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTGTCTGTTGCAGAAAGTGCCGGa                                                                                                                                       >  2:1106875/1‑150 (MQ=255)
                attttatCGTCTATCAACGTAACCACTTTTTGAATCTGACGTTCACGACGAGGTTGAACAGATTCCCGCAGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTc                                                                                                                                 <  1:777508/147‑1 (MQ=255)
                atatTATCGTCTATCAACGTAACCACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCATGattatta                                                                                                                                 >  1:780690/1‑149 (MQ=255)
                 tattatCGTCTATCAACGTAACCACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCa                                                                                                                                <  1:1225812/150‑1 (MQ=255)
                   ttatCGTATATCAACGTAACCACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATg                                                                                                                              <  2:625256/150‑1 (MQ=255)
                     atCGTCTATCAACGTAACCACTTTTTGAAACTGACGTTCACGACGAGGTTGATCAGATTCCCGATGATGAAGTTCCGGGCGCATCAGACACACCATAGAATGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGct                                                                                                                            <  1:84266/150‑1 (MQ=255)
                        gTCTATCAACGTAACCACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCtctc                                                                                                                         >  2:1666657/1‑150 (MQ=255)
                           tATCAACGTAACCACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTg                                                                                                                      >  1:918454/1‑150 (MQ=255)
                            aTCAACTTAACCACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGTAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGt                                                                                                                     <  1:1161077/150‑1 (MQ=255)
                            aTCAACGTAACCACTTTTTGAAACTGACGTTCACGACGTGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGt                                                                                                                     >  2:727384/1‑150 (MQ=255)
                             tCAACGTAACCACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCACCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTa                                                                                                                    >  1:1700035/1‑150 (MQ=255)
                                aCGTAACCACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACa                                                                                                                 <  2:1498747/150‑1 (MQ=255)
                                   tAACCACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGCCACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGAGTATTAATGCTCTCCTGTAACAGGc                                                                                                              >  2:282870/1‑150 (MQ=255)
                                   tAACCACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGc                                                                                                              >  2:1357986/1‑150 (MQ=255)
                                   tAACCACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCAGGAATATTCATGCTCTCCTGTAACAGGc                                                                                                              >  2:1629269/1‑150 (MQ=255)
                                    aaCCACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACTATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCg                                                                                                             >  1:1182184/1‑150 (MQ=255)
                                      ccACTTTTTGAAACTGACGTTCACGACGAGGTAGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTCTTTCAGAAAGTGACGGATTAATCATGCTCTCCTGTAACTGGCGAt                                                                                                           >  1:159124/1‑150 (MQ=255)
                                       cACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGATTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATg                                                                                                          <  2:1047141/150‑1 (MQ=255)
                                         cTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGAc                                                                                                        >  1:1690407/1‑150 (MQ=255)
                                           ttttGAAACTGACGTTCACGACGAGGTTGAACTGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGATTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTg                                                                                                      >  2:842556/1‑150 (MQ=255)
                                              tGAAACTGACGTTCACGACTAGGTTGAACAGATTCCAGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGt                                                                                                   <  2:1324677/150‑1 (MQ=255)
                                               gAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGtt                                                                                                  >  2:1388455/1‑150 (MQ=255)
                                                aaaCTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTg                                                                                                 >  2:797628/1‑150 (MQ=255)
                                                    tGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCAcc                                                                                             <  1:1222579/150‑1 (MQ=255)
                                                      aCGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATg                                                                                                          <  1:1102958/135‑1 (MQ=255)
                                                      aCGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATg                                                                                                          >  2:1102958/1‑135 (MQ=255)
                                                       cGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATg                                                                                          <  2:664075/150‑1 (MQ=255)
                                                       cGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATg                                                                                          >  2:683982/1‑150 (MQ=255)
                                                       cGTTCACGACGAGGTAGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATg                                                                                          >  1:937301/1‑150 (MQ=255)
                                                         ttCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATggg                                                                                        <  1:1645815/150‑1 (MQ=255)
                                                          tCACGACGAGGTTGAACAGATTCACGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGACTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGATCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGt                                                                                       <  1:107828/150‑1 (MQ=255)
                                                            acgacgAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCAGTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTaa                                                                                     <  2:801540/150‑1 (MQ=255)
                                                             cgacgaGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAg                                                                                    >  2:132283/1‑150 (MQ=255)
                                                             cgacgaGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGGTGCACCATGGGTAAg                                                                                    >  1:1263716/1‑150 (MQ=255)
                                                                    gTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCCCCATGGGTAAGTCAGCGt                                                                             >  2:302638/1‑150 (MQ=255)
                                                                      tGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCATCATGGGTAAGTCAGCGTcc                                                                           >  2:1538555/1‑150 (MQ=255)
                                                                      tGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTcc                                                                           <  2:603072/150‑1 (MQ=255)
                                                                        aaCAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTATGTATCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACATGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAg                                                                         >  2:1241781/1‑150 (MQ=255)
                                                                          cAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTc                                                                       >  2:1348063/1‑150 (MQ=255)
                                                                           aGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCt                                                                      <  1:933853/150‑1 (MQ=255)
                                                                            gATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCtt                                                                     >  2:258801/1‑150 (MQ=255)
                                                                              ttCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTc                                                                   <  2:1035612/150‑1 (MQ=255)
                                                                                cccGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAACGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTg                                                                 >  1:1083711/1‑150 (MQ=255)
                                                                                cccGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTg                                                                 >  1:19036/1‑150 (MQ=255)
                                                                                cccGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTg                                                                 <  2:1157060/150‑1 (MQ=255)
                                                                                  cGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTgcg                                                               >  2:1186445/1‑150 (MQ=255)
                                                                                   gCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTgcgc                                                              <  1:392847/150‑1 (MQ=255)
                                                                                     tgatgaAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAg                                                            <  2:918454/150‑1 (MQ=255)
                                                                                     tgatgaAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGAGTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAg                                                            >  2:694294/1‑150 (MQ=255)
                                                                                        tgaAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATa                                                         >  2:1244343/1‑150 (MQ=255)
                                                                                            gTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGtt                                                     <  1:1357986/150‑1 (MQ=255)
                                                                                               ccGGGCGCAGTAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTACAGTCTTTCTGCGCAGATAGGttttt                                                  >  2:881928/1‑150 (MQ=255)
                                                                                                cGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTg                                                 <  2:1290297/150‑1 (MQ=255)
                                                                                                cGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTg                                                 <  2:580841/150‑1 (MQ=255)
                                                                                                  ggCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGAt                                               <  1:258801/150‑1 (MQ=255)
                                                                                                  ggCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACAATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGAt                                               >  2:1661564/1‑150 (MQ=255)
                                                                                                   gcgcAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATg                                              <  1:727384/150‑1 (MQ=255)
                                                                                                    cgcAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATgg                                             <  1:683982/150‑1 (MQ=255)
                                                                                                    cgcAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATgg                                             >  1:1540698/1‑150 (MQ=255)
                                                                                                     gcagcaGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATggg                                            <  2:1709335/150‑1 (MQ=255)
                                                                                                     gcagcaGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATggg                                            >  1:1048779/1‑150 (MQ=255)
                                                                                                       agcagACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGa                                          >  2:361388/1‑150 (MQ=255)
                                                                                                           tacacacCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGTTTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAat                                      <  2:1331697/149‑1 (MQ=255)
                                                                                                           gACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAat                                      >  1:1330783/1‑150 (MQ=255)
                                                                                                               cacCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGAGTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTg                                  <  1:282237/150‑1 (MQ=255)
                                                                                                                acCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGt                                 <  1:797628/150‑1 (MQ=255)
                                                                                                                  cATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGTTc                               <  1:663104/150‑1 (MQ=255)
                                                                                                                   aTCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGTTcc                              <  2:19036/150‑1 (MQ=255)
                                                                                                                   aTCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGTTcc                              <  1:26817/150‑1 (MQ=255)
                                                                                                                       ccTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTTAGTCAGCGTCCTGGCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGTTCcatca                          >  2:840382/1‑150 (MQ=255)
                                                                                                                       ccTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGTTCcagca                          >  2:1100015/1‑150 (MQ=255)
                                                                                                                       ccTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGTTCcagca                          >  1:1242333/1‑150 (MQ=255)
                                                                                                                          gTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGTTCCAGCagag                       <  1:573594/150‑1 (MQ=255)
                                                                                                                             gcgcAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGTTCCAGCAGAGTGc                    >  1:1349506/1‑150 (MQ=255)
                                                                                                                                cAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGTTCCAGCAGAGTGCGTg                 >  1:1607333/1‑150 (MQ=255)
                                                                                                                                cAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGTTCCAGCAGAGTGCGTg                 <  2:146709/150‑1 (MQ=255)
                                                                                                                                cAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGTTCCAGCAGAGTGCGTg                 >  1:1492085/1‑150 (MQ=255)
                                                                                                                                 aGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGTTCCAGCAGAGTGCGTgg                <  1:1666657/150‑1 (MQ=255)
                                                                                                                                 aGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGTTCCAGCAGAGTGCGTgg                <  1:132283/150‑1 (MQ=255)
                                                                                                                                  gCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGTTCCAGCAGAGTGCGTGGc               >  2:1373961/1‑150 (MQ=255)
                                                                                                                                  ccggcACTTTCTGAAACAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGTTCCAGCAGAGTGCGTGGc               >  2:795842/5‑150 (MQ=255)
                                                                                                                                   cggcACTTTCTGAAACAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGTTCCAGCAGAGTGCGTGGCa              >  2:538131/4‑150 (MQ=255)
                                                                                                                                   cggcACTTTCTGAAACAGAAAGTGCCAGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGTTCCAGCAGAGTGCGTGGCa              >  2:1372335/4‑150 (MQ=255)
                                                                                                                                     tCACTTTCTGTTTCATAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGTTCCAGCAGAGTGCGTGGCaaa            <  2:159124/150‑1 (MQ=255)
                                                                                                                                      cACTTTCTGAAACAGAAAGTGCCTGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGTTCCAGCAGAGTGCGTGGCaaaa           <  1:1372335/150‑1 (MQ=255)
                                                                                                                                      cACTTTCTGAAACAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGTTCCAGCAGAGTGCGTGGCaaaa           >  2:1259562/1‑150 (MQ=255)
                                                                                                                                        atttctgtttctgAAAGTGCCGGATAATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGTTCCAGCAGAGTGCGTGGCAAAAGt         <  2:460864/149‑1 (MQ=255)
                                                                                                                                          ttCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGTTCCAGCAGAGTGCGTGGCAAAAGTaa       >  2:741145/1‑150 (MQ=255)
                                                                                                                                          ttCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGTTCCAGCAGAGTGCGTGGCAAAAGTaa       >  2:329485/1‑150 (MQ=255)
                                                                                                                                           tCTGTTTCAGAAAGTGCCGGATTATTAATGCTCTACTGTAACAGGCGTTGACTGAGTTGCTCCATGGGTAAGTCAGCGTACAGTCTTTCTGAGCAGATAGGTTTTTGATGTGGAAAATATTTTTCCAGCAGATTGCGTGGCAAAAGTAAt      >  2:854312/1‑150 (MQ=255)
                                                                                                                                            cTGTTTCCGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGGTTTTGATGGGGAAAATATTGTTCCAGCAGAGTGCGTGGCAAAAGTAATg     <  1:364056/150‑1 (MQ=255)
                                                                                                                                              gTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGTTCCAGCAGAGTGCGTGGCAAAAGTAATGaa   >  2:686659/1‑150 (MQ=255)
                                                                                                                                              gTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGTTCCAGCAGAGTGCGTGGCAAAAGTAATGaa   >  2:864377/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
TAATATCTCTTCGCGAATATTATCGTCTATCAACGTAACCACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGTTCCAGCAGAGTGCGTGGCAAAAGTAATGAAA  >  NC_000913/1446642‑1446934

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