Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,833,205 G→A A260T (GCC→ACC)  xthA → exonuclease III

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,833,2050GA98.4% 222.7 / NA 65A260T (GCC→ACC) xthAexonuclease III
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (30/33);  total (30/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

CTCATGGTTTGATTACCGCTCAAAAGGTTTTGACGATAACCGTGGTCTGCGCATCGACCTGCTGCTCGCCAGCCAACCGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACGCCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTGAGAAAAATCAACAAACTGTCAGTAATGATTTGTTGCCTGCCGTCCTTTGTTATACCGTCTCTGCGTTTTT  >  NC_000913/1833057‑1833350
                                                                                                                                                    |                                                                                                                                                 
ctcATGGTTTGATTACCGCTCAAAAGGTTTTGACGATAACCGTGGTCTGCGCATCGACCTGCTGCTCGCCAGCCAACCGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACAcc                                                                                                                                                 >  2:588116/1‑151 (MQ=255)
 tcATGGTTTGATTACCGCTCAAAAGGATTTGACGATAACCGTGGTCTGCGCATCGACCTGCTGCTCGCCAGCCAACCGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACAccc                                                                                                                                                >  1:825084/1‑151 (MQ=255)
   aTGGTTTGATTACCGCTCAAAAGGTTTTGACGATAACCGTGGTCTGCGCATCGACCTGCTGCTCGCCAGCCAACCGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACAccccc                                                                                                                                              <  2:248097/151‑1 (MQ=255)
     ggTTTGATTACCGCTCAAAAGGTTTTGACGATAACCGTGGTCTGCGCATCGACCTGCTGCTCGCCAGCCAACCGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGt                                                                                                                                            <  1:1097671/151‑1 (MQ=255)
          gATTACCGCTCAAAAGGTTTTGACGATAACCGTGGTCTGCGCATCGACCTGCTGCTCGCCAGCCAACCGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTggg                                                                                                                                       >  1:87575/1‑151 (MQ=255)
             tACCGCTCAAAAGGTTTTGACGATAACCGTGGTCTGCGCATCGACCTGCTGCTCGCCAGCCAACCGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCg                                                                                                                                     >  2:648228/1‑150 (MQ=255)
                   tCAAAAGGTTTTGACGATAACCGTGGTCTGCGCATCGACCTGCTGCTCGCCAGCCAACCGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTcc                                                                                                                              >  1:133383/1‑151 (MQ=255)
                             ttGACGATAACCGTGGTCTGCGCATCGACCTGCTGCTCGCCAGCCAACCGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAAt                                                                                                                     <  2:21384/150‑1 (MQ=255)
                              tGACGATAACCGTGGTCTGCGCATCGACCTGCTGCTCGCCAGCCAACCGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAAttt                                                                                                                   >  1:709275/1‑151 (MQ=255)
                               gACGATAACCGTGGTCTGCGCATCGACCTGCTGCTCGCCAGCCAACCGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAAttt                                                                                                                   >  2:1156084/1‑150 (MQ=255)
                                   aTAACCGTGGTCTGCGCATCGACCTGCTGCTCGCCAGCCAACCGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTagc                                                                                                                >  2:153843/1‑149 (MQ=255)
                                     aaCCGTGGTCTGCGCATCGACCTGCTGCTCGCCAGCCAACCGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTagcagc                                                                                                             >  1:344140/1‑150 (MQ=255)
                                     aaCCGTGGTCTGCGCATCGACCTGCTGCTCGCCAGCCAACCGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGct                                                                                                            >  2:515562/1‑151 (MQ=255)
                                        cGTGGTCTGCGCATCGACCTGCTGCTCGCCAGCCAACCGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCtctc                                                                                                         >  1:662494/1‑151 (MQ=255)
                                         gTGGTCTGCGCATCGACCTGCTGCTCGCCAGCCAACCGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTcc                                                                                                        >  2:527453/1‑151 (MQ=255)
                                                gcgcATCGACCTGCTGCTCGCCAGCCAACCGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCt                                                                                                   <  1:259021/149‑1 (MQ=255)
                                                gcgcATCGACCTGCTGCTCGCCAGCCAACCGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTc                                                                                                  >  1:819176/1‑150 (MQ=255)
                                                gcgcATCGACCTGCTGCTCGCCAGCCAACCGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCa                                                                                                 <  2:139089/151‑1 (MQ=255)
                                                      cGACCTGCTGCTCGCCAGCCAACCGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTgg                                                                                           <  1:700211/151‑1 (MQ=255)
                                                      cGACCTGCTGCTCGCCAGCCAACCGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTg                                                                                            <  1:14603/150‑1 (MQ=255)
                                                            gctgctCGCCAGCCAACCGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCa                                                                                       <  1:113798/149‑1 (MQ=255)
                                                            gctgctCGCCAGCCAACCGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCa                                                                                       <  1:192750/149‑1 (MQ=255)
                                                            gctgctCGCCAGCCAACCGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCa                                                                                       <  1:1115431/149‑1 (MQ=255)
                                                                ctCGCCAGCCAACCGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACGCCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGgaga                                                                                  <  1:1079757/150‑1 (MQ=255)
                                                                cccGCCAGCCAACCGCTGGCAGAATTTTGCGTAGAAACCTGCATCGACTATGAAATCCGCAGCATGGAAAACCCGTCCGATCACCCCCCCGTCCGGGCGATCTTCCGCCCCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGgaga                                                                                  <  1:766786/148‑1 (MQ=255)
                                                                  cGCCAGCCAACCGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGaa                                                                                 <  2:587640/149‑1 (MQ=255)
                                                                  cGCCAGCCAACCGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAAtt                                                                               <  2:1134983/151‑1 (MQ=255)
                                                                   gccagccaACCGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTa                                                                              <  2:546752/151‑1 (MQ=255)
                                                                       gccaACCGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCt                                                                          <  2:422244/151‑1 (MQ=255)
                                                                        ccaACCGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCtt                                                                         >  1:588585/1‑151 (MQ=255)
                                                                        acaACCGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCtt                                                                         >  2:212547/2‑151 (MQ=255)
                                                                         caACCGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTg                                                                        <  2:782910/151‑1 (MQ=255)
                                                                         caACCGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTg                                                                        <  1:90821/151‑1 (MQ=255)
                                                                          aaCCGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTga                                                                       <  1:1094970/151‑1 (MQ=255)
                                                                            ccGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTgaga                                                                     <  1:74896/151‑1 (MQ=255)
                                                                             cGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTgaga                                                                     >  2:94433/1‑150 (MQ=255)
                                                                              gCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTGAGaaa                                                                   <  1:610899/151‑1 (MQ=255)
                                                                              gCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTGAGaa                                                                    >  1:20689/1‑150 (MQ=255)
                                                                                tGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTGAGaaaa                                                                  >  2:1172796/1‑150 (MQ=255)
                                                                                    aGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTGAGAAAAATca                                                              >  1:725901/1‑150 (MQ=255)
                                                                                          ttGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTGAGAAAAATCAACAAACt                                                       >  2:164131/1‑151 (MQ=255)
                                                                                           tGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTGAGAAAAATCAACAAACt                                                       <  1:796645/150‑1 (MQ=255)
                                                                                             cGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTGAGAAAAATCAACAAACTGTc                                                    <  1:71097/151‑1 (MQ=255)
                                                                                                aGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTGAGAAAAATCAACAAACTGTCAg                                                  <  2:442094/150‑1 (MQ=255)
                                                                                                 gAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTGAGAAAAATCAACAAACTGTCAGTa                                                <  2:345738/151‑1 (MQ=255)
                                                                                                     ccGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTGAGAAAAATCAACAAACTGTCAGTAATg                                             >  1:557198/1‑150 (MQ=255)
                                                                                                        gCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTGAGAAAAATCAACAAACTGTCAGTAATGAttt                                         <  1:935240/151‑1 (MQ=255)
                                                                                                            cGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTGAGAAAAATCAACAAACTGTCAGTAATGATttgtt                                      <  2:87599/150‑1 (MQ=255)
                                                                                                                 aTGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTGAGAAAAATCAACAAACTGTCAGTAATGATTTGTtgcctgc                                >  2:281890/1‑151 (MQ=255)
                                                                                                                 aTGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTGAGAAAAATCAACAAACTGTCAGTAATGATTTGTtgcctg                                 <  2:1174834/150‑1 (MQ=255)
                                                                                                                 aTGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTGAGAAAAATCAACAAACTGTCAGTAATGATTTGTtgcctg                                 <  2:191539/150‑1 (MQ=255)
                                                                                                                   gAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTGAGAAAAATCAACAAACTGTCAGTAATGATTTGTTGCCTGCCg                              <  1:592560/151‑1 (MQ=255)
                                                                                                                          gcagcaTGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTGAGAAAAATCAACAAACTGTCAGTAATGATTTGTTGCCTGCCGTCCTTTg                       >  2:346974/1‑151 (MQ=255)
                                                                                                                           cagcaTGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTGAGAAAAATCAACAAACTGTCAGTAATGATTTGTTGCCTGCCGTCCTTTGt                      >  2:126837/1‑151 (MQ=255)
                                                                                                                           cagcaTGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTGAGAAAAATCAACAAACTGTCAGTAATGATTTGTTGCCTGCCGTCCTTTGt                      >  2:721733/1‑151 (MQ=255)
                                                                                                                            agcaTGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTGAGAAAAATCAACAAACTGTCAGTAATGATTTGTTGCCTGCCGTCCTTTGt                      <  2:825271/150‑1 (MQ=255)
                                                                                                                               aTGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTGAGAAAAATCAACAAACTGTCAGTAATGATTTGTTGCCTGCCTTCCTTTGTtat                   <  2:981462/150‑1 (MQ=255)
                                                                                                                                  gAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTGAGAAAAATCAACAAACTGTCAGTAATGATTTGTTGCCTGCCGTCCTTTGTTATACCg               >  2:785633/1‑151 (MQ=255)
                                                                                                                                  gAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTGAGAAAAATCAACAAACTGTCAGTAATGATTTGTTGCCTGCCGTCCTTTGTTATACCg               >  2:791999/1‑151 (MQ=255)
                                                                                                                                    aaaaCCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTGAGAAAAATCAACAAACTGTCAGTAATGATTTGTTGCCTGCCGTCCTTTGTTATACCGtc             <  1:1162520/151‑1 (MQ=255)
                                                                                                                                     aaaCCGTCCGATCCCACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTGAGAAAAATCAACAAACTGTCAGTAATGATTTGTTGCCTGCCGTCCTTTGTTATACCGtct            <  1:200673/151‑1 (MQ=255)
                                                                                                                                       aCCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTGAGAAAAATCAACAAACTGTCAGTAATGATTTGTTGCCTGCCGTCCTTTGTTATACCGTctct          <  1:527321/151‑1 (MQ=255)
                                                                                                                                         cGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTGAGAAAAATCAACAAACTGTCAGTAATGATTTGTTGCCTGCCGTCCTTTGTTATACCGTCTCTGc        >  2:1154674/1‑151 (MQ=255)
                                                                                                                                               aTCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTGAGAAAAATCAACAAACTGTCAGTAATGATTTGTTGCCTGCCGTCCTTTGTTATACCGTCTCTGCGttttt  >  1:244054/1‑151 (MQ=255)
                                                                                                                                                tCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTGAGAAAAATCAACAAACTGTCAGTAATGATTTGTTGCCTGCCGTCCTTTGTTATACCGTCTCTGCGttttt  >  1:1086818/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
CTCATGGTTTGATTACCGCTCAAAAGGTTTTGACGATAACCGTGGTCTGCGCATCGACCTGCTGCTCGCCAGCCAACCGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACGCCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTGAGAAAAATCAACAAACTGTCAGTAATGATTTGTTGCCTGCCGTCCTTTGTTATACCGTCTCTGCGTTTTT  >  NC_000913/1833057‑1833350

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: