Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,213,131 C→T 5.5% noncoding (92/120 nt) rrfE → 5S ribosomal RNA of rrnE operon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,213,1310CT 5.5% 791.0 / 21.6 236noncoding (92/120 nt)rrfE5S ribosomal RNA of rrnE operon
Reads supporting (aligned to +/- strand):  ref base C (100/123);  new base T (11/2);  total (111/125)
Fisher's exact test for biased strand distribution p-value = 7.80e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 248 total aligned reads displayed.

CAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCG  >  NC_000913/4213014‑4213280
                                                                                                                     |                                                                                                                                                      
cAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAa                                                                                                                        >  1:3673050‑M2/1‑149 (MQ=255)
cAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAa                                                                                                                        <  1:743215‑M2/150‑2 (MQ=255)
cAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAa                                                                                                                        >  2:1513327‑M2/1‑149 (MQ=255)
 ttctgcgttCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAat                                                                                                                       <  2:954667‑M2/146‑3 (MQ=255)
  gAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTGGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAata                                                                                                                      <  1:2488449‑M2/150‑4 (MQ=255)
    aGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAg                                                                                                                    <  1:14675/150‑1 (MQ=32)
     gCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAataaaa                                                                                                                   >  1:2424837‑M2/1‑144 (MQ=255)
     gCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAataaaa                                                                                                                   >  2:3595045‑M2/1‑144 (MQ=255)
     gCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGa                                                                                                                   <  1:1290752/150‑1 (MQ=32)
     gCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGa                                                                                                                   <  1:2264303/150‑1 (MQ=32)
     gCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGa                                                                                                                   <  2:3840877/150‑1 (MQ=32)
      cGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGaa                                                                                                                  >  2:2698292/1‑150 (MQ=33)
      cGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGaa                                                                                                                  <  1:3343150/150‑1 (MQ=33)
         tCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAataaaacgaa                                                                                                               >  1:2568939‑M2/1‑140 (MQ=255)
         tCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGaaaaa                                                                                                               >  2:2184025/1‑150 (MQ=35)
           tGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAataaaacgaaag                                                                                                             <  2:1238412‑M2/150‑13 (MQ=255)
           tGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAAcc                                                                                                             <  1:3036302/150‑1 (MQ=38)
           tGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAAcc                                                                                                             <  1:123257/150‑1 (MQ=38)
             aTAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGGGTGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAAcccc                                                                                                           >  2:217110/1‑150 (MQ=255)
             aTAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAataaaacgaaaggc                                                                                                           <  1:1911743‑M2/150‑15 (MQ=255)
             aTAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAAcccc                                                                                                           <  2:4198043/150‑1 (MQ=39)
              tAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCg                                                                                                          >  1:2129596/1‑150 (MQ=40)
                aaaCAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAataaaacgaaaggccca                                                                                                        >  2:655559‑M2/1‑133 (MQ=255)
                aaaCAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACGCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGt                                                                                                        >  2:782423/1‑150 (MQ=255)
                      aaTTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATaa                                                                                                  <  2:2010108/150‑1 (MQ=255)
                            ccTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCgg                                                                                           <  1:954378/150‑1 (MQ=255)
                            ccTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCgg                                                                                           <  1:4133796/150‑1 (MQ=255)
                            ccTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCgg                                                                                           <  1:3419157/150‑1 (MQ=255)
                             cTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGTAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCggg                                                                                          >  1:4194755/1‑150 (MQ=255)
                             cTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCggg                                                                                          <  2:85445/150‑1 (MQ=255)
                             cTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCggg                                                                                          >  1:1866718/1‑150 (MQ=255)
                              tGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCgggg                                                                                         >  1:2083951/1‑150 (MQ=255)
                              tGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCgggg                                                                                         <  1:196124/150‑1 (MQ=255)
                                gcggcAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGATATGCCAGGCATCAAataaaacgaaaggcccagtcttccgactgagcc                                                                                        >  1:2762330‑M2/1‑117 (MQ=255)
                                gcggcAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGtt                                                                                       >  1:2088793/1‑150 (MQ=255)
                                gcggcAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGtt                                                                                       <  2:1384350/150‑1 (MQ=255)
                                gcggcAGGAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGCTGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACGGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGtt                                                                                       <  1:811943/150‑1 (MQ=255)
                                 cggcAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAataaaacgaaaggcccagtcttccgactgagcct                                                                                       >  1:3604765‑M2/1‑116 (MQ=255)
                                 cggcAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACCGCCAGGCATCAAATTAGAAAAACCCCGG‑CCTTATGGACCGGGGttt                                                                                      <  2:2261721/150‑1 (MQ=14)
                                  ggcAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAGataaaacgaaaggcccagtctttcgactgagcctt                                                                                      >  1:3180592‑M2/1‑115 (MQ=255)
                                    cAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAGataacacgaaaggcccagtcttccgactgagcctttc                                                                                    >  2:2482755‑M2/1‑113 (MQ=255)
                                    cAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAGataaaacgaaaggcccagtcttccgactgagcctttc                                                                                    >  2:2550791‑M2/1‑113 (MQ=255)
                                    cAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGtttttt                                                                                   <  1:3987893/150‑1 (MQ=255)
                                    cAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGtttttt                                                                                   <  2:665492/150‑1 (MQ=255)
                                          cgcgGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCatat                                                                             >  1:1422920/1‑150 (MQ=255)
                                           gcgGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATc                                                                            <  1:1090875/150‑1 (MQ=255)
                                               tggtCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAatt                                                                        >  2:1605041/1‑150 (MQ=255)
                                                  tCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAattatt                                                                     <  1:693002/150‑1 (MQ=255)
                                                  tCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAattatt                                                                     <  2:1303658/150‑1 (MQ=255)
                                                    ccACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTg                                                                   <  2:408017/150‑1 (MQ=255)
                                                    ccACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTg                                                                   <  2:3453135/150‑1 (MQ=255)
                                                    ccACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTg                                                                   >  1:3000722/1‑150 (MQ=255)
                                                     cACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAataaaacgaaaggctcagtcttccgactgagcctttcgttttatttgatgcctg                                                                   >  2:3227805‑M2/1‑96 (MQ=255)
                                                     cACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGAGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATACGG‑CCGGGGTTTTTTGCATATCAATTATTTGc                                                                  <  2:3164129/150‑1 (MQ=255)
                                                      aCCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCa                                                                 >  1:329332/1‑150 (MQ=255)
                                                        cTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCatg                                                               >  2:1767158/1‑150 (MQ=255)
                                                         tGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATTGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCatga                                                              <  1:2042420/150‑1 (MQ=255)
                                                          gACCCCATGCAGATCTCAGAAGTGAAACTCCGTAGGGCCGATGGTAGTGTTGGGTCTCCCCATGCGTGAGTAGGTTACTGCCAGGCATCTAATTAGAAAAACCACGGTCCATAAGG‑CCGGGGTTTTTTGCATATCCATTATTTGCatgat                                                             <  2:953201/150‑1 (MQ=255)
                                                             cccATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGaa                                                          >  1:1478972/1‑150 (MQ=255)
                                                             cccATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGaa                                                          >  1:2930780/1‑150 (MQ=255)
                                                              ccATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAg                                                         <  1:856467/150‑1 (MQ=255)
                                                                  gCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGaa                                                     <  1:988487/150‑1 (MQ=255)
                                                                   ccGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAAt                                                    >  1:1702291/1‑150 (MQ=255)
                                                                        cTCAGAAGTGAACCGCCGTAGCGCCGATGGTAGTGTGGGGTCGCCCCATGCGAGAGTAGGGTCCTGCCAGGCATCAAATTAGAAAAACCCCGGGCCAGCCGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCAt                                               <  2:1910333/150‑1 (MQ=25)
                                                                        cTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCAt                                               <  2:3123199/150‑1 (MQ=255)
                                                                        cTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCAt                                               <  2:3111127/150‑1 (MQ=255)
                                                                         tCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCAtt                                              >  1:3894790/1‑149 (MQ=255)
                                                                         tCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAAGTATTTGCATGATGAAGGGAATCTCATg                                              >  2:169310/1‑150 (MQ=255)
                                                                          cAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGt                                             >  2:1085677/1‑150 (MQ=255)
                                                                           aGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTc                                            >  1:214526/1‑150 (MQ=255)
                                                                           aGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCAATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTc                                            >  1:1192737/1‑150 (MQ=255)
                                                                             aaGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAg                                          >  1:1034902/1‑150 (MQ=255)
                                                                             aaGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAg                                          <  1:505629/150‑1 (MQ=255)
                                                                             aaGTGAAACGCCGTAGCGCCGATGGTAGAGTGGGGTCTCCCCATGCGAGAGTAGGGAACAGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAg                                          <  2:1489462/150‑1 (MQ=255)
                                                                                tGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTGTTTTGCATATCAATTATTTGCATGATGAAGGGCATCTCATGTCAGTTc                                       >  2:1213899/1‑150 (MQ=255)
                                                                                 gAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCt                                      <  2:2348167/150‑1 (MQ=255)
                                                                                  aaaCGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTg                                     <  2:2924096/150‑1 (MQ=255)
                                                                                     cGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGtat                                  <  2:1800053/150‑1 (MQ=255)
                                                                                         gTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGGTTTTTGCATATCAATTATTTGCATGATGAAGGGACTCTCATGTCAGTTCTGTATATcc                              >  2:3762142/1‑150 (MQ=255)
                                                                                           aGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCaa                            <  2:4250583/150‑1 (MQ=255)
                                                                                               ccGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATtc                        <  2:381244/150‑1 (MQ=255)
                                                                                                  aTGTTAGTGTGGGGTCTCCCCTTGCGTGTGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCTGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATTATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATtcgtc                     <  1:1435658/150‑1 (MQ=38)
                                                                                                  aTGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATtcgtc                     <  2:1405698/150‑1 (MQ=255)
                                                                                                  aTGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATtcgtc                     <  1:160235/150‑1 (MQ=255)
                                                                                                  aTGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATtcgtc                     <  2:2612713/150‑1 (MQ=255)
                                                                                                  aTGGGAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATCAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATtcgtc                     <  1:1621089/150‑1 (MQ=255)
                                                                                                   tGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATtcgtcg                    <  2:2136777/150‑1 (MQ=255)
                                                                                                      tAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTaa                 <  1:2788744/150‑1 (MQ=255)
                                                                                                      tAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTaa                 <  1:3917272/150‑1 (MQ=255)
                                                                                                      tAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTaa                 <  2:1170541/150‑1 (MQ=255)
                                                                                                        gtgtGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATc               <  1:2698292/150‑1 (MQ=255)
                                                                                                         tgtgGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCa              >  1:3567641/1‑150 (MQ=255)
                                                                                                          gtgGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGGATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCaa             <  1:2402233/150‑1 (MQ=255)
                                                                                                               gTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTAc        >  2:4056265/1‑150 (MQ=255)
                                                                                                               gTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTAc        >  2:2125200/1‑150 (MQ=255)
                                                                                                               gTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTAc        >  1:4023657/1‑150 (MQ=255)
                                                                                                                   ccccATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGtt    >  1:2976880/1‑150 (MQ=255)
                                                                                                                     ccATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTcg  >  2:332350/1‑150 (MQ=255)
                                                                                                                     ccATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTcg  <  1:685170/150‑1 (MQ=255)
                                                                                                                     ccATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTcg  <  1:4219787/150‑1 (MQ=255)
                                                                                                                     |                                                                                                                                                      
CAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCG  >  NC_000913/4213014‑4213280

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