Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,105,693 C→A 7.4% T58N (ACC→AAC)  hisD → histidinol dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,105,6930CA 7.4% 206.2 / 5.7 68T58N (ACC→AAC) hisDhistidinol dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (31/32);  new base A (3/2);  total (34/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.50e-01

GTCTCAACGCCGCTGATCCGGATTTCGCCCGACATCTGGATCATCTGCTGAGCTGGGAAAGTGTGTCCGATGACGCGGTCAATCAGCGCGTGCTCGACATCATCAAGGCCGTGCGTGAACGCGGTGATGCGGCGCTGGTGGAATTCACCCAGCGTTTCGATGGTGTTGACGCCAAGTCGATTGACGACCTGATCCTCGATCGCGCGCGCCTGGAACTGGCGCTGACCCGTATTACCCCGGTCCAGCGCGAGGCCCTGGAAAAGGCCGCCAACCGTGTGCGCATTTATCACGAGCGG  >  NC_002947/1105546‑1105841
                                                                                                                                                     |                                                                                                                                                     
gTCTCAACGCCGCTGATCCGGATTTCGCCCGACATCTGGATCATCTGCTGAGCTGGGAAAGTGTGTCCGATGACGCGGTCAATCAGCGCGTGCTCGACATCATCAAGGCCGTGCGTGAACGCGGTGATGCGGCGCTGGTGGAATTCAcc                                                                                                                                                      <  4:31664/149‑1 (MQ=255)
   tcAACGCCGCTGATCCGGATTTCGCCCGACATCTGGATCATCTGCTGAGCTGGGAAAGTGTGTCCGATGACGCGGTCAATCAGCGCGTGCTCGACATCATCAAGGCCGTGCGTGAACGCGGTGATGCGGCGCTGGTGGAATTCACCCAg                                                                                                                                                   >  5:217978/1‑149 (MQ=255)
     aaCGCCGCTGATCCGGATTTCGCCCGACATCTGGATCATCTGCTGAGCTGGGAAAGTGTGTCCGATGACGCGGTCAATCAGCGCGTGCTCGACATCATCAAGGCCGTGCGTGAACGCGGTGATGCGGCGCTGGTGGAATTCACCCAGCg                                                                                                                                                 >  1:166108/1‑149 (MQ=255)
     aaCGCCGCTGATCCGGATTTCGCCCGACATCTGGATCATCTGCTGAGCTGGGAAAGTGTGTCCGATGACGCGGTCAATCAGCGCGTGCTCGACATCATCAAGGCCGTGCGTGAACGCGGTGATGCGGCGCTGGGTGAATTCCACCAGCg                                                                                                                                                 >  2:115108/1‑149 (MQ=255)
       cgccgcTGATCCGGATTTCGCCCGACATCTGGATCATCTGCTGAGCTGGGAAAGTGTGTCCGATGACGCGGTCAATCAGCGCGTGCTCGACATCATCAAGGCCGTGCGTGAACGCGGTGATGCGGCGCTGGGTGAATTCCACCAGCgcc                                                                                                                                               <  1:115108/149‑3 (MQ=255)
        gccgcTGATCCGGATTTCGCCCGACATCTGGATCATCTGCTGAGCTGGGAAAGTGTGTCCGATGACGCGGTCAATCAGCGCGTGCTCGACATCATCAAGGCCGTGCGTGAACGCGGTGATGCGGCGCTGGTGGAATTCACCCAGCGttt                                                                                                                                              >  6:109193/1‑149 (MQ=255)
        gccgcTGATCCGGATTTCGCCCGACATCTGGATCATCTGCTGAGCTGGGAAAGTGTGTCCGATGACGCGGTCAATCAGCGCGTGCTCGACATCATCAAGGCCGTGCGTGAACGCGGTGATGCGGCGCTGGTGGAATTCACCCAGCGttt                                                                                                                                              >  3:195428/1‑149 (MQ=255)
            cTGATCCGGAGTTCGCCCGCCCTCTCGCTCCTCTGCTGCGCTGGGAACGTGTGTCCGTTGACGCGGTCAATCAGCGCGTGCTCGACATCATCAAGGCCGTGCGTGAACGCGGTGATGCGGCGCTGGTGGAATTCACCCAGCGTTTCGAt                                                                                                                                          <  2:166108/149‑1 (MQ=255)
                tCCGGATTTCGCCCGACATCTGGATCATCTGCTGAGCTGGGAAAGTGTGTCCGATGACGCGGTCAATCAGCGCGTGCTCGACATCATCAAGGCCGTGCGTGAACGCGGTGATGCGGCGCTGGTGGAATTCACCCAGCGTTTCGATGgtg                                                                                                                                      <  1:172345/149‑1 (MQ=255)
                  cGGATTTCGCCCGACATCTGGATCATCTGCTGAGCTGGGAAAGTGTGTCCGATGACGCGGTCAATCAGCGCGTGCTCGACATCATCAAGGCCGTGCGTGAACGCGGTGATGCGGCGCTGGTGGAATTCACCCAGCGTTTCGATg                                                                                                                                         <  8:267919/144‑1 (MQ=255)
                  cGGATTTCGCCCGACATCTGGATCATCTGCTGAGCTGGGAAAGTGTGTCCGATGACGCGGTCAATCAGCGCGTGCTCGACATCATCAAGGCCGTGCGTGAACGCGGTGATGCGGCGCTGGTGGAATTCACCCAGCGTTTCGATg                                                                                                                                         >  7:267919/1‑144 (MQ=255)
                  cGGATTTCGCCCGACATCTGGATCATCTGCTGAGCTGGGAAAGTGTGTCCGATGACGCGGTCAATCAGCGCGTGCTCGACATCATCAAGGCCGTGCGTGAACGCGGTGATGCGGCGCTGGTGGAATTCACCCAGCGTTTCGATGgtgt                                                                                                                                     <  3:78888/148‑1 (MQ=255)
                     aTTTCGCCCGACATCTGGATCATCTGCTGAGCTGGGAAAGTGTGTCCGATGACGCGGTCAATCAGCGCGTGCTCGACATCATCAAGGCCGTGCGTGAACGCGGTGATGCGGCGCTGGTGGAATTCACCCAGCGTTTCGATGGTGTTGAc                                                                                                                                 >  7:16818/1‑149 (MQ=255)
                         cGCCCGACATCTGGATCATCTGCTGAGCTGGGAAAGTGTGTCCGATGACGCGGTCAATCAGCGCGTGCTCGACATCATCAAGGCCGTGCGTGAACGCGGTGATGCGGCGCTGGTGGAATTCACCCAGCGTTTCGATGGTGTTGACGcc                                                                                                                              <  4:195428/148‑1 (MQ=255)
                            ccGACATCTGGATCATCTGCTGAGCTGGGAAAGTGTGTCCGATGACGCGGTCAATCAGCGCGTGCTCGACATCATCAAGGCCGTGCGTGAACGCGGTGATGCGGCGCTGGTGGAATTCACCCAGCGTTTCGATGGTGTTGACGCCAAGt                                                                                                                          <  8:181447/149‑1 (MQ=255)
                            ccGACATCTGGATCATCTGCTGAGCTGGGAAAGTGTGTCCGATGACGCGGTCAATCAGCGCGTGCTCGACATCATCAAGGCCGTGCGTGAACGCGGTGATGCGGCGCTGGTGGAATTCACCCAGCGTTTCGATGGTGTTGACGCCAAGt                                                                                                                          >  3:46782/1‑149 (MQ=255)
                            ccGACATCTGGATCATCTGCTGAGCTGGGAAAGTGTGTCCGATGACGCGGTCAATCAGCGCGTGCTCGACATCATCAAGGCCGTGCGTGAACGCGGTGATGCGGCGCTGGTGGAATTCACCCAGCGTTTCGATGGTGTTGACGCCAAGt                                                                                                                          >  4:157287/1‑149 (MQ=255)
                            ccGACATCTGGATCATCTGCTGAGCTGGGAAAGTGTGTCCGATGACGCGGTCAATCAGCGCGTGCTCGACATCATCAAGGCCGTGCGTGAACGCGGTGATGCGGCGCTGGGTGAATTCCACCAGCgccgc                                                                                                                                             <  5:156052/130‑5 (MQ=255)
                            ccGACATCTGGATCATCTGCTGAGCTGGGAAAGTGTGTCCGATGACGCGGTCAATCAGCGCGTGCTCGACATCATCAAGGCCGTGCGTGAACGCGGTGATGCGGCGCTGGGTGAATTCCACCAGCgccgc                                                                                                                                             >  6:156052/1‑126 (MQ=255)
                             cGACATCTGGATCATCTGCTGAGCTGGGAAAGTGTGTCCGATGACGCGGTCAATCAGCGCGTGCTCGACATCATCAAGGCCGTGCGTGAACGCGGTGATGCGGCGCTGGTGGAATTCACCCAGCGTTTCGATGGTGTTGACGCCAAGTc                                                                                                                         >  7:127919/1‑149 (MQ=255)
                                 aTCTGGATCATCTGCTGAGCTGGGAAAGTGTGTCCGATGACGCGGTCAATCAGCGCGTGCTCGACATCATCAAGGCCGTGCGTGAACGCGGTGATGCGGCGCTGGTGGAATTCACCCAGCGTTTCGATGGTGTTGACGCCAAGTCGAtt                                                                                                                     <  8:127919/149‑1 (MQ=255)
                                 aTCTGGATCATCTGCTGAGCTGGGAAAGTGTGTCCGATGACGCGGTCAATCAGCGCGTGCTCGACATCATCAAGGCCGTGCGTGAACGCGGTGATGCGGCGCTGGTGGAATTCACCCAGCGTTTCGATGGTGTTGACGCCAAGTCGAtt                                                                                                                     <  4:46782/149‑1 (MQ=255)
                                 aTCAGGATCATCTGCTGAGCTGGGAAAGTGTGTCCGATGACGCGGTCAATCAGCGCGTGCTCGACATCATCAAGGCCGTGCGTGAACGCGGTGATGCTGCGCTGGTGGAATTCACCCAGCGTTTCGATGGTGTTTACGCCAAGTCGAtt                                                                                                                     <  4:277416/149‑1 (MQ=255)
                                      gATCATCTGCTGAGCTGGGAAAGTGTGTCCGATGACGCGGTCAATCAGCGCGTGCTCGACATCATCAAGGCCGTGCGTGAACGCGGTGATGCGGCGCTGGTGGAATTCACCCAGCGTTTCGATGGTGTTGACGCCAAGTCGATTgac                                                                                                                  <  7:142676/147‑1 (MQ=255)
                                       atcatcTGCTGAGCTGGGAAAGTGTGTCCGATGACGCGGTCAATCAGCGCGTGCTCGACATCATCAAGGCCGTGCGTGAACGCGGTGATGCGGCGCTGGTGGAATTCACCCAGCGTTTCGATGGTGTTGACGCCAAGTCGATTgacgac                                                                                                               <  4:71304/149‑1 (MQ=255)
                                        tcatcTGCTGAGCTGGGAAAGTGTGTCCGATGACGCGGTCAATCAGCGCGTGCTCGACATCATCAAGGCCGTTCGTGAACGCGGTGATGCGGCGCTGGTGGAATTCACCCAGCGTTTCGATGGTGTTGACGCCAAGTCGATTGACGAcc                                                                                                              <  2:130688/149‑1 (MQ=255)
                                         catcTGCTGAGCTGGGAAAGTGTGTCCGATGACGCGGTCAATCAGCGCGTGCTCGACATCATCAAGGCCGTGCGTGAACGCGGTGATGCGGCGCTGGTGGAATTCACCCAGCGTTTCGATGGTGTTGACGCCAAGTCGATTGACGACCt                                                                                                             >  2:47094/1‑149 (MQ=255)
                                          atcTGCTGAGCTGGGAAAGTGTGTCCGATGACGCGGTCAATCAGCGCGTGCTCGACATCATCAAGGCCGTGCGTGAACGCGGTGATGCGGCGCTGGTGGAATTCACCCAGCGTTTCGATGGTGTTGACGCCAAGTCGATTGACGACCt                                                                                                             <  2:285698/148‑1 (MQ=255)
                                             tgctgAGCTGGGAAAGTGTGTCCGATGACGCGGTCAATCAGCGCGTGCTCGACATCATCAAGGCCGTGCGTGAACGCGGTGATGCGCCGCTGGTGGAATTCACCCAGCGTTTCGATGgt                                                                                                                                       >  7:177804/1‑119 (MQ=255)
                                             tgctgAGCTGGGAAAGTGTGTCCGATGACGCGGTCAATCAGCGCGTGCTCGACATCATCAAGGCCGTGCGTGAACGCGGTGATGCGACGCTGGTGGAATTCACCCAGCGTTTCGATGgt                                                                                                                                       <  8:177804/119‑1 (MQ=255)
                                              gctgAGCTGGGAAAGTGTGTCCGATGACGCGGTCAATCAGCGCGTGCTCGACATCATCAAGGCCGTGCGTGAACGCGGTGATGCGGCGCTGGTGGAATTCCCCCAGCGTTTCGAGGGTGTTGACGCCAAGTCGATTGACGACCTGATcc                                                                                                        >  6:95733/1‑149 (MQ=255)
                                               ctgAGCTGGGAAAGTGTGTCCGATGACGCGGTCAATCAGCGCGTGCTCGACATCATCAAGGCCGTGCGTGAACGCGGTGATGCGGCGCTGGTGGAATTCACCCAGCGTTTCGATGGTGTTGACGCCAAGTCGATTGACGACCTGATCCt                                                                                                       >  8:145470/1‑149 (MQ=255)
                                                       ggAAAGTGTGTCCGATGACGCGGTCAATCAGCGCGTGCTCGACATCATCAAGGCCGTGCGTGAACATCGGTGATGCGGCGCTGGTGGAATTCACCCAGCGTTTCGATGGTGTTGACGCCAAGTCGATTGACGACCTGATCCTCGATcgc                                                                                                >  1:94652/1‑149 (MQ=255)
                                                               tgtCCGATGACGCGGTCAATCAGCGCGTGCTCGACATCATCAAGGCCGTGCGTGAACGCGGTGATGCGGCGCTGGTGGAATTCACCCAGCGTTTCGATGGTGTTGACGCCAAGTCGATTGACGACCTGATCCTCGATCGCGCGCGCCt                                                                                       <  8:16818/148‑1 (MQ=255)
                                                                 tCCGATGACGCGGTCAATCAGCGCGTGCTCGACATCATCAAGGCCGTGCGTGAACGCGGTGATGCGGCGCTGGTGGAATTCACCCAGCGTTTCGATGGTGTTGACGCCAAGTc                                                                                                                         >  8:191351/1‑113 (MQ=255)
                                                                 tCCGATGACGCGGTCAATCAGCGCGTGCTCGACATCATCAAGGCCGTGCGTGAACGCGGTGATGCGGCGCTGGTGGAATTCACCCAGCGTTTCGATGGTGTTGACGCCAAGTc                                                                                                                         <  7:191351/113‑1 (MQ=255)
                                                                 tCCGATGACGCGGTCAATCAGCGCGTGCTCGACATCATCAAGGCCGTGCGTGAACGCGGTGATGCGGCGCTGGTGGAATTCACCCAGCGTTTCGATGGTGTTGACGCCAAGTCGATTGACGACCTGATCCTCGATCGCGCGCGCCTgg                                                                                     >  6:270829/1‑148 (MQ=255)
                                                                    gATGACGCGGTCAATCAGCGCGTGCTCGACATCATCAAGGCCGTGCGTGAACGCGGTGATGCGGCGCTGGTGGAATTCACCCAGCGTTTCGATGGTGTTGACGCCAAg                                                                                                                           >  5:179713/1‑108 (MQ=255)
                                                                    gATGACGCGGTCAATCAGCGCGTGCTCGACATCATCAAGGCCGTGCGTGAACGCGGTGATGCGGCGCTGGTGGAATTCACCCAGCGTTTCGATGGTGTTGACGCCAAg                                                                                                                           <  6:179713/108‑1 (MQ=255)
                                                                     aTGACGCGGTCAATCAGCGCGTGCTCGACATCATCAAGGCCGTGCGTGAACGCGGTGATGCGGCGCTGGTGGAATTCACCCAGCGTTTCGATGGTGTTGACGCCAAGTCGATTga                                                                                                                   >  7:16192/1‑115 (MQ=255)
                                                                     aTGACGCGGTCAATCAGCGCGTGCTCGACATCATCAAGGCCGTGCGTGAACGCGGTGATGCGGCGCTGGTGGAATTCACCCAGCGTTTCGATGGTGTTGACGCCAAGTCGATTga                                                                                                                   <  8:16192/115‑1 (MQ=255)
                                                                     aTGACGCGGTCAATCAGCGCGTGCTCGACATCATCAAGGCCGTGCGTGAACGCGGTGATGCGGCGCTGGTGGAATTCACCCAGCGTTTCGATGGTGTTGACGCCAAGTCGATTGACGACCTGATCCTCGATCGCGCGCGCCTGGAACTg                                                                                <  6:217978/149‑1 (MQ=255)
                                                                         cgcgGTCAATCAGCGCGTGCTCGACATCATCAAGGCCGTGCGTGAACGCGGAGATGCGGCGCTGGTGGAATTCACCCAGCGTTTCGATGGTGATGACGCCAAGTCGATTGACGACCTGATCCTCGATCGCGCGCGCCTGGAACTGgcgc                                                                            >  5:218361/1‑149 (MQ=255)
                                                                             gTCAATCAGCGCGTGCTCGACATCATCAAGGCCGTGCGTGAACGCGGTGATGCGGCGCTGGTGGAATTCACCCAGCGTTTCGATGGTGTTGACGCCAAGTCGATTGACGACCTGATCCTCGATCGCGCGCGCCTGGAACTGGCGCTGa                                                                         >  4:247184/1‑148 (MQ=255)
                                                                             gTCAATCAGCGCGTGCTCGACATCATCAAGGCCGTGCGTGAACGCGGTGATGCGGCGCTGGTGGAATTCACCCAGCGTTTCGATGGTGTTGACGCCAAGTCGATTGACGACCTGATCCTCGATCGCGCGCGCCTGGAACTGGCGCTGa                                                                         <  3:247184/148‑1 (MQ=255)
                                                                             gTCAATCAGCGCGTGCTCGACATCATCAAGGCCGTGCGTGAACGCGGTGATGCGGCGCTGGTGGAATTCACCCAGCGTTTCGATGGTGTTGACGCCAAGTCGATTGACGACCTGATCCTCGATCGCGCGCGCCTGGAACTGGCGCTGAc                                                                        <  7:145470/149‑1 (MQ=255)
                                                                                 aTCAGCGCGTGCTCGACATCATCAAGGCCGTGCGTGAACGCGGTGATGCGGCGCTGGTGGAATTCACCCAGCGTTTCGATGGTGTTGACGCCAAGTCGATTGACGACCTGATCCTCGATCGCGCGCGCCTGGAACTGGCGCTGACCCGt                                                                    <  8:28714/149‑1 (MQ=255)
                                                                                           gCTCGACATCATCAAGGCCGTGCGTGAACGCGGTGATGCGGCGCTGGTGGAATTCACCCAGCGTTTCGATGGTGTTGACGCCAAGTCGATTGACGACCTGATCCTCGATCGCGCGCGCCTGGAACTGGCGCTGACCCGTATTACCCCgg                                                          >  2:204013/1‑149 (MQ=255)
                                                                                            cTCGACATCATCAAGGCCGTGCGTGAACGCGGTGATGCGGCGCTGGTGGAATTCACCCAGCGTTTCGATGGTGTTGACGCCAAGTCGATTGACGACCTGATCCTCGATCGCGACGCGCCTGGAACTGGCGCTGACCCGTATTACCCCg                                                           <  1:47094/148‑1 (MQ=255)
                                                                                             tCGACATCATCAAGGCCGTGCGTGAACGCGGTGATGCGGCGCTGGTGGAATTCACCCAGCGTTTCGATGGTGTTGACGCCAAGTCGATTGACGACCTGATCCTCGATCGCGCGCGCCTGGAACTGGCGCTGACCCGTATTACCCCGGTc                                                        <  3:157287/149‑1 (MQ=255)
                                                                                               gACATCATCAAGGCCGTGCGTGAACGCGGTGATGCGGCGCTGGTGGAATTCACCCAGCGTTTCGATGGTGTTGACGCCAAGTCGATTGACGACCTGATCCTCGATCGCGCGCGCCt                                                                                       <  6:178442/116‑1 (MQ=255)
                                                                                               gACATCATCAAGGCCGTGCGTGAACGCGGTGATGCGGCGCTGGTGGAATTCACCCAGCGTTTCGATGGTGTTGACGCCAAGTCGATTGACGACCTGATCCTCGATCGCGCGCGCCt                                                                                       >  5:178442/1‑116 (MQ=255)
                                                                                                      tcaAGGCCGTGCGTGAACGCGGTGATGCGGCGCTGGTTGGAATTCACCCAGCGTTTCGATGGTGTTGACGCCAAGTCGATTGACGACCTGATCCTCGATCGCGCGCGCCTGGAACTGGCGCTGACCCGTATTACCCCGGTCCAgcgc                                                  >  1:56290/1‑147 (MQ=255)
                                                                                                       caAGGCCGTGCGTGAACGCGGTGATGCGGCGCTGGTGGAATTCACCCAGCGTTTCGATGGTGTTGAc                                                                                                                                 <  7:155985/67‑1 (MQ=255)
                                                                                                       caAGGCCGTGCGTGAACGCGGTGATGCGGCGCTGGTGGAATTCACCCAGCGTTTCGATGGTGTTGAc                                                                                                                                 >  8:155985/1‑67 (MQ=255)
                                                                                                       caAGGCCGTGCGTGAACGCGGTGATGCGGCGCTGGTGGAATTCACCCAGCGTTTCGATGGTGTTGACGCCAAGTCGATTGACGACCTGATCCTCGATCGCGCGCGCCTGGAACTGGCGCTGACCCGTATTACCCCGGTCCAGCGCGAgg                                              >  7:248670/1‑149 (MQ=255)
                                                                                                        aaGGCCGTGCGTGAACGCGGTGATGCGGCGCTGGTGGAATTCACCCAGCGTTTCGATGGTGTTGACGCCAAGTCGATTGACGACCTGATCCTCGATCGCGCGCGCCTGGAACTGGCGCTGACCCGTATTACCCCGGTCCAGCGCGAGGc                                             >  7:227329/1‑149 (MQ=255)
                                                                                                            ccGTGCGTGAACGCGGTGATGCGGCGCTGGTGGAATTCACCCAGCGTTTCGATGGTGTTGACGCCAAGTCGATTGACGACCTGATCCTCGATCGCGCGCGCCTGGAACTGGCGCTGACCCGTATTACCCCGGTCCAGCGCGAGGCCCTg                                         <  1:204013/149‑1 (MQ=255)
                                                                                                             cgtgcgtgAACGCGGTGATGCGGCGCTGGTGGAATTCACCCAGCGTTTCGATGGTGTTGACGCCAAGTCGATTGACGACCTGATCCTCGATCGCGCGCGCCTGGAACTGGCGCTGACCCGTATTACCCCGGTCCAGCGCGAGGCCCTg                                         >  7:120801/1‑148 (MQ=255)
                                                                                                                 cgtgAACGCGGTGATGCGGCGCTGGTGGAATTCACCCAGCGTTTCGATGGTGTTGACGCCAAGTCGATTGACGACCTGATCCTCGATCGCGCGCGCCTGGAACTGGCGCTGACCCGTATTACCCCGGTCCAGCGCGAGGCCCTGGaaaa                                    <  8:120801/149‑1 (MQ=255)
                                                                                                                   tgAACGCGGTGATGCGGCGCTGGTGGAATTCACCCAGCGTTTCGATGGTGTTGACGCCAAGTCGATTGACGACCTGATCCTCGATCGCGCGCGCCTGGAACTGGCGCTGACCCGTATTACCCCGGTCCAGCGCGAGGCCCTGGAAAAgg                                  >  7:204301/1‑149 (MQ=255)
                                                                                                                    gAACGCGGTGATGCGGCGCTGGTGGAATTCACCCAGCGTTTCGATGGTGTTGACGCCAAGTCGATTGACGACCTGATCCTCGATCGCGCGCGCCTGGAACTGGCGCTGACCCGTATTACCCCGGTCCAGCGCGAGGCCCTGGAAAAGgc                                 >  3:102154/1‑149 (MQ=255)
                                                                                                                           agTGATGCGGCGCTGGTGGAATTCACCCAGCGTTTCGATGGTGTTGACGCCAAGTCGATTGACGACCTGATCCTCGATCGCGCGCGCCTGGAACTGGCGCTGACCCGTATTACCCCGGTCCAGCGCGAGGCCCTGGAAAAGGCCGCCaa                           >  2:49908/2‑149 (MQ=255)
                                                                                                                             tGATGCGGCGCTGGTGGAATTCACCCAGCGTTTCGATGGTGTTGACGCCAAGTCGATTGACGACCTGATCCTCGATCGCGCGCGCCTGGAACTGGCGCTGACCCGTATTACCCCGGTCCAGCGCGAGGCCCTGGAAAAGGCCGCCAAcc                         <  8:227329/149‑1 (MQ=255)
                                                                                                                                  aaacgcTGGGTGAATTCCACCAGCGTTTCGATGGTGTTGACGCCAAGTCGATTGACGACCTGATCCTCGATCGCGCGCGCCTGGAACTGGCGCTGACCCGTATTACCCCGGTCCAGCGCGAGGCCCTGGAAAAGGCCGCCAACCGtgtg                    >  8:192432/4‑149 (MQ=255)
                                                                                                                                              gAATTCACCCAGCGTTTCGATGGTGTTGACGCCAAGTCGATTGACGACCTGATCCTCGATCGCGCGCGCCTGGAACTGGCGCTGACCCGTATTACCCCGGTCCAGCGCGAGGCCCTGGAAAAGGCCGCCAACCGTGTGCGCATTTATc          <  2:56290/148‑1 (MQ=255)
                                                                                                                                              gAATTCACCCAGCGTTTCGATGGTGTTGACGCCAAGTCGATTGACGACCTGATCCTCGATCGCGCGCGCCTGGAACTGGCGCTGACCCGTATTACCCCGGTCCAGCGCGAGGCCCTGGAAAAGGCCGCCAACCGTGTGCGCATTTATc          <  8:204301/148‑1 (MQ=255)
                                                                                                                                                    aCCCAGCGTTTCGATGGTGTTGACGCCAAGTCGATTGACGACCTGATCCTCGATCGCGCGCGCCTGGAACTGGCGCTGACCCGTATTACCCCGGTCCAGCGCGAGGCCCTGGAAAAGGCCGCCAACCGTGTGCGCATTTATc          <  3:54988/142‑1 (MQ=255)
                                                                                                                                                    aCCCAGCGTTTCGATGGTGTTGACGCCAAGTCGATTGACGACCTGATCCTCGATCGCGCGCGCCTGGAACTGGCGCTGACCCGTATTACCCCGGTCCAGCGCGAGGCCCTGGAAAAGGCCGCCAACCGTGTGCGCATTTATc          >  4:54988/1‑142 (MQ=255)
                                                                                                                                                     cccAGCGTTTCGATGGTGTTGACTCCAAGTCGATTGACGACCTGATCCTCGATCGCGCGCGCCTGGAACTGGCGCTGACCCGTATTACCCCGGTCCAGCGCGAGGCCCTGGAAAAGGCCGCCAACCGTGTGCGCATTTATCACGAGCgg  <  5:95733/149‑1 (MQ=255)
                                                                                                                                                     cccAGCGTTTCGATGGTGTTGACGCCAAGTCGATTGACGACCTGATCCTCGATCGCGCGCGCCTGGAACTGGCGCTGACCCGTATTACCCCGGTCCAGCGCGAGGCCCTGGAAAAGGCCGCCAACCGTGTGCGCATTTATCACGAGCgg  <  5:109193/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTCTCAACGCCGCTGATCCGGATTTCGCCCGACATCTGGATCATCTGCTGAGCTGGGAAAGTGTGTCCGATGACGCGGTCAATCAGCGCGTGCTCGACATCATCAAGGCCGTGCGTGAACGCGGTGATGCGGCGCTGGTGGAATTCACCCAGCGTTTCGATGGTGTTGACGCCAAGTCGATTGACGACCTGATCCTCGATCGCGCGCGCCTGGAACTGGCGCTGACCCGTATTACCCCGGTCCAGCGCGAGGCCCTGGAAAAGGCCGCCAACCGTGTGCGCATTTATCACGAGCGG  >  NC_002947/1105546‑1105841

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