Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,172,005 T→C 5.3% intergenic (+249/‑52) rrfB → / → murB 5S ribosomal RNA of rrnB operon/UDP‑N‑acetylenolpyruvoylglucosamine reductase, FAD‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,172,0050TC 5.3% 233.0 / 3.1 75intergenic (+249/‑52)rrfB/murB5S ribosomal RNA of rrnB operon/UDP‑N‑acetylenolpyruvoylglucosamine reductase, FAD‑binding
Reads supporting (aligned to +/- strand):  ref base T (37/34);  new base C (2/2);  total (39/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.15e-01

CGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCCACAGATACGGTAAACTAGCCTCGTTTTTGCATCAGGAAAGCAGCTATGAACCACTCCTTAAAACCCTGGAACACATTTGGCATTGATCATAATGCTCAGCACATTGTATGTGCCGAAGACGAACAACAATTACTCAATGCCT  >  NC_000913/4171864‑4172153
                                                                                                                                             |                                                                                                                                                    
cGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCca                                                                                                                                               <  6:180037/149‑1 (MQ=255)
  ttGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGTCGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCcaca                                                                                                                                             <  4:323921/149‑1 (MQ=255)
      gAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCCACAGATa                                                                                                                                         >  4:27753/1‑149 (MQ=255)
       aaGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCCACAGATAc                                                                                                                                        >  6:211146/1‑149 (MQ=255)
             cGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTACTGTCGTCATATCTACAAGCCATCCCCCCACAGATACGGTaa                                                                                                                                   >  8:12876/1‑148 (MQ=255)
                 ccGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCCACAGATACGGTAAACTAg                                                                                                                              <  1:183822/149‑1 (MQ=255)
                 ccGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCCACAGATACGGTAAACTAg                                                                                                                              >  4:237244/1‑149 (MQ=255)
                     aGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACACGCCATCCCCCCACAGATACGGTAAACTAGCCTc                                                                                                                          >  3:72109/1‑149 (MQ=255)
                                 aGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCCACAGATACGGTAAACTAGCCTCGTTTTTGCATCa                                                                                                              <  5:211146/149‑1 (MQ=255)
                                   gACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCCACAGATACGGTAAACTAGCCTCGTTTTTGCATCAg                                                                                                             >  4:99596/1‑148 (MQ=255)
                                        ccgccATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCCACAGATACGGTAAACTAGCCTCGTTTTTGCATCAGGAAag                                                                                                        >  1:6545/1‑148 (MQ=255)
                                            cATAAACTGCGAGGCATCAAATTAAGCAGAAGGCCATCCTGCCGGATGGGCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATCTCGACACGCCCTCCGCCCGCAGATACGGTGACCGGGCCGCGGGTTTGCATCCGGCAagcagc                                                                                                    >  2:113269/1‑148 (MQ=18)
                                            cATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCCACAGATACGGTAAACTAGCCTCGTTTTTGCATCAGGAAagcagc                                                                                                    >  8:193664/1‑148 (MQ=255)
                                                aaCTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCCACAGATACGGTAAACTAGCCTCGTTTTTGCATCAGGAAAGCAGCTATg                                                                                                >  2:62141/1‑148 (MQ=255)
                                                 aCTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCCACAGATACGGTAAACTAGCCTCGTTTTTGCATCAGGAAAGCAGCTATGaa                                                                                              <  1:62141/149‑1 (MQ=255)
                                                 aCTGCCAGGCATCAAATTAAGCAGAAGGCCATACAGACGGATGGCCTTTTTGCGTTACTAAAAACTATTCCTGTAGTCATAAATACAAACCATCCCCCCAAAGATAAGGTAAACAAGCCACGTTTTTGAATCAGGAAAGCACCAATGaa                                                                                              >  8:216155/1‑149 (MQ=21)
                                                  cTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCCACAGATACGGTAAACTAGCCTCGTTTTTGCATCAGGAAAGCAGCTATGAAc                                                                                             <  2:6545/149‑1 (MQ=255)
                                                     ccAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCCACAGATACGGTAAACTAGCCTCGTTTTTGCATCAGGAAAGCAGCTATGAACCAc                                                                                          <  7:12876/149‑1 (MQ=255)
                                                       aGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCAGTTTTGCGGTTCTACAAACTCTTCCTGTCGTCATATCTCCAAGCCATCCCCCCACAGAGACGGTAAACTAGCCTCGTTTTTGCATAAGGCAAGCACCTATGAACCACt                                                                                         >  3:235566/1‑148 (MQ=38)
                                                        ggCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCCACAGATACGGTAAACTAGCCTCGTTTTTGCATCAGGAAAGCAGCTATGAACCACTcc                                                                                       <  3:237244/149‑1 (MQ=255)
                                                         gCATCAAATTAAGC‑‑AAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCACCCCCCCACAGATACGGTAAACTAGCCTCGTTTTTGCATCAGGaaa                                                                                                         >  5:51305/1‑128 (MQ=255)
                                                         gCATCAAATTAAGC‑‑AAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCACCCCCCCACAGATACGGTAAACTAGCCTCGTTTTTGCATCAGGaaa                                                                                                         <  6:51305/128‑1 (MQ=255)
                                                          cATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCCACAGATACGGTAAACTAGCCTCGTTTTTGCATCAGGAAAGCAGCTATGAACCACTCCtt                                                                                     <  3:27753/149‑1 (MQ=255)
                                                           aTCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTATTACAGTCGTCATATCTACAAGCCATCCCCCAACAGATACGGTAAACTAGCCTCGTTTTTGCATCAGGAAAGCAGCTATGAACCCCTCCtt                                                                                     >  3:114647/1‑148 (MQ=38)
                                                            tCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCCACAGATACGGTAAACTAGCCTCGTTTTTGCATCAGGAAAGCAGCTATGAACCACTCCTTaa                                                                                   <  3:335570/149‑1 (MQ=255)
                                                              aaaTTAATCATAAGGGCATCATTACGTATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCCACAGATACGGTAAACTAGCCTCGTTTTTGCATCAGGAAAGCAGCTATGAACCACTCCTTaaaa                                                                                 <  7:216155/149‑1 (MQ=255)
                                                              aaaTTAAGCAGAAGGCCATCCTGACGGATGGTCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCCACAGATACGGTAAACTAGCCTCGTTTTTGCATCAGGAAAGCAGCTATGAACCACTCCTTaaaa                                                                                 <  1:97330/149‑1 (MQ=255)
                                                                aTTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCCACAGATACGGTAAACTAGCCTCGTTTTTGCATCAGGAAAGCAGCTATGAACCACTCCTTAAAAcc                                                                               >  5:205995/1‑149 (MQ=255)
                                                                    aGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCCACAGATACGGTAAACTAGCCTCGTTTTTGCATCAGGAAAGCAGCTATGAACCACTCCTTAAAACCCTg                                                                            <  4:72109/148‑1 (MQ=255)
                                                                    aGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCCACAGATACGGTAAACTAGCCTCGTTTTTGCATCAGGAAAGCAGCTATGAACCACTCCTTAAAACCCTg                                                                            <  3:99596/148‑1 (MQ=255)
                                                                       aGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCCACAGATACGGTAAACTAGCCTCGTTTTTGCATCAGGAAAGCAGCTATGAACCACTCCTTAAAACCCTGGaa                                                                         >  8:228248/1‑148 (MQ=255)
                                                                       aGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCCACAGATACGGTAAACTAGCCTCGTTTTTGCATCAGGAAAGCAGCTATGAACCACTCCTTAAAACCCTGGaa                                                                         >  5:171141/1‑148 (MQ=255)
                                                                       aGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCACCCCCCCACAGATACGGTAAACTAGCCTCGTTTTTGCATCAGGAAAGCAGCTATGAACCACTCCTTAAAACCCTgg                                                                           <  2:84920/146‑1 (MQ=255)
                                                                       aGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCACCCCCCCACAGATACGGTAAACTAGCCTCGTTTTTGCATCAGGAAAGCAGCTATGAACCACTCCTTAAAACCCTGGaag                                                                        >  1:84920/1‑148 (MQ=255)
                                                                        gAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCCACAGATACGGTAAACTAGCCTCGTTTTTGCATCAGGAAAGCAGCTATGAACCACTCCTTAAAACCCTGGAaca                                                                       >  3:149169/1‑149 (MQ=255)
                                                                                tCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCCACAGATACgg                                                                                                                                      <  2:113150/78‑1 (MQ=255)
                                                                                tCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCCACAGATACgg                                                                                                                                      >  1:113150/1‑78 (MQ=255)
                                                                                tCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCCACAGATACGGTAAACTAGCCTCGTTTTTGCATCAGGAAAGCAGCTATGAACCACTCCTTAAAACCCTGGAACACATTTGGc                                                               <  7:228248/149‑1 (MQ=255)
                                                                                     aCGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCCACAGATACGGTAAACTAGCCTCGTTTTTGCATCAGGAAAGCAGCTATGAACCACTCCTTAAAACCCTGGAACACATTTGGCATTGa                                                          <  6:156497/149‑1 (MQ=255)
                                                                                     aCGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCCACAGATACGGTAAACTAGCCTCGTTTTTGCATCAGGAAAGCAGCTATGAACCACTCCTTAAAACCCTGGAACACATTTGGCATTGa                                                          <  7:193664/149‑1 (MQ=255)
                                                                                         aTGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCCACAGATACGGTAAACTAGCCTCGTTTTTGCATCAGGAAAGCAGCTATGAACCACTCCTTAAAACCCTGGAACACATTTGGCATTGATCa                                                       >  4:142405/1‑148 (MQ=255)
                                                                                            gCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCCACAGATACGGTAAACTAGCCTCGTTTTTGCAGCAGGAAAGCAGCTATGAACCACTCCTTAAAACCCTGGAACACATTTGGCATTGATCATaa                                                    <  1:113269/148‑1 (MQ=255)
                                                                                            gCCGTTTTGTGTTGCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCCACAGATACGGTAAACTAGCCTCGTTTTTGCATCAGGAAAGCAGCTATGAACCACTCCTTAAAACCCTGGAACACATTTGGCATTGATCATaa                                                    <  4:15761/148‑1 (MQ=255)
                                                                                              cgttttGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCCACAGATACGGTAAACTAGCCTCGTTTTTGCATCAGGAAAGCAGCTATGAACCACTCCTTAAAACCCTGGAACACATTTGGCATTGATCATAATGc                                                 >  6:48775/3‑149 (MQ=255)
                                                                                              cTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCCACAGATACGGTAAACTAGCCTCGTTTTTGCATCAGGAAAGCAGCTATGAACCACTCCTTAAAACCCTGGAACACATTTGGCATTGATCATAATGc                                                 <  6:171141/149‑1 (MQ=255)
                                                                                                          cTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCCACAGATACGGTAAACTAGCCTCGTTTTTGCATCAGGAAAGCAGCTATGAACCACTCCTTAAAACCCTGGAACACATTTGGCATTGATCATAATGCTCAGCACATTg                                      >  4:98827/1‑148 (MQ=255)
                                                                                                           tACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCCACAGATACGGTAAACTAGCCTCGTTTTTGCATCAGGAAAGCAGCTATGAACCACTCCTTAAAACCCTGGAACACATTTGGCATTGATCATAATGCTCAGCACAtt                                       <  4:347519/146‑1 (MQ=255)
                                                                                                           tACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCCACAGATACGGTAAACTAGCCTCGTTTTTGCATCAGGAAAGCAGCTATGAACCACTCCTTAAAACCCTGGAACACATTTGGCATTGATCATAATGCTCAGCACAtt                                       >  3:347519/1‑146 (MQ=255)
                                                                                                              aaaCTCTTCCTGTCGTCATATCTACAAGCCATCCCCCCACAGATACGGTAAACTAGCCTCGTTTTTGCATCAGGAAAGCAGCTATGAACCACTCCt                                                                                      >  4:289546/1‑96 (MQ=255)
                                                                                                              aaaCTCTTCCTGTCGTCATATCTACAAGCCATCCCCCCACAGATACGGTAAACTAGCCTCGTTTTTGCATCAGGAAAGCAGCTATGAACCACTCCt                                                                                      <  3:289546/96‑1 (MQ=255)
                                                                                                              aaaCTCTTCCTGTCGTCATATCTACAAGCCATCCCCCCACAGATACGGTAAACTAGCCTCGTTTTTGCATCAGGAAAGCAGCTATGAACCACTCCTTAAAACCCTGGAACACATTTg                                                                 <  8:153287/117‑1 (MQ=255)
                                                                                                              aaaCTCTTCCTGTCGTCATATCTACAAGCCATCCCCCCACAGATACGGTAAACTAGCCTCGTTTTTGCATCAGGAAAGCAGCTATGAACCACTCCTTAAAACCCTGGAACACATTTg                                                                 >  7:153287/1‑117 (MQ=255)
                                                                                                              aaaCTCTTCCTGTCGTCATATCTACAAGCCATCCCCCCACAGATACGGTAAACTAGCCTCGTTTTTGCATCAGGAAAGCAGCTATGAACCACTCCTTAAAACCCTGGAACACATTTGGCATTGATCATAATGCTCAGCACATTGTAtgt                                 >  4:351525/1‑149 (MQ=255)
                                                                                                                aCTCTTCCTGTCGTCATATCTACAAGCCATCCCCCCACAGATACGGTAAACTAGCCTCGTTTTTGCATCAGGAAAGCAGCTATGAACCACTCCTTAAAACCCTGGAACACATTTGGCATTGATCATAATGc                                                 <  8:299175/131‑1 (MQ=255)
                                                                                                                aCTCTTCCTGTCGTCATATCTACAAGCCATCCCCCCACAGATACGGTAAACTAGCCTCGTTTTTGCATCAGGAAAGCAGCTATGAACCACTCCTTAAAACCCTGGAACACATTTGGCATTGATCATAATGc                                                 >  7:299175/1‑131 (MQ=255)
                                                                                                                      ccTGTCGTCATATCTACAAGCCATCCCCCCACAGATACGGTAAACTAGCCTCGTTTTTGCATCAGGAAAGCAGCTATGAACCACTCCTTAAAACCCTGGAACACATTTGGCATTGATCATAATGCTCAGCACATTGTATGTGCCGAAg                          <  6:205995/148‑1 (MQ=255)
                                                                                                                          tcgtcATATCTACAAGCCATCCCCCCACAGATACGGTAAACTAGCCTCGTTTTTGCATCAGGAAAGCAGCTATGAACCACTCCTTAAAACCCTGGAACACATTTGGCATTGATCATAATGCTCAGCACATTGTATGTGCCGAAGACGa                      <  3:98827/148‑1 (MQ=255)
                                                                                                                            gtcATATCTACAAGCCATCCCCCCACAGATACGGTAAACTAGCCTCGTTTTTGCATCAGGAAAGCAGCTATGAACCACTCCTTAAAACCCTGGAACACATTTGGCATTGATCATAATGCTCAGCACATTGTATGTGCCGAAGACGaac                    <  3:142405/148‑1 (MQ=255)
                                                                                                                             tcATATCTACAAGCCATCCCCCCACAGATACGGTAAACTAGCCTCGTTTTTGCATCAGGAAAGCAGCTATGAACCACTCCTTAAAACCCTGGAACACATTTGGCATTGATCATAATGCTCAGCACATTGTATGTGCCGAAGACGaacaa                  <  4:149169/149‑1 (MQ=255)
                                                                                                                             tcATATCTACAAGCCATCCCCCCACAGATACGGTAAACTAGCCTCGTTTTTGCATCAGGAAAGCAGCTATGAACCACTCCTTAAAACCCTGGAACACATTTGGCATTGATCATAATGCTCAGCACATTGTATGTGCCGAAGACGaacaa                  >  6:40488/1‑149 (MQ=255)
                                                                                                                              cATATCTACAAGCCATCCCCCCACAGATACGGTAAACTAGCCTCGTTTTTGCATCAGGAAAGCAGCTATGAACCACTCCTTAAAACCCTGGACCACATTTGGCATTGATCATAATGCTCAGCACATTg                                      <  1:333481/128‑1 (MQ=255)
                                                                                                                              cATATCTACAAGCCATCCCCCCACAGATACGGTAAACTAGCCTCGTTTTTGCATCAGGAAAGCAGCTATGAACCACTCCTTAAAACCCTGGAACACATTTGGCATTGATCATAATGCTCAGCACATTg                                      >  2:333481/1‑128 (MQ=255)
                                                                                                                                tatCTACAAGCCATCCCCCCACAGATACGGTAAACTAGCCTCGTTTTTGCATCAGGAAAGCAGCTATGAACCACTCCTTAAAACCCTGGAACAGATTTGGCATTGATCATAATGCTCCGCCCATTGTATGTGCCGAAGCCGACcgacaa               >  5:322688/1‑149 (MQ=255)
                                                                                                                                    tACAAGCCATCCCCCCACAGATACGGTAAACTAGCCTCGTTTTTGCATCAGGAAAGCAGCTATGAACCACTCCTTAAAACCCTGGAACACATTTGGCATTGATCATAATGCTCAGCACAtt                                       <  4:34864/121‑1 (MQ=255)
                                                                                                                                    tACAAGCCATCCCCCCACAGATACGGTAAACTAGCCTCGTTTTTGCATCAGGAAAGCAGCTATGAACCACTCCTTAAAACCCTGGAACACATTTGGCATTGATCATAATGCTCAGCACAtt                                       >  3:34864/1‑121 (MQ=255)
                                                                                                                                    tACAAGCCATCCCCCCACAGATACGGTAAACTAGCCTCGTTTTTGCATCAGGAAAGCAGCTATGAACCACTCCTTAAAACCCTGGAACACATTTGGCATTGATCATAATGCTCAGCACATTGTATGTGCCGAAGACGAACAACAATTAc           >  4:252271/1‑149 (MQ=255)
                                                                                                                                     aCAAGCCATCCCCCCACAGATACGGTAAACTAGCCTCGTTTTTGCATCAGGAAAGCAGCTATGAACCACTCCTTAAAACCCTGGAACACATTTGGCATTGATCATAATGCTCAGCACATTGTATGTGCCGAAGACGAACAACAATTAc           >  2:363074/1‑148 (MQ=255)
                                                                                                                                     aCAAGCCATCCCCCCACAGATACGGTAAACTAGCCTCGTTTTTGCATCAGGAAAGCAGCTATGAACCACTCCTTAAAACCCTGGAACACATTTGGCATTGATCATAATGCTCAGCACATTGTATGTGCCGAAGACGAACAACAATTAc           <  1:363074/148‑1 (MQ=255)
                                                                                                                                         gCCATCCCCCCACAGATACGGTAACCTAGCCGCGTTTTTGCATCAGGAAAGCAGCTATGACCCACTCCGTAAAACCCTGGAACACATTTGGCAGTGATCATAATGCTCAGCACATTGTATGTGCCGAAGACGAACAACATTTACTCaag      >  6:332328/1‑148 (MQ=255)
                                                                                                                                         gCCATCCCCCCACAGATACGGTAAACTAGCCTCGTTTTTGCATCAGGAAAGCAGCTATGAACCACTCCTTAAAACCCTGGAACACATTTGGCATTGATCATAATGCTCAGCACATTGTATGTGCCGAAGa                         <  2:177262/130‑1 (MQ=255)
                                                                                                                                         gCCATCCCCCCACAGATACGGTAAACTAGCCTCGTTTTTGCATCAGGAAAGCAGCTATGAACCACTCCTTAAAACCCTGGAACACATTTGGCATTGATCATAATGCTCAGCACATTGTATGTGCCGAAGa                         >  1:177262/1‑130 (MQ=255)
                                                                                                                                         gCCATCCCCCCACAGATACGGTAAACTAGCCTCGTTTTTGCATCAGGAAAGCAGCTATGAACCACTCCTTAAAACCCTGGAACACATTTGGCATTGATCATAATGCTCAGCACATTGTATGTGCCGAAGACGAACAACAATTACTCAAt      <  3:252271/149‑1 (MQ=255)
                                                                                                                                         gCCATCCCCCCACAGATACGGTAAACTAGCCTCGTTTTTGCATCAGGAAAGCAGCTATGAACCACTCCTTAAAACCCTGGAACACATTTGGCATTGATCATAATGCTCAGCACATTGTATGTGCCGAAGACGAACAACAATTACTCAAt      >  7:49998/1‑149 (MQ=255)
                                                                                                                                          ccATCCCCCCACAGATACGGTAAACTAGCCTCGTTTTTGCATCAGGAAAGCAGCTATGAACCACTCCTTAAAACCCTGGAACACATTTGGCATTGATCATAATGCTCAGCACATTGTATGTGCCGAAGACGAACAACAATTACTCAATg     >  7:146605/1‑149 (MQ=255)
                                                                                                                                            aTCCCCCCACAGATACGGTAAACTAGCCTCGTTTTTGCATCAGGAAAGCAGCTATGAACCACTCCTTAGAACCCTGGAACACATTTGGCATTGATCATAATGCTCAGCACATTGTATGTGCCGAAGACGAACAACAATTACTCAATGcc   >  6:136592/1‑149 (MQ=255)
                                                                                                                                             tCCCCCCACAGATACGGTAAACTAGCCTCGTTTTTGCATCAGGAAAGCAGCTATGAACCACTCCTTAAAACCCTGGAACACATTTGGCATTGATCATAATGCTCAGCACATTGTATGTGCCGAAGACGAACAACAATTACTCAATGCCt  >  7:218409/1‑149 (MQ=255)
                                                                                                                                             tCCCCCCACAGATACGGTAAACTAGCCTCGTTTTTGAATCAGGAAAGCAGCTATGAACCACTCCTTAAAACCCTGGAACACATTTGGCATTGATCATAATGCTCAGCACATTGTATGTGCCGAAGACGAACAACAATTACTCAATGCCt  >  6:106292/1‑149 (MQ=255)
                                                                                                                                             |                                                                                                                                                    
CGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCCACAGATACGGTAAACTAGCCTCGTTTTTGCATCAGGAAAGCAGCTATGAACCACTCCTTAAAACCCTGGAACACATTTGGCATTGATCATAATGCTCAGCACATTGTATGTGCCGAAGACGAACAACAATTACTCAATGCCT  >  NC_000913/4171864‑4172153

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: