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.5% 236.9 / NA 67A260T (GCC→ACC) xthAexonuclease III
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (29/37);  minor base C (0/1);  total (29/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 = 1.00e+00

CATGGTTTGATTACCGCTCAAAAGGTTTTGACGATAACCGTGGTCTGCGCATCGACCTGCTGCTCGCCAGCCAACCGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACGCCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTGAGAAAAATCAACAAACTGTCAGTAATGATTTGTTGCCTGCCGTCCTTTGTTATACCGTCTCTGCGTTTTTAG  >  NC_000913/1833059‑1833352
                                                                                                                                                  |                                                                                                                                                   
cATGGTTTGATTACCGCTCAAAAGGTTTTGACGATAACCGTGGTCTGCGCATCGACCTGCTGCTCGCCAGCCAACCGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACAcccc                                                                                                                                                 <  2:830618/151‑1 (MQ=255)
         aTTACCGCTCAAAAGGTTTTGACGATAACCGTGGTCTGCGCATCGACCTGCTGCTCGCCAGCCAACCGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTggg                                                                                                                                         <  1:355825/150‑1 (MQ=255)
           tACCGCTCAAAAGGTTTTGACGATAACCGTGGTCTGCGCATCGACCTGCTGCTCGCCAGCCAACCGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCg                                                                                                                                       >  1:1020755/1‑150 (MQ=255)
           tACCGCTCAAAAGGTTTTGACGATAACCGTGGTCTGCGCATCGACCTGCTGCTCGCCAGCCAACCGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCg                                                                                                                                       <  2:1181747/150‑1 (MQ=255)
                  cAAAAGGTTTTGACGATAACCGTGGTCTGCGCATCGACCTGCTGCTCGCCAGCCAACCGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTccg                                                                                                                               <  1:188924/151‑1 (MQ=255)
                   aaaaGGTTTTGACGATAACCGTGGTCTGCGCATCGACCTGCTGCTCGCCAGCCAACCGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCtt                                                                                                                                  <  2:183633/147‑1 (MQ=255)
                       ggTTTTGACGATAACCGTGGTCTGCGCATCGACCTGCTGCTCGCCAGCCAACCGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCt                                                                                                                          >  1:179603/1‑151 (MQ=255)
                       ggTTTTGACGATAACCGTGGTCTGCGCATCGACCTGCTGCTCGCCAGCCAACCGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCt                                                                                                                          <  1:311255/151‑1 (MQ=255)
                             gACGATAACCGTGGTCTGCGCATCGACCTGCTGCTCGCCAGCCAACCGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAAttt                                                                                                                     >  1:427347/1‑150 (MQ=255)
                             gACGATAACCGTGGTCTGCGCATCGACCTGCTGCTCGCCAGCCAACCGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTa                                                                                                                    <  2:647310/151‑1 (MQ=255)
                                gATAACCGTGGTCTGCGCATCGACCTGCTGCTCGCCAGCCAACCGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTagc                                                                                                                  >  1:590342/1‑150 (MQ=255)
                                   aaCCGTGGTCTGCGCATCGACCTGCTGCTCGCCAGCCAACCGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGct                                                                                                              >  1:485538/1‑151 (MQ=255)
                                       gTGATCTGCGCATCGACCTGCTGCTCGCCAGCCAACCGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTcc                                                                                                          <  1:1014870/151‑1 (MQ=255)
                                        tGGTCTGCGCATCGACCTGCTGCTCGCCAGCCAACCGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCt                                                                                                         >  2:710363/1‑151 (MQ=255)
                                         ggTCTGCGCATCGACCTGCTGCTCGCCAGCCAACCGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCt                                                                                                         <  2:14824/150‑1 (MQ=255)
                                          gTCTGCGCATCGACCTGCTGCTCGCCAGCCAACCGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTgg                                                                                                       <  2:228732/151‑1 (MQ=255)
                                                 cATCGACCTGCTGCTCGCCAGCCAACCGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAAc                                                                                                >  2:985115/1‑151 (MQ=255)
                                                 cATCGACCTGCTGCTCGCCAGCCAACCGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAAc                                                                                                <  1:64989/151‑1 (MQ=255)
                                                  aTCGACCTGCTGCTCGCCAGCCAACCGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACt                                                                                               >  1:1160469/1‑151 (MQ=255)
                                                    cGACCTGCTGCTCGCCAGCCAACCGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTgg                                                                                             >  2:45767/1‑151 (MQ=255)
                                                     gACCTGCTGCTCGCCAGCCAACCGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACCCCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTggg                                                                                            <  1:709037/151‑1 (MQ=255)
                                                     gACCTGCTGCTCGCCAGCCAACCGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTggg                                                                                            >  1:1017327/1‑151 (MQ=255)
                                                             gctCGCCAGCCAACCGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGgag                                                                                     <  1:144040/150‑1 (MQ=255)
                                                             gctCGCCAGCCAACCGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGga                                                                                      <  1:435268/149‑1 (MQ=255)
                                                                   cagccaACCGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAAc                                                                              >  1:777369/1‑151 (MQ=255)
                                                                     gccaACCGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCt                                                                            >  1:842603/1‑151 (MQ=255)
                                                                      ccaACCGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCtt                                                                           <  1:443722/151‑1 (MQ=255)
                                                                       caACCGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTg                                                                          <  1:1098930/151‑1 (MQ=255)
                                                                        aaCCGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTga                                                                         <  2:805093/151‑1 (MQ=255)
                                                                        aaCCGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTg                                                                          <  2:353481/150‑1 (MQ=255)
                                                                        aaCCGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTg                                                                          >  1:847484/1‑150 (MQ=255)
                                                                           cGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTgaga                                                                       >  2:141358/1‑150 (MQ=255)
                                                                           cGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTGAGaa                                                                      <  1:521873/151‑1 (MQ=255)
                                                                            gCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTGAGaa                                                                      >  1:671391/1‑150 (MQ=255)
                                                                              tGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTGAGaaaa                                                                    >  2:645848/1‑150 (MQ=255)
                                                                               ggCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTGAGaaaaa                                                                   <  1:507852/150‑1 (MQ=255)
                                                                               ggCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTGAGaaaa                                                                    <  1:53070/149‑1 (MQ=255)
                                                                               ggCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTGAGAAAAAt                                                                  <  1:1173476/151‑1 (MQ=255)
                                                                                     aTGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTGAGAAAAATcaacaa                                                            <  2:210307/151‑1 (MQ=255)
                                                                                        ttGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTGAGTAAAATCAACAAAc                                                          >  2:1150273/1‑150 (MQ=255)
                                                                                         tGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTGAGAAAAATCAACAAACt                                                         <  2:1111383/150‑1 (MQ=255)
                                                                                         tGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTGAGAAAAATCAACAAACTg                                                        >  1:1158331/1‑151 (MQ=255)
                                                                                          gCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTGAGAAAAATCAACAAACTGt                                                       <  2:90844/151‑1 (MQ=255)
                                                                                            gTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTGAGAAAAATCAACAAACTGTCa                                                     <  1:834593/151‑1 (MQ=255)
                                                                                             tAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTGAGAAAAATCAACAAACTGTc                                                      >  1:712523/1‑149 (MQ=255)
                                                                                              aGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTGAGAAAAATCAACAAACTGTCAGt                                                   >  1:643626/1‑151 (MQ=255)
                                                                                              aGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTGAGAAAAATCAACAAACTGTCAGt                                                   >  1:879475/1‑151 (MQ=255)
                                                                                              aGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTGAGAAAAATCAACAAACTGTCAGt                                                   >  1:933016/1‑151 (MQ=255)
                                                                                               gAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTGAGAAAAATCAACAAACTGTCAGTa                                                  <  2:1229624/151‑1 (MQ=255)
                                                                                                aaaCCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTGAGAAAAATCAACAAACTGTCAGTa                                                  >  2:1166253/1‑150 (MQ=255)
                                                                                                  aCCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTGAGAAAAATCAACAAACTGTCAGTAAt                                                >  2:326763/1‑150 (MQ=255)
                                                                                                  aCCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTGAGAAAAATCAACAAACTGTCAGTAATg                                               <  2:442522/151‑1 (MQ=255)
                                                                                                   ccGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTGAGAAAAATCAACAAACTGTCAGTAATg                                               <  2:68816/150‑1 (MQ=255)
                                                                                                    cGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTGAGAAAAATCAACAAACTGTCAGTAATGAt                                             <  2:1231626/151‑1 (MQ=255)
                                                                                                        aTCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTGAGAAAAATCAACAAACTGTCAGTAATGATttgt                                         >  2:102798/1‑151 (MQ=255)
                                                                                                          cGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTGAGAAAAATCAACAAACTGTCAGTAATGATttgttg                                       <  1:776756/151‑1 (MQ=255)
                                                                                                            aCTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTGAGAAAAATCAACAAACTGTCAGTAATGATTTGTtgcc                                     >  2:844633/1‑151 (MQ=255)
                                                                                                                      ccGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTGAGAAAAATCAACAAACTGTCAGTAATGATTTGTTGCCTGCCGTCCtt                           <  1:906911/151‑1 (MQ=255)
                                                                                                                        gcagcaTGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTGAGAAAAATCAACAAACTGTCAGTAATGATTTGTTGCCTGCCGTCCTTTg                         <  2:200062/151‑1 (MQ=255)
                                                                                                                            caTGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTGAGAAAAATCAACAAACTGTCAGTAATGATTTGTTGCCTGCCGTCCTTTGTtat                     >  2:161518/1‑151 (MQ=255)
                                                                                                                            caTGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTGAGAAAAATCAACAAACTGTCAGTAATGATTTGTTGCCTGCCGTCCTTTGTtat                     <  1:114234/151‑1 (MQ=255)
                                                                                                                            caTGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTGAGAAAAATCAACAAACTGTCAGTAATGATTTGTTGCCTGCCGTCCTTTGTtat                     <  1:966248/151‑1 (MQ=255)
                                                                                                                                      ccGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTGAGAAAAATCAACAAACTGTCAGTGATGATTTGTTGCCTGCCGTCCTTTGTTATACCGTCTCTg           <  2:1160487/151‑1 (MQ=255)
                                                                                                                                       cGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTGAGAAAAATCAACAAACTGTCAGTAATGATTTGTTGCCTGCCGTCCTTTGTTATACCGTCTCTGc          >  2:954105/1‑151 (MQ=255)
                                                                                                                                        gTCCGATCACACCCCCGGCGGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTGAGAAAAATCAACAAACTGTCAGTAATGATTTGTTGCCTGCCGTCCTTTGTTATACCGTCTCTGc          <  2:344699/150‑1 (MQ=255)
                                                                                                                                          ccGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTGAGAAAAATCAACAAACTGTCAGTAATGATTTGTTGCCTGCCGTCCTTTGTTATACCGTCTCTGCGtt       >  2:86393/1‑151 (MQ=255)
                                                                                                                                          ccGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTGAGAAAAATCAACAAACTGTCAGTAATGATTTGTTGCCTGCCGTCCTTTGTTATACCGTCTCTGCGt        <  2:111471/150‑1 (MQ=255)
                                                                                                                                               cacaCCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTGAGAAAAATCAACAAACTGTCAGTAATGATTTGTTGCCTGCCGTCCTTTGTTATACCGTCTCTGCGTTTTTa   <  2:608111/150‑1 (MQ=255)
                                                                                                                                               cacaCCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTGAGAAAAATCAACAAACTGTCAGTAATGATTTGTTGCCTGCCGTCCTTTGTTATACCGTCTCTGCGTTTTTAg  <  2:881439/151‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CATGGTTTGATTACCGCTCAAAAGGTTTTGACGATAACCGTGGTCTGCGCATCGACCTGCTGCTCGCCAGCCAACCGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACGCCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTGAGAAAAATCAACAAACTGTCAGTAATGATTTGTTGCCTGCCGTCCTTTGTTATACCGTCTCTGCGTTTTTAG  >  NC_000913/1833059‑1833352

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