Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,426,592 T→A 7.0% intergenic (+35/+36) rplI → / ← yjfZ 50S ribosomal subunit protein L9/uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,426,5920TA 7.0% 415.8 / 15.8 153intergenic (+35/+36)rplI/yjfZ50S ribosomal subunit protein L9/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base T (77/65);  new base A (7/4);  total (84/69)
Fisher's exact test for biased strand distribution p-value = 7.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.10e-01

Only 100 of 153 total aligned reads displayed.

GCGAAGTTCGTCTGCCGAACGGCGTTCTGCGTACCACTGGCGAACACGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCAACTGATGTCTCGCTCCTGCGGCACTCGC  >  NC_000913/4426445‑4426739
                                                                                                                                                   |                                                                                                                                                   
gCGAAGTTCGTCTGCCGAACGGCGTTCTGCGTACCACTGGCGAACACGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTg                                                                                                                                                   >  1:1738934/1‑150 (MQ=255)
    aGTTCGTCTGCCGAACGGCGTTCTGCGTACCACTGGCGAACACGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCg                                                                                                                                               >  1:907556/1‑150 (MQ=255)
      ttCGTCTGCCGAACGGCGTTCTGCGTACCACTGGCGATCACGAAGTGAGCTTCCAGGTTCACATCGAAGTATTCGCGAAAGTGATCGTAAACGAAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGc                                                                                                                                             <  1:85817/150‑1 (MQ=255)
      ttCGTCTGCCGAACGGCGTTCTGCGTACCACTGGCGAACACGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCAATGGTCGGc                                                                                                                                             <  2:354779/150‑1 (MQ=255)
       tCGTCTGCCGAACGGCGTTCTGCGTACCACTGGCGAACACGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCg                                                                                                                                            <  2:939786/150‑1 (MQ=255)
        cGTCTGCCGAACGGCGTTCTGCGTACCACTGGCGAACACGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGGCGTAAAAACGCCGACCAATGGTCGGCGt                                                                                                                                           <  2:1804106/150‑1 (MQ=255)
         gtatgacGAACGGCGTTCTGCGTACCACTGGCGAACAAGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGtt                                                                                                                                          <  2:122953/144‑1 (MQ=255)
         gTCTGCCGAACGGCGTTCTGCGTACCACTGGCGAACACGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGtt                                                                                                                                          >  2:104905/1‑150 (MQ=255)
         gTCTGCCGAACGGCGTTCTGCGTACCACTGGCGAACACGAAGGGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGGTATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGtt                                                                                                                                          <  1:1682533/150‑1 (MQ=255)
           cTGCCGAACGGCGTTCTGCGTCCCACTGGCGAACACGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTCGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGCCCATTGGGCGGCGtttt                                                                                                                                        <  1:899938/150‑1 (MQ=255)
            tGCCGAACGGCGTTCTGCGTACCACTGGCGAACACGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTg                                                                                                                                       >  1:815850/1‑150 (MQ=255)
             gCCGAACGGCGTTCTGCGTACCACTGGCGAACACGAAGTGATCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGc                                                                                                                                      <  2:223956/150‑1 (MQ=255)
             gCCGAACGGCGTTCTGCGTACCACTGGCGAACACGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGc                                                                                                                                      >  1:1534287/1‑150 (MQ=255)
               cGAACGGCGTTCTGCGTACCACTGGCGAACACGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCtt                                                                                                                                    >  2:64837/1‑150 (MQ=255)
               cGAACGGCGTTCTGCGTACCACTGGCGAACACGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGGCGGCGTTTTGCtt                                                                                                                                    <  1:505486/150‑1 (MQ=255)
                 aaCGGCGTTCTGCGTACCACTGGCGAACACGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTc                                                                                                                                  <  1:1257246/150‑1 (MQ=255)
                 aaCGGCGTTCTGCGTACCACTGGCGAACACGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTc                                                                                                                                  <  1:1247763/150‑1 (MQ=255)
                             cGTACCACTGGCGAACACGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCa                                                                                                                      >  1:992396/1‑150 (MQ=255)
                              ggaCCACTGGCTAACACGACGTGAGCGTCAAGGTGCACAGCGAAGTATTCGCGAATGTGATCGTAAACGTAGTAGCTGAATTATTCGTTATTCAACGAGACGTAAAAACGCCGACCAATGGTCGGCGTTTTGCTTTCTATTTTTCGTCAg                                                                                                                     <  2:964/148‑1 (MQ=255)
                              gTACCACTGTCGAACACGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAg                                                                                                                     <  1:644925/150‑1 (MQ=255)
                               tACCACTGGCGAACACGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAgg                                                                                                                    <  2:1294132/150‑1 (MQ=255)
                                  cACTGGCGAACACGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGTTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTAt                                                                                                                 >  1:1096959/1‑150 (MQ=255)
                                    cTGGCGAACACGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTa                                                                                                               >  2:1715749/1‑150 (MQ=255)
                                        cGAACACGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGttt                                                                                                           >  2:1191074/1‑150 (MQ=255)
                                        cGAACACGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGAACGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGatt                                                                                                           >  1:107235/1‑150 (MQ=255)
                                          aaCACGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGGGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCg                                                                                                         <  2:890904/150‑1 (MQ=255)
                                          aaCACGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCg                                                                                                         <  2:1544630/150‑1 (MQ=255)
                                            cacGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCa                                                                                                       <  2:1032502/150‑1 (MQ=255)
                                            cacGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCa                                                                                                       >  1:183158/1‑150 (MQ=255)
                                            cacGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCa                                                                                                       >  2:835869/1‑150 (MQ=255)
                                              cGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAg                                                                                                     >  2:1649283/1‑150 (MQ=255)
                                              cGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAg                                                                                                     >  2:27046/1‑150 (MQ=255)
                                               gAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGt                                                                                                    >  1:864894/1‑150 (MQ=255)
                                                aaGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTa                                                                                                   >  2:502067/1‑150 (MQ=255)
                                                   tGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCTGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGAt                                                                                                >  1:296485/1‑150 (MQ=255)
                                                    gAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATc                                                                                               >  2:1642132/1‑150 (MQ=255)
                                                     aGCTTCCGGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATcc                                                                                              >  1:1211492/1‑150 (MQ=255)
                                                         tCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTc                                                                                          <  2:1242/150‑1 (MQ=255)
                                                          ccAGGTTCACAGCGAAGTATTCCCGAAAGTGATCGTAAACGTTGTAGCTGAATATTTAGTTATTCAAGGAGAGGTATAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTATCGCAAGTTGATCCATTCa                                                                                         <  1:309851/150‑1 (MQ=255)
                                                          ccAGGTTCACAGAGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACAATTGGTCGGCGTTATGATTTCTATTTTTCGTCAGGTATTAGATTCGCAAGTAGATCCATTCa                                                                                         >  1:1390975/1‑150 (MQ=255)
                                                           cAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCAATTCAc                                                                                        >  1:1131421/1‑150 (MQ=255)
                                                                   caGCGACGTAGTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACCTAAAAACCCCGACCATAGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCACGTAGATCCATTCACCATTctct                                                                                <  2:299480/150‑1 (MQ=255)
                                                                     gCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTtg                                                                              <  1:1296915/150‑1 (MQ=255)
                                                                     gCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTtg                                                                              <  2:1114449/150‑1 (MQ=255)
                                                                       gAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTtgtg                                                                            >  1:1493114/1‑150 (MQ=255)
                                                                         aGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTc                                                                          <  1:511582/150‑1 (MQ=255)
                                                                          gTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTcc                                                                         >  1:1417576/1‑150 (MQ=255)
                                                                           tATTCGCGAAAGTGATCGTATACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAATACGCCGACCATTGGTCGGCGTTTTGCTTTCTTTTTTTCGTCAGGTATTAGTTTCGCAAGTTGATGCATTCACCATTCTGTTGTGTCCa                                                                        >  2:191113/1‑150 (MQ=255)
                                                                            tttCGCGAAAGTGATCGTATACGTCGTAGCTGAATAATTCGTTATTCAACGAGCAGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCAt                                                                       <  2:308903/149‑1 (MQ=255)
                                                                               cgcgAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGCTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACGATTCTCTTGTGTCCATaat                                                                    >  2:1378465/1‑149 (MQ=255)
                                                                                  gAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCa                                                                 <  2:337000/150‑1 (MQ=255)
                                                                                       tGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCTGTTTGCGTTCGATTTTTCGTCAGGTATTTGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATTAGTCAggtag                                                            >  2:458038/1‑147 (MQ=255)
                                                                                        gATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTTTTAGTTTCGCAATTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATgg                                                           >  2:61036/1‑150 (MQ=255)
                                                                                           cGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGttt                                                        >  2:1839369/1‑150 (MQ=255)
                                                                                           cGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTGTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTTGATCCATTCACAATTCTCTTGTGTCCATAAGTCAGGATGGttt                                                        >  2:346045/1‑150 (MQ=255)
                                                                                            gTAAACGTAGTAGCTGAATAATTAGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTg                                                       >  1:1018335/1‑150 (MQ=255)
                                                                                             tAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGc                                                      <  1:1007530/150‑1 (MQ=255)
                                                                                              aaaCGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTACATAAGTCAGGATGGTTTGCTg                                                    >  2:1147003/1‑150 (MQ=255)
                                                                                                aCGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTgg                                                   >  2:846101/1‑150 (MQ=255)
                                                                                                  gtagtaGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCg                                                 <  2:1132213/150‑1 (MQ=255)
                                                                                                  gtagtaGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCg                                                 <  1:1282819/150‑1 (MQ=255)
                                                                                                     gtagCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCTATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTa                                              <  2:1305596/150‑1 (MQ=255)
                                                                                                     gtagCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTa                                              <  1:1462017/150‑1 (MQ=255)
                                                                                                      tagCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTTGTCGGCGTTTTGCTTTCTATTTTTCGCCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAg                                             <  2:1322968/150‑1 (MQ=255)
                                                                                                      tagCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAg                                             <  1:430512/150‑1 (MQ=255)
                                                                                                      tagCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCAACATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAg                                             >  1:226333/1‑150 (MQ=255)
                                                                                                      tagCTGAATAATTCGTTATTCAACGAGACGTAAAAACACCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAg                                             <  2:1467716/150‑1 (MQ=255)
                                                                                                       agCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCAt                                                                                     >  2:1097275/1‑109 (MQ=255)
                                                                                                        gCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCTTAAGTCAGGATGGTTTGCTGGCGTTAGat                                           >  2:95311/1‑150 (MQ=255)
                                                                                                        gCTGAATAATTCGTTATTCAACCAGATGTAAAAACGACGACCATTGGTCGGCGTTTTTCTTTATATTTTTCGTCAGGTATAAGTTTAGCAAGTAGATACATTCACCATTATCTTATGTCCATAAGTAATGATGGTTTGCTGGAGTTAGat                                           >  1:1381909/1‑150 (MQ=255)
                                                                                                         cTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGata                                          <  2:1107071/150‑1 (MQ=255)
                                                                                                                atTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGa                                   >  1:1670603/1‑150 (MQ=255)
                                                                                                                atTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCAATATGTCAGGATGGTTTGCTGGCGTTAGATACATACgc                                   >  1:1148975/1‑149 (MQ=255)
                                                                                                                   cGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCCTTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACTAgcg                                >  2:109337/1‑150 (MQ=255)
                                                                                                                   cGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAgcg                                >  1:964508/1‑150 (MQ=255)
                                                                                                                     ttATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCa                              <  2:34699/150‑1 (MQ=255)
                                                                                                                       aTTCAACGAGACGTAAAAACTCCGACCATTGGTCGGCGTTTTGGTTTCTATTTTTCGTCAGGAATTAGTTTCGCAAGTAGATCCATTCACCATTATCTTGAGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCAAc                            >  1:1865227/1‑150 (MQ=255)
                                                                                                                         tCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCAACTg                          <  2:1838033/150‑1 (MQ=255)
                                                                                                                          cAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCAACTGa                         >  2:1285731/1‑150 (MQ=255)
                                                                                                                          cAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCAACTGa                         >  1:179503/1‑150 (MQ=255)
                                                                                                                            aagcttcGTAAAAACGCAGACCTTTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCAACTGATg                       <  2:787209/144‑1 (MQ=255)
                                                                                                                                gaCGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCAACTGATGtctc                   <  1:94080/150‑1 (MQ=255)
                                                                                                                                 aCGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCAACTGATGTCTCg                  <  2:666978/150‑1 (MQ=255)
                                                                                                                                  cGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCAACTGATGTCTCGc                 >  1:990383/1‑150 (MQ=255)
                                                                                                                                    gcaaaaCGCCGACCAATGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCAACTGATGTCTCGCTc               >  2:1351530/3‑150 (MQ=255)
                                                                                                                                     caaaaCGCCGACCAATGGTCGGCGTTTTGATTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTTGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGTATGGTTTGCTGGCGTTAGATTCATACGCGCGCAACTGATGTCTCGCTcc              >  2:415680/2‑150 (MQ=255)
                                                                                                                                     aaaaaCGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCAACTGATGTCTCGCTcc              <  1:761268/150‑1 (MQ=255)
                                                                                                                                          cGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATACACAATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCAACTGAAGTCTCGCTCCTGCgg         >  2:802582/1‑150 (MQ=255)
                                                                                                                                           gCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATTCATACGAGCGCCACTGATGTCTCGCTCCTGCggg        >  2:363397/1‑149 (MQ=255)
                                                                                                                                           gCCGACCAATGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCCCCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCAACTGATGTCTCGCTCCTGCGGc        >  2:42616/1‑150 (MQ=255)
                                                                                                                                            ccGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCAACTGATGTCTCGCTCCTGCGGcc       >  2:675623/1‑149 (MQ=255)
                                                                                                                                            ccGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCAACTGATGTCTCGCTCCTGCGGCa       >  2:984572/1‑150 (MQ=255)
                                                                                                                                            ccGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCAACTGATGTCTCGCTCCTGCGGCa       >  1:371280/1‑150 (MQ=255)
                                                                                                                                             cGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCAACTGATGTCTCGCTCCTGCGGCAc      >  1:1417846/1‑150 (MQ=255)
                                                                                                                                             cGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCAACTGATGTCTCGCTCCTGCGGCAc      <  1:27046/150‑1 (MQ=255)
                                                                                                                                               cccATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTAGTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGGTAGATACATACGAGCGCAACTGATGTCTCGCTCCTGCGGCActc    <  2:1554554/149‑1 (MQ=255)
                                                                                                                                               atcATTGGTCGGCGTTTTGATTTCTATTTTTATTCAGGTATTAGTTTCGCAAGTAGATCCATTCAACATTCTCTTGTTTCCATAAGTCAGGATTGTTTGCTGGCGTTAGATACATACTAGCGCAACTGATTTCTCGCTCCTGCGGCActc    <  2:1471955/148‑1 (MQ=255)
                                                                                                                                               aCCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCAACTGATGTCTCGCTCCTGCGGCActc    >  2:485692/1‑150 (MQ=255)
                                                                                                                                               aCCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCAACTGATGTCTCGCTCCTGCGGCActc    <  2:1744286/150‑1 (MQ=255)
                                                                                                                                                 cATTGGTCGGCGTTTTGCTGTCTATTTTTCGTCAGGTATTAGATTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCAACTGAAGTCTCGCTCCTGCGGCActcgc  >  1:1128088/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GCGAAGTTCGTCTGCCGAACGGCGTTCTGCGTACCACTGGCGAACACGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCAACTGATGTCTCGCTCCTGCGGCACTCGC  >  NC_000913/4426445‑4426739

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