Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,029,397 T→A 6.7% K191N (AAA→AAT)  rlmI ← 23S rRNA m(5)C1962 methyltransferase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,029,3970TA 6.7% 182.3 / 3.7 60K191N (AAA→AAT) rlmI23S rRNA m(5)C1962 methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base T (28/28);  new base A (2/2);  total (30/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.88e-01

CCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGTTTTGTGTCCGTGCTGAATATCCACCAGCAGTTTCATTCCGTGTTCTTCAATCGGCAGCAGGGCAGGTGGCAACTCGCCGGTGACGGGGCCCTGGGTCAGCTCCATTCCTTCTTTTTTACGTACCGCGACGTCGCTGCGATCGTAAATCGAACATTCCGGGTACAGCGTTTGCAGGGCA  >  NC_000913/1029249‑1029544
                                                                                                                                                    |                                                                                                                                                   
ccATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGTTTTGTGTCCGTGCTGAATATCCACCAGCAGt                                                                                                                                                     >  5:166628/1‑149 (MQ=255)
     aGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGTTTTGTGTCCGTGCTGAATATCCACCAGCAGTTTCa                                                                                                                                                 <  6:59488/148‑1 (MQ=255)
     aGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGTTTTGTGTCCGTGCTGAATATCCACCAGCAGTTTCAt                                                                                                                                                <  7:259136/149‑1 (MQ=255)
      gTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGTTTTGTGTCCGTGCTGAATATCCCACAGCAGTTTCAtt                                                                                                                                               >  2:151835/1‑149 (MQ=255)
         ccGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGTTTTGTGTCCGTGCTGAATATCCACCAGCAGTTTCATTCCg                                                                                                                                            >  4:151907/1‑149 (MQ=255)
          cGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGTTTTGTGTCCGTGCTGAATATCCACCAGCAGTTTCATTCCgt                                                                                                                                           <  6:166628/149‑1 (MQ=255)
               cGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGTTTTGTGTCCGTGCTGAATATCCACCAGCAGTTTCATTCCGTGttct                                                                                                                                      <  4:182942/149‑1 (MQ=255)
                 gCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGTTTTGTGTCCGTGCTGAATATCCACCAGCAGTTTCATTCCGTGttctt                                                                                                                                     >  1:318722/1‑148 (MQ=255)
                         accGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGTTTTGTGTCCGTGCTGAATATCCACCAGCAGTTTCATTCCGTGTTCTTCAATCGgca                                                                                                                            <  7:276363/149‑1 (MQ=255)
                              tataGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGTTTTGTGTCCGTGCTGAATATCCACCAGCAGTTTCATTCCGTGTTCTTCAATCGGCAgcagg                                                                                                                       <  5:59705/149‑1 (MQ=255)
                              tataGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGTTTTGTGTCCGTGCTGAATATCCACCAGCAGTTTCATTCCGTGTTCTTCAATCGGCAgcagg                                                                                                                       >  8:241014/1‑149 (MQ=255)
                                  ggAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGTTTTGTGTCCGTGCTGAATATCCACCAGCAGTTTCATTCCGTGTTCTTCAATCGGCAgcagggcag                                                                                                                   >  2:232487/1‑149 (MQ=255)
                                  ggAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGTTTTGTGTCCGTGCTGAATATCCACCAGCAGTTTCATTCCGTGTTCTTCAATCGGCAgcagggca                                                                                                                    <  6:298421/148‑1 (MQ=255)
                                   gagaAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGTTTTGTGTCCGTGCTGAATATCCACCAGCAGTTTCATTCCGTGTTCTTCAATCGGCAgcagggcagg                                                                                                                  >  1:44075/1‑149 (MQ=255)
                                      aaaCAGTTCAGCACACGTTTATTTTCATCGTAGAGGCGGTTAGCCAGGCGTCTATCACGCTGGTCCAGGTAGTAGCCCGTTTTGTGTCCGTGCTGAATATCCACCAGCAGTTTCATTCCGTGTTCTTCAATCGGCAGCAGGGCAGGTgg                                                                                                               <  1:246198/149‑1 (MQ=255)
                                               aGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGTTTTGTGTCCGTGCTGAATATCCCCCAGCAGTTTCATTCCGTGTTCTTCAATCGGCCGCAGGGCAGGTGGCAACTCGcc                                                                                                      >  2:207291/1‑149 (MQ=255)
                                               aGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGTTTTGTGTCCGTGCTGAATATCCACCAGCAGTTTCATTCCGTGTTCTTCAAt                                                                                                                                 <  7:4190/122‑1 (MQ=255)
                                               aGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGTTTTGTGTCCGTGCTGAATATCCACCAGCAGTTTCATTCCGTGTTCTTCAAt                                                                                                                                 >  8:4190/1‑122 (MQ=255)
                                                gCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGTTTTGTGGCCGTGCTGAATATCCACCAGCAGTTTCATTCCGTGTTCTTCAATCGGCAGCAGGGCAGGTGGCAACTCGCCg                                                                                                     <  1:151835/149‑1 (MQ=255)
                                                    acGTTTATTTTCAACGTAGCGGC‑GGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGTTTTGTGTCCGTGCTGAATATCCACCAGCAGTTTCATTCCGTGTTCTTCAATCGgcagca                                                                                                                         <  8:303798/124‑1 (MQ=255)
                                                    acGTTTATTTTCAACGTAGCGGC‑GGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGTTTTGTGTCCGTGCTGAATATCCACCAGCAGTTTCATTCCGTGTTCTTCAATCGgcagca                                                                                                                         >  7:303798/1‑124 (MQ=255)
                                                      gTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGTTTTGTGTCCGTGCTGAATATCCACCAGCAGTTTCATTCCGTGTTCTTCAATCGGCAGCAGGGCAGGTGGCAACTCGCCGGTGACg                                                                                               >  2:37896/1‑149 (MQ=255)
                                                        ttATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGTTTTGTGTCCGTGCTGAATATCCACCAGCAGATTCATTCCGTGTTCTTCAATCGGCAGCAGGGCAGGt                                                                                                                 >  4:307538/1‑129 (MQ=255)
                                                        ttATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGTTTTGTGTCCGTGCTGAATATCCACCAGCAGATTCATTCCGTGTTCTTCAATCGGCAGCAGGGCAGGt                                                                                                                 <  3:307538/129‑1 (MQ=255)
                                                          aTTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGTTTTGTGTCCGTGCTGAATATCCACCAGCAGTTTCATTCCGTGTTCTTCAATCGGCAgcagg                                                                                                                       <  1:97374/121‑1 (MQ=255)
                                                          aTTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGTTTTGTGTCCGTGCTGAATATCCACCAGCAGTTTCATTCCGTGTTCTTCAATCGGCAgcagg                                                                                                                       >  2:97374/1‑121 (MQ=255)
                                                              tCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGTTTTGTGTCCGTGCTGAATATCCACCAGCAGTTTCATTCCGTGTTCTTCAATCGGCAGCAGGGCAGGTGGCAACTCGCCGGTGACGGGGCCCt                                                                                        <  8:114951/148‑1 (MQ=255)
                                                                aaCGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGTTTTGTGTCCGTGCTGAATATCCACCAGCAGTTTCATTCCGTGTTCTTCAATCGGCAGCAGGGCAGGTGGCAACTCGCCGGTGACGGGGCCCTgg                                                                                      <  7:138224/148‑1 (MQ=255)
                                                                aaCGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGTTTTGTGTCCGTGCTGAATATCCACCAGCAGTTTCATTCCGTGTTCTTCAATCGGCAGCAGGGCAGGTGGCAACTCGCCGGTGACGGGGCCCTgg                                                                                      <  3:151907/148‑1 (MQ=255)
                                                                     aGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGGTTTGTGTCCGTGCGGAATATCCACCAGCAGTTTCATTCCGTGGTCGTCAATCGGCCGCAGGGCAGGTGGCTACTCGCCGGTGACGGGGCCCTGGGTCAGc                                                                                >  5:110510/1‑149 (MQ=255)
                                                                      gcggcgGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGTTTTGTGTCCGTGCTGAATATCCACCAGCAGTTTCATTCCGTGTTCTTCAATCGGCAGCAGGGCAGGTGGCAACTCGCCGGTGACGGGGCCCTGGGTCAGc                                                                                <  1:276259/148‑1 (MQ=255)
                                                                                  cAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGTTTTGTGTCCGTGCTGAATATCCACCAGCAGTTTCATTCCGTGTTCTTCAATCGGCAGCAGGGCAGGTGGCAACTCGCCGGTGACGGGGCCCTTGGTCAGCTCCATTCcttctt                                                                   >  4:38148/1‑149 (MQ=255)
                                                                                   aGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGTTTTGTGTCCGTGCTGAATATCCACCAGCAGTTTCATTCCGTGTTCTTCAATCGGCAGCAGGGCAGGTGGCAACTCGCCGGTGACGGGGCCCTGGGTCAGCTCCATTCCTTCttt                                                                  >  2:153008/1‑149 (MQ=255)
                                                                                       ggcTATCACGCTGGTCCAGGTAGTAGCCCGTTTTGTGTCCGTGCTGAATATCCACCAGCAGTTTCATTCCGTGTTCTTCAATCGGCAGCAGGGCAGGTGGCAACTCGCCGGTGACGGGGCCCTGGGTCAGCTCCATTCCTTCTTTTTTa                                                              >  8:79944/1‑149 (MQ=255)
                                                                                       ggcTATCACGCTGGTCCAGGTAGTAGCCCGTTTTGTGTCCGTGCTGAATATCCACCAGCAGTTTCATTCCGTGTTCTTCAATCGGCAGCAGGGCAGGTGGCAACTCGCCGGTGACGGGGCCCTGGGTCAGCTCCATTCCTTCTTTTTTa                                                              >  1:356752/1‑149 (MQ=255)
                                                                                          tATCACGCTGGTCCAGGTAGTAGCCCGTTTTGTGTCCGTGCTGAATATCCACCCGCAGTTTCATTCCCTGTTCTTCAAACGCCAGCAGGGCAGGTGGCAACTCGCCGGTGCCGGGGCCCTGGGTCAGCTCCAATCCTTCTTTTTTACg                                                            >  2:346269/1‑148 (MQ=255)
                                                                                              aCGCTGGTCCAGGTAGTAGCCCGTTTTGTGTCCGTGCTGAATATCCACCAGCAGTTTCATTCCGTGTTCTTCAATCGGCAGCAGGGCAGGTGGCAACTCGCCGGTGACGGGGCCCTGGGTCAGCTCCa                                                                            <  7:363053/128‑1 (MQ=255)
                                                                                              aCGCTGGTCCAGGTAGTAGCCCGTTTTGTGTCCGTGCTGAATATCCACCAGCAGTTTCATTCCGTGTTCTTCAATCGGCAGCAGGGCAGGTGGCAACTCGCCGGTGACGGGGCCCTGGGTCAGCTCCa                                                                            >  8:363053/1‑128 (MQ=255)
                                                                                              aCGCTGGTCCAGGTAGTAGCCCGTTTTGTGTCCGTGCTGAATATCCACCAGCAGTTTCATTCCGTGTTCTTCAATCGGCAGCAGGGCAGGTGGCAACTCGCCGGTGACGGGGCCCTGGGTCAGCTCCATTCCTTCTTTTTTACGTACcg                                                       >  6:270707/1‑149 (MQ=255)
                                                                                                gCTGGTCCAGGTAGTAGCCCGTTTTGTGTCCGTGCTGAATATCCACCAGCAGTTTCATTCCGTGTTCTTCAATCGGCAGCAGGGCAGGTGGCAACTCGCCGGTGACGGGGCCCTGGGTCAGCTCCATTCCTTCTTTTTTACGTACcgc                                                      >  4:221748/1‑148 (MQ=255)
                                                                                                gCTGGTCCAGGTAGTAGCCCGTTTTGTGTCCGTGCTGAATATCCACCAGCAGTTTCATTCCGTGTTCTTCAATCGGCAGCAGGGCAGGTGGCAACTCGCCGGTGACGGGGCCCTGGGTCAGCTCCATTCCTTCTTTTTTACGTACcgc                                                      >  5:35472/1‑148 (MQ=255)
                                                                                                          gtagtaGCCCGTTTTGTGTCCGTGCTGAATATCCACCAGCAGTTTCATTCCGTGTTCTTCAATCGGCAGCAGGGCAGGTGGCAACTCGCCGGTGACGGGGCCCTGGGTCAGCTCCATTCCTTCTTTTTTACGTACCGCGACGTCGCTGc                                           >  5:152201/1‑149 (MQ=255)
                                                                                                          gtagtaGCCCGTTTTGTGTCCGTGCTGAATATCCACCAGCAGTTTCATTCCGTGTTCTTCAATCGGCAGCAGGGCAGGTGGCAACTCGCCGGTGACGGGGCCCTGGGTCAGCTCCATTCCTTCTTTTTTACGTACCGCGACGTCGCTGc                                           <  2:318722/149‑1 (MQ=255)
                                                                                                           tagtagCCCGTTTTGTGTCCGTGCTGAATATCCACCAGCAGTTTCATTCCGTGTTCTTCAATCGGCAGCAGGGCAGGTGGAAACTCGCCGGTGACGTGGCCCTGGGTCAGCTCCATTCCTTCTTTTTTACGTAACGCGACGTCGCTGc                                           >  6:184149/1‑148 (MQ=255)
                                                                                                                gCCCGTTTTGTGTCCGTGCTGAATATCCACCAGCAGTTTCATTCCGTGTTCTTCAATCGGCAGCAGGGCAGGTGGCAACTCGCCGGTGACGGGGCCCTGGGTCAGCTCCATTCCTTCTTTTTTACGTACCGCGACGTCGCTGCGATCg                                      >  3:360764/1‑148 (MQ=255)
                                                                                                                gCCCGTTTTGTGTCCGTGCTGAATATCCACCAGCAGTTTCATTCCGTGTTCTTCAATCGGCAGCAGGGCAGGTGGCAACTCGCCGGTGACGGGGCCCTGGGTCAGCTCCATTCCTTCTTTTTTACGTACCGCGACGTCGCTGCGATCg                                      <  7:79944/148‑1 (MQ=255)
                                                                                                                 cccGTTTTGTGTCCGTGCTGAATATCCACCAGCAGTTTCATTCCGTGTTCTCCAATCGGCAGCAGGGCAGGTGGCAACTCGCCGGTGACGGGGCCCTGGGTCAGCTCCATTCCTTCTTTTTTACGTACCGCGACGTCGCTGCGATCGt                                     >  3:160633/1‑148 (MQ=255)
                                                                                                                    gTTTTGTGTCCGTGCTGAATATCCACCAGCAGTTTCATTCCGTGTTCTTCAATCGGCAGCAGGGCAGGTGGCAACTCGCCGGTGACGGGGCCCTGGGTCAGCTCCATTCCTTCTTTTTTACGTACCGCGACGTCGCTGCGATCGTAAAt                                 <  2:356752/149‑1 (MQ=255)
                                                                                                                    gTTTTGTGTCCGTGCTGAATATCCACCAGCAGTTTCATTCCGTGTTCTTCAATCGGCAGCAGGGCAGGTGGCAACTCGCCGGTGACGGGGCCCTGGGTCAGCTCCATTCCTTCTTTTTTACGTACCGCGACGTCGCTGCGATCGTAAAt                                 <  1:37896/149‑1 (MQ=255)
                                                                                                                           gtCCGTGCTGAATATCCACCAGCAGATTCATTCCGTGTTCTTCAATCGGCAGCAGGGCAGGTGGCAACTCGCCGGTGACGGGGCCCTGGGTCAGCTCCATTCCTTCTTTTTTACGTACCGCGa                                                    <  3:69969/123‑1 (MQ=255)
                                                                                                                           gtCCGTGCTGAATATCCACCAGCAGATTCATTCCGTGTTCTTCAATCGGCAGCAGGGCAGGTGGCAACTCGCCGGTGACGGGGCCCTGGGTCAGCTCCATTCCTTCTTTTTTACGTACCGCGa                                                    >  4:69969/1‑123 (MQ=255)
                                                                                                                            tCCGTGCTGAATATCCACCAGCAGTTTCATTCCGTGTTCTTCAATCGGCAGCAGGGCAGGTGGCAACTCGCCGGTGACGGGGCACTGGGTCAGCTCCATTCCTTCTTTTTTACGTACCGCGACGTCGCTGCGATCGTAAATCGAACAtt                         >  5:5778/1‑149 (MQ=255)
                                                                                                                             ccGTGCTGAATATCCACCAGCAGTTTCATTCCGTGTTCTTCAATCGGCAGCAGGGCAGGTGGCAACTCGCCGGTGACGGGGCCCTGGGTCAGCTCCATTCCTTCTTTTTTACGTACCGCGACGTCGCTGCGATCGTAAATCGAACATTc                        <  2:44075/149‑1 (MQ=255)
                                                                                                                              cGTGCTGAATATCCACCAGCAGTTTCATTCCGTGTTCTTCAATCGGCAGCAGGGCAGGTGGCAACTCGCCGGTGACGGGGCCCTGGGTCAGCTCCATTCCTTCTTTTTTACGTACCGCGACGTCGCTGCGATCGTAAATCGAACATTcc                       <  6:64729/149‑1 (MQ=255)
                                                                                                                                  cTGAATATCCACCAGCAGTTTCATTCCGTGTTCTTCAATCGGCAGCAGGGCAGGTGGCAACTCGCCGGTGACGGGGCCCTGGGTCAGCTCCATTCCTTCTTTTTTACGTACCGCGACGTCGCTGCGATCGTAAATCGAACATTCCggg                    <  6:35472/148‑1 (MQ=255)
                                                                                                                                      atatCCACCAGCAGTTTCATTCCGTGTTCTTCAATCGGCAGCAGGGCAGGTGGCAACTCGCCGGTGACGGGGCCCTGGGTCAGCTCCATTCCTTCTTTTTTACGTACCGCGACGTCGCTGCGATCGTAAATCGAACATTCCGGGTACAg               >  4:235084/1‑149 (MQ=255)
                                                                                                                                          ccaccaGCAGTTTCATTCCGTGTTCTTCAATCGGCAGCAGGGCAGGTGGCAACTCGCCGGTGACGGGGCCCTGGGTCAGCTCCATTCCTTCTTTTTTACGTACCGCGACGTCGCTGCGATCGTAAATCGAACATTCCGGGTACAGCGtt           >  1:70675/1‑149 (MQ=255)
                                                                                                                                           caccaGCAGTTTCATTCCGTGTTCTTCAATCGGCAGCAGGGCAGGTGGCAACTCGCCGGTGACGGGGCCCTGGGTCAGCTCCATTCCTTCTTTTTTACGTACCGCGACGTCGCTGCGATCGTAAATCGAACATTCCGGGTACAGCGtt           <  1:232487/148‑1 (MQ=255)
                                                                                                                                                 cagTTTCATTCCGTGTTCTTCAATCGGCAGCAGGGCAGGTGGCAACTCGCCGGTGACGGGGCCCTGGGTCAGCTCCATTCCTTCTTTTTTACGTACCGCGACGTCGCTGCGATCGTAAATCGAACATTCCGGGTACAGCGTTTGCAggg    <  2:70675/149‑1 (MQ=255)
                                                                                                                                                   gTTTCATTCCGTGTTCTTCAATCGGCAGCAGGGCAGGTGGCAACTCGCCGGTGACGGGGCCCTGGGTCAGCTCCATTCCTTCTTTTTTACGTACCGCGACGTCGCTGCGATCGTAAATCGAACATTCCGGGTACAGCGTTTGCAGGGCa  <  3:38148/149‑1 (MQ=255)
                                                                                                                                                   gTTTCATTCCGTGTTCTTCAATCGGCAGCAGGGCAGGTGGCAACTCGCCGGTGACGGGGCCCTGGGTCAGCTCCATTCCTTCTTTTTTACGTACCGCGACGTCGCTGCGATCGTAAATCGAACATTCCGGGTACAGCGTTTGCAGGGCa  <  1:153008/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGTTTTGTGTCCGTGCTGAATATCCACCAGCAGTTTCATTCCGTGTTCTTCAATCGGCAGCAGGGCAGGTGGCAACTCGCCGGTGACGGGGCCCTGGGTCAGCTCCATTCCTTCTTTTTTACGTACCGCGACGTCGCTGCGATCGTAAATCGAACATTCCGGGTACAGCGTTTGCAGGGCA  >  NC_000913/1029249‑1029544

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