Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,112,355 C→T 5.4% L19L (CTG→CTA rluB ← 23S rRNA pseudouridylate synthase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,112,3550CT 5.4% 228.0 / 2.1 74L19L (CTG→CTArluB23S rRNA pseudouridylate synthase B
Reads supporting (aligned to +/- strand):  ref base C (33/37);  new base T (2/2);  total (35/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.94e-01

TGCCGTCGACGGCGATGGCGTCGTGCAGGTCGACGCGCTGGCCAAGGGTGGCCTCGACGCCGTTGACCTTGATGCGACCCTGGCTGATCCAGGCCTCGACGTCACGCCGCGAGCCCACGCCAATGCGCGCCAGCACTTTCTGCAGCTTTTCGCCGGATGGCGGGGTGATTTCGGTTTCTTTCAGGTCTTGCTCACTCATCTGGGCACCTCCCGGTGTAGGAATGAAAACAGGTTTCTGGCGACGAACGCGCCAAAAGGCCGCGCATCATACGCGGTTCATCAGG  >  NC_002947/5112213‑5112496
                                                                                                                                                 |                                                                                                                                               
tGCCGTCGACGGCGATGGCGTCGTGCAGGTCGACGCGCTGGCCAAGGGTGGCCTCGACGCCGTTGACCTTGATGCGACCCTGGCTGATCCAGGCCTCGACGTCACGCCGCGAGCCCACGCCAATGCGCGCCAGCACTTTCTGCAGCttt                                                                                                                                           >  4:328763/1‑149 (MQ=255)
 gCCGTCGACGGCGATGGCGTCGTGCAGGTCGACGCGCTGGCCAAGGGTGGCCTCGACGCCGTTGACCTTGATGCGACCCTGGCTGATCCAGGCCTCGACGTCACGCCGCGAGCCCACGCCAATGCGCGCCAGCACTTTCTGCAGCtttt                                                                                                                                          >  5:75095/1‑149 (MQ=255)
   cGTCGACGGCGATGGCGTCGTGCAGGTCGACGCGCTGGCCAAGGGTGGCCTCGACGCCGTTGACCTTGATGCGACCCTGGCTGATCCAGGCCTCGACGTCACGCCGCGAGCCCACGCCAATGCGCGCCAGCACTTTCTGCAGCTTTTCg                                                                                                                                        >  7:250803/1‑149 (MQ=255)
       gACGGCGATGGCGTCGTGCAGGTCGACGCGCTGGCCAAGGGTGGCCTCGACGCCGTTGACCTTGATGCGACCCTGGCTGATCCAGGCCTCGACGTCACGCCGCGAGCCCACGCCAATGCGCGCCAGCACTTTCTGCAGCTTTTCGCCg                                                                                                                                     >  3:310453/1‑148 (MQ=255)
       gACGGCGATGGCGTCGTGCAGGTCGACGCGCTGGCCAAGGGTGGCCTCGACGCCGTTGACCTTGATGCGACCCTGGCTGATCCAGGCCTCGACGTCACGCCGCGAGCCCACGCCAATGCGCGCCAGCACTTTCTGCAGCTTTTCGCCg                                                                                                                                     <  1:128467/148‑1 (MQ=255)
       gACGGCGATGGCGTCGTGCAGGTCGACGCGCTGGCCAAGGGTGGCCTCGACGCCGTTGACCTTGATGCGACCCTGGCTGATCCAGGCCTCGACGTCACGCCGCGAGCCCACGCCAATGCGCGCCAGCACTTTCTGCAGCTTTTCGCCg                                                                                                                                     >  2:67731/1‑148 (MQ=255)
             gATGGCGTCGTGCAGGTCGACGCGCTGGCCAAGGGTGGCCTCGACGCCGTTGACCTTGATGCGACCCTGGCTGATCCAGGCCTCGACGTCACGCCGCGAGCCCACGCCAATGCGCGCCAGCACTTTCTGCAGCTTTTCGCCGGATGGCg                                                                                                                              >  3:364530/1‑149 (MQ=255)
             gATGGCGTCGTGCAGGTCGACGCGCTGGCCAAGGGTGGCCTCGACGCCGTTGACCTTGATGCGACCCTGGCTGATCCAGGCCTCGACGTCACGCCGCGAGCCCACGCCAATGCGCGCCAGCACTTTCTGCAGCTTTTCGCCGGATGGCg                                                                                                                              >  1:252495/1‑149 (MQ=255)
             gATGGCGTCGTGCAGGGCGACGCGCTGGCCAATGGTGGCCTCGACGCCGTTGCATCCTTGTAGCCCCCCCACGCTGATCCAGGCCTCGACGTCACGCCGCGAGCCCACGCCAATGCGCGCCAGCACTTTCTGCAGCTTTTCGCCGGATg                                                                                                                                 <  1:122244/149‑1 (MQ=255)
                   gtcgtGCAGGTCGACGCGCTGGCCAAGGGTGGCCTCGACGCCGTTGACCTTGATGCGACCCTGGCTGATCCAGGCCTCGACGTCACGCCGCGAGCCCACGCCAATGCGCGCCAGCACTTTCTGCAGCTTTTCGCCGGATGGCGGGGTGa                                                                                                                        >  8:292872/1‑149 (MQ=255)
                     cgtGCAGGTCGACGCGCTGGCCAAGGGTGGCCTCGACGCCGTTGACCTTGATGCGACCCTGGCTGATCCAGGCCTCGACGTCACGCCGCGAGCCCACGCCAATGCGCGCCAGCACTTTCTGCAGCTTTTCGCCGGATGGCGGGGTGAtt                                                                                                                      <  1:246827/149‑1 (MQ=255)
                         cAGGTCGACGCGCTGGCCAAGGGTGGCCTCGACGCCGTTGACCTTGATGCGACCCTGGCTGATCCAGGCCTCGACGTCACGCCGCGAGCCCACGCCAATGCGCGCCAGCACTTTCTGCAGCTTTTCGCCGGATGGCGGGGTGATTTCgg                                                                                                                  <  8:250803/149‑1 (MQ=255)
                          aGGTCGACGCGCTGGCCAAGGGTGGCCTCGACGCCGTTGACCTTGATGCGACCCTGGCTGATCCAGGCCTCGACGTCACGCCGCGAGCCCACGCCAATGCGCGCCAGCACTTTCTGCAGCTTTTCGCCGGATGGCGGGGTGATTTCGGt                                                                                                                 >  6:195510/1‑149 (MQ=255)
                           ggTCGACGCGCTGGCCAAGGGTGGCCTCGACGCCGTTGACCTTGATGCGACCCTGGCTGATCCAGGCCTCGACGTCACGCCGCGAGCCCACGCCAATGCGCGCCAGCACTTTCTGCAGCTTTTCGCCGGATGGCGGGGTGATTTCGGtt                                                                                                                <  7:302441/149‑1 (MQ=255)
                               gACGCGCTGGCCAAGGGTGGCCTCGACGCCGTTGACCTTGATGCGACCCTGGCTGATCCAGGCCTCGACGTCACGCCGCGAGCCCACGCCAATGCGCGCCAGCACTTTCTGTAGCTTTTCGCCGGATGGCGGGGTGATTTCGGtttctt                                                                                                            >  7:229263/1‑149 (MQ=255)
                               gACGCGCTGGCCAAGGGTGGCCTCGACGCCGTTGACCTTGATGCGACCCTGGCTGATCCAGGCCTCGACGTCACGCCGCGAGCCCACGCCAATGCGCGCCAGCACTTTCTGTAGCTTTTCGCCGGATGGCGGGGTGATTTCGGtttctt                                                                                                            >  7:223447/1‑149 (MQ=255)
                                 cgcgCTGGCCAAGGGTGGCCTCGACGCCGTTGACCTTGATGCGACCCTGGCTGATCCAGGCCTCGACGTCACGCCGCGAGCCCACGCCAATGCGCGCCAGCACTTTCTGCAGCTTTTCGCCGGATGGCGGGGTGATTTCGGtttctttc                                                                                                         <  3:328763/149‑1 (MQ=255)
                                 cgcgCTGGCCAAGGGTGGCCTCGACGCCGTTGACCTTGATGCGACCCTGGCTGATCCAGGCCTCGACGTCACGCCGCGAGCCCACGCCAATGCGCGCCAGCACTTTCTGCAGCTTTTCGCCGGATGGCGGGGTGATTTCGGtttctttc                                                                                                         <  5:48431/149‑1 (MQ=255)
                                   cgcTGGCCAAGGGTGGCCTCGACGCCGTTGACCTTGATGCGACCCTGGCTGATCCAGGCCTCGACGTCACGCCGCGAGCCCACGCCAATGCGCGCCAGCACTTTCTGCAGCTTTTCGCCGGATGGCGGGGTGATTTTGGTTTCTTTCAg                                                                                                       <  3:299254/149‑1 (MQ=255)
                                   cgcTGGCCAAGGGTGGCCTCGACGCCGTTGACCTTGATGCGACCCTGGCTGATCCAGGCCTCGACGTCACGCCGCGAGCCCACGCCAATGCGCGCCAGCACTTTCTGCAGCTTTTCGCCGGATGGCGGGGTGATATCGGTTTCTTTCAg                                                                                                       <  4:50861/149‑1 (MQ=255)
                                    gcTGGCCAAGGGTGGCCTCGACGCCGTTGACCTTGATGCGACCCTGGCTGATCCAGGCCTCGACGTCACGCCGCGAGCCCACGCCAATGCGCGCCAGCACTTTCTGCAGCTTTTc                                                                                                                                         >  2:116783/1‑115 (MQ=255)
                                    gcTGGCCAAGGGTGGCCTCGACGCCGTTGACCTTGATGCGACCCTGGCTGATCCAGGCCTCGACGTCACGCCGCGAGCCCACGCCAATGCGCGCCAGCACTTTCTGCAGCTTTTc                                                                                                                                         <  1:116783/115‑1 (MQ=255)
                                    gcTGGCCAAGGGTGGCCTCGACGCCGTTGACCTTGATGCGACCCTGGCTGATCCAGGCCTCGACGTCACGCCGCGAGCCCACGCCAATGCGCGCCAGCACTTTCTGCAGCTTTTCGCCGGATGGCGGGGTGATTTCGGtttctttc                                                                                                         >  5:7959/1‑146 (MQ=255)
                                    gcTGGCCAAGGGTGGCCTCGACGCCGTTGACCTTGATGCGACCCTGGCTGATCCAGGCCTCGACGTCACGCCGCGAGCCCACGCCAATGCGCGCCAGCACTTTCTGCAGCTTTTCGCCGGATGGCGGGGTGATTTCGGtttctttc                                                                                                         <  6:7959/146‑1 (MQ=255)
                                    gcTGGCCAAGGGTGGCCTCGACGCCGTTGACCTTGATGCGACCCTGGCTGATCCAGGCCTCGACGTCACGCCGCGAGCCCACGCCAATGCGCGCCAGCACTTTCTGCAGCTTTTCGCCGGATGGCGGGGTGATTTCGGtttcttt                                                                                                           <  8:12843/145‑1 (MQ=255)
                                    gcTGGCCAAGGGTGGCCTCGACGCCGTTGACCTTGATGCGACCCTGGCTGATCCAGGCCTCGACGTCACGCCGCGAGCCCACGCCAATGCGCGCCAGCACTTTCTGCAGCTTTTCGCCGGATGGCGGGGTGATTTCGGtttcttt                                                                                                           >  7:12843/1‑145 (MQ=255)
                                        gCCAAGGGTGGCCTCGACGCCGTTGACCTTGATGCGACCCTGGCTGATCCAGGCCTCGACGTCACGCCGCGAGCCCACGCCAATGCGCGCCAGCACTTTCTGCAGCTTTTCGCCGGATGGCGGGGTGATTTCGGTTTCTTTCAGGTCtt                                                                                                  <  7:292872/149‑1 (MQ=255)
                                          cAAGGGTGGCCTCGACGCCGTTGACCTTGATGCGACCCTGGCTGATCCAGGCCTCGACGTCACGCCGCGAGCCCACGCCAATGCGCGCCAGCACTTTCTGCAGCTTTTCGCCGGATGGCGGGGTGATTTCGGTTTCTTTCAGGTCTTGc                                                                                               >  1:192080/1‑149 (MQ=255)
                                          cAAGGGTGGCCTCGACGCCGTTGACCTTGATGCGACCCTGGCTGATCCAGGCCTCGACGTCACGCCGCGAGCCCACGCCAATGCGCGCCAGCACTTTCTGCAGCTTTTCGCCGGATGGCGGGGTGATTTCGGTTTCTTTCAGGTCTTGc                                                                                               >  2:160925/1‑149 (MQ=255)
                                          cAAGGGTGGCCTCGACGCCGTTGACCTTGATGCGACCCTGGCTGATCCAGGCCTCGACGTCACGCCGCGAGCCCACGCCAATGCGCGCCAGCACTTTCTGCAGCTTTTCGCCGGATGGCGGGGTGATTTCGGTTTCTTTCAGGTCTTGc                                                                                               <  8:115991/149‑1 (MQ=255)
                                            aGGGTGGCCTCGACGCCGTTGACCTTGATGCGACCCTGGCTGATCCAGGCCTCGACGTCACGCCGCGAGCCCACGCCAATGCGCGCCAGCACTTTCTGCAGCTTTTCGCCGGATGGCGGGGTGATTTCGGTTTCTTTCAGGTCTTGct                                                                                              >  2:253518/1‑148 (MQ=255)
                                                  gCCTCGACGCCGTTGACCTTGATGCGACCCTGGCTGATCCAGGCCTCGACGTCACGCCGCGAGCCCACGCCAATGCGCGCCAGCACTTTCTGCAGCTTTTCGCCGGATGGCGGGGTGATTTCGGTTTCTTTCAGGTCTTGctcactca                                                                                        >  7:123677/1‑148 (MQ=255)
                                                  gCCTCGACGCCGTTGACCTTGATGCGACCCTGGCTGATCCAGGCCTCGACGTCACGCCGCGAGCCCACGCCAATGCGCGCCAGCACTTTCTGCAGCTTTTCGCCGGATGGCGGGGTGATTTCGGTTTCTTTCAGGTCTTGCTCACTCAt                                                                                       >  3:369011/1‑149 (MQ=255)
                                                   ccTCGACGCCGTTGACCTTGATGCGACCCTGGCTGATCCAGGCCTCGACGTCACGCCGCGAGCCCACGCCAATGCGCGCCAGCACTTTCTGCAGCTTTTCGCCGGATGGCGGGGTGATTTCGGTTTCTTTCAGGTCTTGctca                                                                                            <  8:290849/143‑1 (MQ=255)
                                                   ccTCGACGCCGTTGACCTTGATGCGACCCTGGCTGATCCAGGCCTCGACGTCACGCCGCGAGCCCACGCCAATGCGCGCCAGCACTTTCTGCAGCTTTTCGCCGGATGGCGGGGTGATTTCGGTTTCTTTCAGGTCTTGctca                                                                                            >  7:290849/1‑143 (MQ=255)
                                                   ccTCGACGCCGTTGACCTTGATGCGACCCTGGCTGATCCAGGCCTCGACGTCACGCCGCGAGCCCACGCCAATGCGCGCCAGCACTTTCTGCAGCTTTTCGCCGGATGGCGGGGTGATTTCGGTTTCTTTCAGGTCTTGCTCACTCATc                                                                                      <  4:289222/149‑1 (MQ=255)
                                                     tCGACGCCGTTGACCTTGATGCGACCCTGGCTGATCCAGGCCTCGACGTCACGCCGCGAGCCCACGCCAATGCGCGCCAGCACTTTCTGCAGCTTTTCGCCGGATGGCGGGGTGATTTCGGTTTCTTTCAGGTCTTGCTCACTCATCTg                                                                                    >  1:217003/1‑149 (MQ=255)
                                                       gACGCCGTTGACCTTGATGCGACCCTGGCTGATCCAGGCCTCGACGTCACGCCGCGAGCCCACGCCAATGCGCGCCAGCACTTTCTGCAGCTTTTCGCCGGATGGCGGGGTGATTTCGGTTTCTTTCAGGTCTTGCTCACTCATCTggg                                                                                  <  6:245562/149‑1 (MQ=255)
                                                       gACGCCGTTGACCTTGATGCGACCCTGGCTGATCCAGGCCTCGACGTCACGCCGCGAGCCCACGCCAATGCGCGCCAGCACTTTCTGCAGCTTTTCGCCGGATGGCGGGGTGATTTCGGTTTCTTTCAGGTCTTGCTCACTCATCTgg                                                                                   <  5:375610/148‑1 (MQ=255)
                                                          gCCGTTGACCTTGATGCGACCCTGGCTGATCCAGGCCTCGACGTCACGCCGCGAGCCCACGCCAATGCGCGCCAGCACTTTCTGCAGCTTTTCGCCGGATGGCGGGGTGATTTCGGTTTCTTTCAGGTCTTGCTCACTCATCTGGGCAc                                                                               <  2:252495/149‑1 (MQ=255)
                                                          gCCGTTGACCTTGATGCGACCCTGGCTGATCCAGGCCTCGACGTCACGCCGCGAGCCCACGCCAATGCGCGCCAGCACTTTCTGCAGCTTTTCGCCGGATGGCGGGGTGATTTCGGTTTCTTTCAGGTCTTGCTCACTCATCTGGGCAc                                                                               >  2:113271/1‑149 (MQ=255)
                                                          cccGTGGACCTGGACGCGCCCCTGCCTGATCCAGGCCTCGACGTCACGCCGCGAGCCCACGCCAATGCGCGCCAGCACTTTCTGCAGCTTTTCGCCGGATGGCGGGGTGATTTCGGTTTCTTTCAGGTCTTGCTCACTCATCTGGGCAc                                                                               <  4:310453/148‑1 (MQ=255)
                                                            cGTTGACCTTGATGCGACCCTGGCTGATCCAGGCCTCGACGTCACGCCGCGAGCCCACGCCAATGCGCGCCAGCACTTTCTGCAGCTTTTCGCCGGATGGCGGGGTGATTTCGGTTTCTTTCa                                                                                                        <  2:241249/123‑1 (MQ=255)
                                                            cGTTGACCTTGATGCGACCCTGGCTGATCCAGGCCTCGACGTCACGCCGCGAGCCCACGCCAATGCGCGCCAGCACTTTCTGCAGCTTTTCGCCGGATGGCGGGGTGATTTCGGTTTCTTTCa                                                                                                        >  1:241249/1‑123 (MQ=255)
                                                                     cTTGATGCGACCCTGGCTTATCCAGGCCTCGACGTCACGCCGCGAGCCCACGCCAATGCGCGCCAGCACTTTCAGCAGCTTTTCGCCGGATGGCGGGGTGATTTCGGTTTCTTTCAGGTCTTGCTCACTCATCTGGGCACCTCCCGgt                                                                       >  5:175456/1‑148 (MQ=255)
                                                                        gATGCGACCCTGGCTGATCCAGGCCTCGACGTCACGCCGCGAGCCCACGCCAATGCGCGCCAGCACTTTCTGCAGCTTTTCGCCGGATGGCGGGGTGATTTCGGTTTCTTTCAGGTCTTGCTCACTCATCTGGGCACCTCCCGGTGTa                                                                    <  6:75095/148‑1 (MQ=255)
                                                                        gATGCGACCCTGGCTGATCCAGGCCTCGACGTCACGCCGCGAGCCCACGCCAATGCGCGCCAGCACTTTCTGCAGCTTTTCGCCGGATGGCGGGGTGATTTCGGTTTCTTTCAGGTCTTGCTCACTCATCTGGGCACCTCCCGGTGTa                                                                    <  4:123225/148‑1 (MQ=255)
                                                                        gATGCGACCCTGGCTGATCCAGGCCTCCACGGCACGCCGCGAGCCCACGCCAATGCGCGCCAGCACTTTCTGCAGCTTTTCGCCGGATGGCGGGGTGATTTCGGTTTCTTTCCAGGTCTTGCTCACTCATCTGGGCACCTCCCGgtgt                                                                     <  1:253518/148‑1 (MQ=255)
                                                                           tGCGACCCTGGCTGATCCAGGCCTCGACGTCACGCCGCGAGCCCACGCCAATGCGCGCCAGCACTTTCTGCAGCTTTTCGCCGGATGGCGGGGTGATTTCGGTTTCTTTCAGGTCTTGCTCACTCATCTGGGCACCTCCCGGTGTAgg                                                                  >  5:292774/1‑148 (MQ=255)
                                                                            gCGACCCTGGCTGATCCAGGCCTCGACGTCACGCCGCGAGCCCACGCCAATGCGCGCCAGCACTTTCTGCAGCTTTTCGCCGGATGGCGGGGTGATTTCGGTTTCTTTCAGGTCTTGCTCACTCATCTGGGCACCTCCCGGTGTAGGaa                                                                >  6:189993/1‑149 (MQ=255)
                                                                              gACCCTGGCTGATCCAGGCCTCGACGTCACGCCGCGAGCCCACGCCAATGCGCGCCAGCACTTTCTGCAGCTTTTCGCCGGATGGCGGGGTGATTTCGGTTTCTTTCAGGTCTTGCTCACTCATCTGGGCACCTCCCGGTGTAGGAATg                                                              <  6:64797/149‑1 (MQ=255)
                                                                              gACCCTGGCTGATCCAGGCCTCGACGTCACGCCGCGAGCCAACGCCAATGCGCGCCAGCACTTTCTGCAGCTTTTCGCCGGATGGCGGGGTGATTTCGGTTTCTTTCAGGTCTTGCTCACTCATCTGGGCACCTCCCGGTGTAGGAAt                                                               >  8:139064/1‑148 (MQ=255)
                                                                               aCCCTGGCTGATCCAGGCCTCGACGTCACGCCGCGAGCCCACGCCAATGCGCGCCAGCACTTTCTGCAGCTTTTCGCCGGATGGCGGGGTGATTTCGGTTTCTTTCAGGTCTTGCTCACTCATCTGGGCACCTCCCgg                                                                        >  5:226755/1‑138 (MQ=255)
                                                                               aCCCTGGCTGATCCAGGCCTCGACGTCACGCCGCGAGCCCACGCCAATGCGCGCCAGCACTTTCTGCAGCTTTTCGCCGGATGGCGGGGTGATTTCGGTTTCTTTCAGGTCTTGCTCACTCATCTGGGCACCTCCCgg                                                                        <  6:226755/138‑1 (MQ=255)
                                                                                 ccTGGCTGATCCAGGCCTCGACGTCACGCCGCGAGCCCACGCCAATGCGCGCCAGCACTTTCTGCAGCTTTTCGCCGGATGGCGGGGTGATTTCGGTTTCTTTCAGGTCTTGCTCACTCATCTGGGCACCTCCCGGTGTAGGAATGaa                                                            >  8:200322/1‑148 (MQ=255)
                                                                                     gCTGATCCAGGCCTCGACGTCACGCCGCGAGCCCACGCCAATGCGCGCCAGCACTTTCTGTAGCTTTTCGCCGGATGGCGGGGTGATTTCGGTTTCTTTCAGGTCTTGCTCACTCATCTGGGCACCTCCCGGTGTAGGAATGAAAACa                                                        <  8:223447/148‑1 (MQ=255)
                                                                                     gCTGATCCAGGCCTCGACGTCACGCCGCGAGCCCACGCCAATGCGCGCCAGCACTTTCTGTAGCTTTTCGCCGGATGGCGGGGTGATTTCGGTTTCTTTCAGGTCTTGCTCACTCATCTGGGCACCTCCCGGTGTAGGAATGAAAACa                                                        <  8:229263/148‑1 (MQ=255)
                                                                                     gCTGATCCAGGCCTCGACGTCACGCCGCGAGCCCACGCCAATGCGCGCCAGCACTTTCTGCAGCTTTTCGCCGGATGGCGGGGTGATTTCGGTTTCTTTCAGGTCTTGCTCACTCATCTGGGCACCTCCCGGTGTAGGAATGAAAAc                                                         <  1:160925/147‑1 (MQ=255)
                                                                                           ccAGGCCTCGACGTCACGCCGCGAGCCCACGCCAATGCGCGCCAGCACTTTCTGCAGCTTTTCGCCGGATGGCGGGGTGATTTCGGttt                                                                                                               >  5:133792/1‑89 (MQ=255)
                                                                                           ccAGGCCTCGACGTCACGCCGCGAGCCCACGCCAATGCGCGCCAGCACTTTCTGCAGCTTTTCGCCGGATGGCGGGGTGATTTCGGttt                                                                                                               <  6:133792/89‑1 (MQ=255)
                                                                                             aGGCCTCGACGTCACGCCGCGAGCCCACGCCAATGCGCGCCAGCACTTTCTGCAGCTTTTCGCCGGATGGCGGGGTGATTTCGGTTTCTTTCAGGTCTTGCTCACTCATCTGGGCACCTCCCGGTGTAGGAATGAAAACAGGTTTCTgg                                               <  5:195510/149‑1 (MQ=255)
                                                                                                 cTCGACGTCACGCCGCGAGCCCACGCCAATGCGCGCCAGCACTTTCTGCAGCTTTTCGCCGGATGGCGGGGTGATTTCGGTTTCTTTCAGGTCTTGCTCACTCATCTGGGCACCTCCCGGTGTAGGAAt                                                               <  1:109404/129‑1 (MQ=255)
                                                                                                 cTCGACGTCACGCCGCGAGCCCACGCCAATGCGCGCCAGCACTTTCTGCAGCTTTTCGCCGGATGGCGGGGTGATTTCGGTTTCTTTCAGGTCTTGCTCACTCATCTGGGCACCTCCCGGTGTAGGAATGAAAACAGGTTTCTGGcgac                                           >  2:289266/1‑149 (MQ=255)
                                                                                                 cTCGACGTCACGCCGCGAGCCCACGCCAATGCGCGCCAGCACTTTCTGCAGCTTTTCGCCGGATGGCGGGGTGATTTCGGTTTCTTTCAGCTCTTGCTCACTCATCTGGGCACCTCCCGGTGTAGGAAt                                                               >  2:109404/1‑129 (MQ=255)
                                                                                                            gccgcGAGCCCACGCCAATGCGCGCCAGCGCTTTCTGCAGCTTTTCGCCGGATGGCGGGGTGATTTCGGTTTGTTTCAGGTCTGGCTCCCTCATCTGGGCACCTGCCGGTGTAGGAATGAAAACAGGTTTCTGGCGCCGCCCGCGCCa                                 >  4:301871/1‑148 (MQ=255)
                                                                                                             ccgcTAGCCCACGCCAATGCGCGCCAGCACTTTCTGCAGCTTTTCGCCGGATGGCGGGGTGATTTCGGTTTCTTTCAGGTCTTGACTCACTCATCTGGGCACCTCCCGGTGTAGGAATGAAAACAGGTTTCTGGCGACGAACGCGCCaa                                <  1:113271/149‑1 (MQ=255)
                                                                                                                cgAGCCCACGCCAATGCGCGCCAGCACTTTCTGCAGCTTTTCGCCGGATGGCGGGGTGATTTCGGTTTCTTTCAGGTCTTGCTCACTCATCTGGGCACCTCCCGGTGTAGGAATGAAAACAGGTTTCTGGCGACGAACGCGCCAAAAgg                            <  4:369011/149‑1 (MQ=255)
                                                                                                                    cccACGCCAATGCGCGCCAGCACTTTCTGCAGCTTTTCGCCGGATGGCGGGGTGATTTCGGTTTCTTTCAGGTCTTGCTCACTCATCTGGGCACCTCCCGgt                                                                       <  4:166528/102‑1 (MQ=255)
                                                                                                                    cccACGCCAATGCGCGCCAGCACTTTCTGCAGCTTTTCGCCGGATGGCGGGGTGATTTCGGTTTCTTTCAGGTCTTGCTCACTCATCTGGGCACCTCCCGgt                                                                       >  3:166528/1‑102 (MQ=255)
                                                                                                                     ccACGCCAATGCGCGCCAGCACTTTCTGCAGCTTTTCGGCGGATGGCGGGGTGATTTCGGTTTCTTTCAGGTGTTGCTCACTCATCTCGGCACCTCCCGGTGTAGGAATGAAAACAGGTTTCTGGCGCCGAACGCGCCAAAAGGCcgcg                       <  4:364530/149‑1 (MQ=255)
                                                                                                                     ccACGCCAATGCGCGCCAGCACTTTCTGCAGCTTTTCGCCGGATGGCGGGGTGATTTCGGTTTCTTTCAGGTCTTGCTCACTCATCTGGGCACCTCCCGGTGTAGGAATGAAAACAGGTTTCTGGCGACGAACGCGCCAAAAGGCcgcg                       <  5:189993/149‑1 (MQ=255)
                                                                                                                                 gcgcCAGCACTTTCTGCAGCTTTTCGCCGGATGGCGGGGTGATTTCGGTTTCTTTCAGGTCTTGCTCACTCATCTGGGCACCTCCCGGTGTa                                                                    >  6:153929/1‑92 (MQ=255)
                                                                                                                                 gcgcCAGCACTTTCTGCAGCTTTTCGCCGGATGGCGGGGTGATTTCGGTTTCTTTCAGGTCTTGCTCACTCATCTGGGCACCTCCCGGTGTa                                                                    <  5:153929/92‑1 (MQ=255)
                                                                                                                                          cTTTCTGCAGCTTTTCGCCGGATGGCGGGGTGATTTCGGTTTCTTTCAGGTCTTGCTCACTCATCTGGGCACCTCCCGGTGTAGGAATGAAAACAGGTTTCTGGCGACGAACGCGCCAAAAGGCCGCGCATCATACGCGGTTCATCAgg  <  2:192080/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
TGCCGTCGACGGCGATGGCGTCGTGCAGGTCGACGCGCTGGCCAAGGGTGGCCTCGACGCCGTTGACCTTGATGCGACCCTGGCTGATCCAGGCCTCGACGTCACGCCGCGAGCCCACGCCAATGCGCGCCAGCACTTTCTGCAGCTTTTCGCCGGATGGCGGGGTGATTTCGGTTTCTTTCAGGTCTTGCTCACTCATCTGGGCACCTCCCGGTGTAGGAATGAAAACAGGTTTCTGGCGACGAACGCGCCAAAAGGCCGCGCATCATACGCGGTTCATCAGG  >  NC_002947/5112213‑5112496

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