Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,029,516 A→T 7.2% S152T (TCG→ACG)  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,5160AT 7.2% 200.6 / 4.2 64S152T (TCG→ACG) rlmI23S rRNA m(5)C1962 methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base A (31/28);  new base T (2/3);  total (33/31)
Fisher's exact test for biased strand distribution p-value = 6.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCGTGCTGAATATCCACCAGCAGTTTCATTCCGTGTTCTTCAATCGGCAGCAGGGCAGGTGGCAACTCGCCGGTGACGGGGCCCTGGGTCAGCTCCATTCCTTCTTTTTTACGTACCGCGACGTCGCTGCGATCGTAAATCGAACATTCCGGGTACAGCGTTTGCAGGGCACTAATTAATGCCGCGCGCTGATATTCTGCCCCAGCACTCAGCAGTTGCAGCACCAGAAAATTACCGAAACGATCGATAGTAATACCCGGCAGGCCATCAGATTCTCCGGCGATTAAACG  >  NC_000913/1029373‑1029663
                                                                                                                                               |                                                                                                                                                   
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                                                                                                                         <  1:153008/149‑1 (MQ=255)
                       gTTTCATTCCGTGTTCTTCAATCGGCAGCAGGGCAGGTGGCAACTCGCCGGTGACGGGGCCCTGGGTCAGCTCCATTCCTTCTTTTTTACGTACCGCGACGTCGCTGCGATCGTAAATCGAACATTCCGGGTACAGCGTTTGCAGGGCa                                                                                                                         <  3:38148/149‑1 (MQ=255)
                         ttCATTCCGTGTTCTTCAATCGGCAGCAGGGCAGGTGGCAACTCGCCGGTGACGGGGCCCTGGGTCAGATCCATTCCTTCTTTTTTACGTACCGCGACGTCGCTGCGATCGTAATTCGTACATTCCGGGTACAGCGTTTGCAGTGCAc                                                                                                                        <  4:62028/148‑1 (MQ=255)
                           cATTCCGTGTTCTTCAATCGGCAGCAGGGCAGGTGGCAACTCGCCGGTGACGGGGCCCTGGGTCAGCTCCATTCCTTCTTTTTTACGTACCGCGACGTCGCTGCGATCGTAAATCGAACATTCCGGGTACAGCGTTTGCAGGGCACtaa                                                                                                                     >  4:331264/1‑149 (MQ=255)
                            aTTCCGTGTTCTTCAATCGGCAGCAGGGCAGGTGGCAACTCGCCGGTGACGGGGCCCTGGGTCAGCTCCATTCCTTCTTTTTTACGTACCGCGACGTCGCTGCGATCGTAAATCGAACATTCCGGGTACAGCGTTTGCAGGGCACtaa                                                                                                                     <  5:270707/148‑1 (MQ=255)
                             ttCCGTGTTCTTCAATCGGCAGCAGGGCAGGTGGCAACTCGCCGGTGACGGGGCCCTGGGTCAGCTCCATTCCTTCTTTTTTACGTACCGCGACGTCGCTGCGATCGTAAATCGAACATTCCGGGTACAGCGTTTGCAGGGCACtaat                                                                                                                    >  1:51891/1‑148 (MQ=255)
                               ccGTGTTCTTCAATCGGCAGCAGGGCAGGTGGCAACTCGCCGGTGACGGGGCCCTGGGTCAGCTCCATTCCTTCTTTTTTACGTACCGCGACGTCGCTGCGATCGTAAATCGAACATTCCGGGTACAGCGTTTGCAGGGCACtaattaa                                                                                                                 >  4:108386/1‑149 (MQ=255)
                                    ttcttcAATCGGCAGCAGGGCAGGTGGCAACTCGCCGGTGACGGGGCCCTGGGTCAGCTCCATTCCTTCTTTTTTACGTACCGCGACGTCGCTGCGATCGTAAATCGAACATTCCGGGTACAGCGTTTGCAGGGCACTAATTAATGCcg                                                                                                            <  1:207291/149‑1 (MQ=255)
                                          aaTCGGCAGCAGGGCAGGTGGCAAATCGCCGGTGCCTTGGCCATGGGTCAGCTCCATTCCTTCTTTTTTAATTACCGCGACGTCGCTGCGATCGTAAATCGAACATTACGGGTAAAGCGTTTGAAGGGCACTAATTAATGCCGCGAGCt                                                                                                      >  8:280442/1‑149 (MQ=255)
                                          aaTCGGCAGCAGGGCAGGTGCCAACTCGCCGGTGAAGGGGCCCTGGGACAGCTCCAATCCTACTTTTTTACGTACCGCGACGTCGCTGCGATCGTAAATCGAACATTCCGGGTACAGCGTTTGCAGGGCACTAAATAATGCCGCGCGCt                                                                                                      >  4:154863/1‑149 (MQ=255)
                                           aTCGGCAGCAGGGCAGGTGGCAACTCGCCGGTGACGGGGCCCTGGGTCAGCTCCATTCCTTCTTTTTTACGTACCGCGACGTCGCTGCGATCGTAAATCGAACATTCCGGGTACAGCGTTTGCAGGGCACTAATTAATGCCGCGCGCTg                                                                                                     <  3:331264/149‑1 (MQ=255)
                                            tCGGCAGCAGGGCAGGTGGCAACTCGCCGGTGACGGGGCCCTGGGTCAGCTCCATTCCTTCTTTTTTACGTACCGCGACGTCGCTGCGATCGTAAATCGAACATTCCGGGTACAGCGTTTGCAGGGCACTAATTAATGCCGCGCGCTGa                                                                                                    <  3:235084/149‑1 (MQ=255)
                                              ggCAGCAGGGCAGGTGGCAACTCGCCGGTGACGGGGCCCTGGGGCAGCTCCATTCCTTCTTTTTTACGTAACGCTACCTCGCTGCTATCGTAAATCGAAAATTCCGGGTACAGCGTTTTCAGGGCACTAATAAATGCCGCGCGATGata                                                                                                  >  5:144253/1‑149 (MQ=255)
                                                ccgcagggcaggTGGCTCCGCGCCGGGGACGGGGCCCTGGGTCAGCTCCCTGCCGTCTTTTTTACGTCCCGCGACCGCGCGGCGAGCGTACATCGAACATTCCGGGTACAGCGTGTGCAGGGCACTAATTAATGCCGCGCGCTGATAtt                                                                                                <  6:110510/147‑1 (MQ=255)
                                                   cagggcaggTGGCAACTCGCCGGTGACGGGGCCCTGGGTCAGCTCCATTCCTTCTTTTTTACGTACCGCGACGTCGCTGCGATCGTAAATCGAACATTCCGGGTACAGCGTTTGCAGGGCACTAATTAATGCCGCGCGCTGATATTCTg                                                                                             >  2:15032/1‑149 (MQ=255)
                                                        caggTGGCAACTCGCCGGTGACGGGGCCCTGGGTCAGCTCCATTCCTTCTTTTTTACGTACCGCGACGTCGCTGCGATCGTAAATCGAACATTCCGGGTACAGCGTTTGCAGGGCACTAATTAATGCCGCGCGCTGATATTCTGCCCCa                                                                                        <  3:108386/149‑1 (MQ=255)
                                                           gTGGCTACTCGCCGGTGACGGGGCCCTGGGTCAGCTCCATTCCTTCTTTTTTACGTACCGCGACGTCGCTGCGATCGTAAATCGaa                                                                                                                                                    <  1:380594/86‑1 (MQ=255)
                                                           gTGGCAACTCGCCGGTGACGGGGCCCTGGGTCAGCTCCATTCCTTCTTTTTTACGTACCGAGACGTCGCTGCGAACCTAAATCGaa                                                                                                                                                    >  2:380594/1‑86 (MQ=255)
                                                              gCAACTCGCCGGTGACGGGGCCCTGGGTCAGCTCCATTCCTTCTTTTTTACGTACCGCGACGTCGCTGCGATCGTAAATCGAACATTCCGGGTACAGCGTTTGCAGGGCACTAATTAATGCCGCGCGCTGATATTCTGCCCCAGCACt                                                                                   >  6:75590/1‑148 (MQ=255)
                                                               cAACTCGCCGGTGACGGGGCCCTGGGTCAGCTCCATTCCTTCTTTTTTACGTACCGCGACGTCGCTGCGATCGTAAATCGAACATTCCGGGTACAGCGTTTGCAGGGCACTAATTAATGCCGCGCGCTGATATTCTGCCCCAGCACTca                                                                                 <  2:51891/149‑1 (MQ=255)
                                                                aaCTCGCCGGTGACGGGGCCCTGGGTCAGCTCCATTCCTTCTTTTTTACGTACCGCGACGTCGCTGCGATCGTAAATCGAACATTCCGGGTACAGCGTTTGCAGGGCACTAATTAATGCCGCGCGCTGATATTCTGCCCCAGCACTcag                                                                                <  5:75590/149‑1 (MQ=255)
                                                                aaCTCGCCGGTGACGGGGCCCTGGGTCAGCTCCATTCCTTCTTTTTTACGTACCGCGACGTCGCTGCGATCGTAAATCGAACATTCCGGGTACAGCGTTTGCAGGGCACTAATTAATGCCGCGCGCTGATATTCTGCCCCAGCACTcag                                                                                <  5:184149/149‑1 (MQ=255)
                                                                     gCCGGTGACGGGGCCCTGGGTCAGCTCCATTCCTTCTTTTTTACGTACCGCGACGTCGCTGCGATCGTAAATCGAACATTCCGGGTACAGCGTTTGCAGGGCACTAATTAATGCCGCGCGCTGATATTCTGCCCCAGCACTCAGCAGtt                                                                           >  7:43693/1‑149 (MQ=255)
                                                                       cGGCGACGGGGCCCTGGGTCAGCTCCCTTCCTTCTTTTTTACGGACCGCGACGGCGCTGCGATCGTAAATCGAACCTTCCGGGGACAGCCTTTGCAGGGCACTAATTAATGCCGCGCGCTGATATTCTGCCCCAGCACTCAGCAGTTgc                                                                         <  6:5778/149‑1 (MQ=255)
                                                                         gTGACGGGGCCCTGGGTCAGCTCCATTCCTTCTTTTTTACGTACCGCGACGTCGCTGCGATCGTAAATCGAACATTCCGGGTACAGCGTTTGCAGGGCACTAATTAATGTCGCGCGCTGATATTCTGCCCCAGCACTCAGCAGTTgcag                                                                       >  3:372703/1‑149 (MQ=255)
                                                                            aCGGGGCCCTGGGTCAGCTCCATTCCTTCTTTTTTACGTACCGCGACGGCGCTGCGATCGTAAATCGAACATTCCGGGTACAGCGTTTGCAGGGCACTAATTAATGCCGCGCGCTGATATTCTGCCCCAGCACTCAGCAGTTGCAGCAc                                                                    <  7:241014/149‑1 (MQ=255)
                                                                              ggggCCCTGGGTCAGCTCCATTCCTTCTTTTTTACGTACCGCGACGTCGCTGCGATCGTAAATCGAACATTCCGGGTACAGCGTTTGCAGGGCACTAATTAATGCCGCGCGCTGATATTCTGCCCCAGCACTCAGCAGTTGCAGCAc                                                                    >  8:397488/1‑147 (MQ=255)
                                                                                ggCCCTGGGTCAGCTCCATTCCTTCTTTTTTACGTACCGCGACGTCGCTGCGATCGTAAATCGAACATTCCGGGTACAGCGTTTGCAGGGCACTAATTAATGCCGCGCGCTGATATTCTGCCCCAGCACTCAGCAGTTGCAGCACCAGa                                                                <  4:160633/149‑1 (MQ=255)
                                                                                    cTGGGTCAGCTCCATTCCTTCTTTTTTACGTACCGCGACGTCGCTGCGATCGTAAATCGAACATTCCGGGTACAGCGTTTGCAGGGCACTAATTAATGCCGCGCGCTGATATTCTGCCCCAGCACTCAGCAGTTGCAGCACCAGAAAAt                                                            <  7:397488/149‑1 (MQ=255)
                                                                                       ggTGAGCTCCATTGCTTCTTTTTTACGTACCGCGACGTCGCTGGGATCGTAAATCGAACATTCCGGGTACAGCGTTTGCAGGGCACTAATTAATGCAGCGCGCTGATATTCTGCCCCAGCAGTCAGCTGTTGCAGCACCAGAAAATTa                                                          <  1:346269/148‑1 (MQ=255)
                                                                                       ggTCAGCTCCATTCCTTCTTTTTTACGTACCGCGACGTCGCTGCGATCGTAAATCGAACATTCCGGGTACAGCGTTTGCAGGGCACTAATTAATGCCGCGCGCTGATATTCTGCCCCAGCACTCAGCAGTTGCAGCACCAGAAAATTa                                                          <  3:221748/148‑1 (MQ=255)
                                                                                           aGCTCCATTCCTTCTTATTTACGTACCGCGACGTCGCTGCGATCGTAAATCGTACATTCCGGGTACAGCGTTTGCAGGGCACTAATTAATGCCGCGCGCTGATATTCTGCCCCAGCACTCAGCAGTTGCAGCACCAGAAAATTACCGa                                                      >  5:269992/1‑148 (MQ=255)
                                                                                               ccATTCCTTCTTTTTTACGTACCGCGACGTCGCTGCGATCGTAAATCGAACATTCCGGGTACAGCGTTTGCAGGGCACTAATTAATGCCGCGCGCTGATATTCTGCCCCAGCACTCAGCAGTTGCAGCACCAGAAAATTACCGAAAcga                                                 >  5:315178/1‑149 (MQ=255)
                                                                                                 aTTCCTTCTTTTTTACGTACCGCGACGTCGCTGCGATCGTAAATCGAACATTCCGGGTACAGCGTTTGCAGGGCACTAATTAATGCCGCGCGCTGATATTCTGCCCCAGCACTCAGCAGTTGCAGCACCAGAAAATTACCGAAAcgatc                                               >  2:38734/1‑149 (MQ=255)
                                                                                                  ttCCTTCTTTTTTACGTACCGCGACGTCGCTGCGATCGTAAATCGAACATTCCGGGTACAGCGTTTGCAGGGCACTAATTAATGCCGCGCGCTGATATTCTGCCCCAGCACTCAGCAGTTGCAGCACCAGAAAATTACCGAAAcgatc                                               >  2:51735/1‑148 (MQ=255)
                                                                                                      ttcttTTTTACGTACCGCGACGTCGCTGCGATCGTAAATCGAACATTCCGGGTACAGCGTTTGCAGGGCACTAATTAATGCCGCGCGCTGATATTCTGCCCCAGCACTCAGCAGTTGCAGCACCAGAAAATTACCGAAACGATCGATAg                                          >  8:381869/1‑149 (MQ=255)
                                                                                                          tttttACGTACCGCGACGTCGCTGCGATCGTAAATCGAACATTCCGGGTACAGCGTTTGCAGGGCACTAATTAATGCCGCGCGCTGATATTCTGCCCCAGCACTCAGCAGTTGCAGCACCAGAAAATTACCGAAACGATCGATAGTaa                                       >  4:169919/1‑148 (MQ=255)
                                                                                                          tttttACGTACCGCGACGTCGCTGCGATCGTAAATCGAACATTCCGGGTACAGCGTTTGCAGGGCACTAATTAATGCCGCGCGCTGATATTCTGCCCCAGCACTCAGCAGTTGCAGCACCAGAAAATTACCGAAACGATCGATAGTaa                                       >  8:228374/1‑148 (MQ=255)
                                                                                                                 gTACCGCGACGTCGCTGCGATCGTAAATCGAACATTCCGGGTACAGCGTTTGCAGGGCACTAATTAATGCCGCGCGCTGATATTCTGCCCCAGCACTCAGCAGTTGCAGCACCAGAAAATTACCGAAACGATCGATAGTAATACCCgg                                >  2:5291/1‑148 (MQ=255)
                                                                                                                   aCCGCGACGTCGCTGCGATCGTAAATCGAACATTCCGGGTACAGCGTTTGCAGGGCACTAATTAATGCCGCGCGCTGATATTCTGCCCCAGCACTCAGCAGTTGCAGCACCAGAAAATTACCGAAACGATCGATAGTAATACCCGGCAg                             >  8:386293/1‑149 (MQ=255)
                                                                                                                    ccGCGACGTCGCTGCGATCGTAAATCGAACATTCCGGGTACAGCGTTTGCAGGGCACTAATTAATGCCGCGCGCTGATATTCTGCCCCAGCACTCAGCAGTTGCAGCACCAGAAAATTACCGAAACGATCGATAGTAATACCCGGCAgg                            >  2:83332/1‑149 (MQ=255)
                                                                                                                     cgcgACGTCGCTGCGATCGTAAATCGAACATTCCGGGTACAGCGTTTGCAGGGCACTAATTAATGCCGCGCGCTGATATTCTGCCCCAGCACTCAGCAGTTGCAGCACCAGAAAATTACCGAAACGATCGATAGTAATACCCGGCAGGc                           >  3:81398/1‑149 (MQ=255)
                                                                                                                     cgcgACGTCGCTGCGATCGTAAATCGAACATTCCGGGTACAGCGTTTGCAGGGCACTAATTAATGCCGCGCGCTGATATTCTGCCCCAGCACTCAGCAGTTGCAGCACCAGAAAATTACCGAAACGATCGATAGTAATACCCGGCAGGc                           >  2:78434/1‑149 (MQ=255)
                                                                                                                      gcgACGTCGCTGCGATCGTAAATCGAACATTCCGGGTACAGCGTTTGCAGGGCACTAATTAATGCCGCGCGCTGATATTCTGCCCCAGCACTCAGCAGTTGCAGCACCAGAAAATTACCGAAACGATCGATAGTAATACCCGGCAGGc                           >  1:192424/1‑148 (MQ=255)
                                                                                                                       cgACGTCGCTGCGATCGTAAATCGTACATTCCGGGTACAGCGTTTGCAGGGCACTAATTAATGCCGCGCGCTGATATTCTGCCCCAGCACTCAGCAGTTGCAGCACCAGAAAATTACCGAAACGATCGATAGTAATACCCGGCAGGCCa                         >  7:187298/1‑149 (MQ=255)
                                                                                                                         aCGTCGCTGCGATCGTAAATCGAACATTCCGGGTACAGCGTTTGCAGGGCACTAATTAATGCCGCGCGCTGATATTCTGCCCCAGCACTCAGCAGTTGCAGCACCAGAAAATTACCGAAACGATCGATAGTAATACCCGGCAGGCCATc                       <  1:15032/149‑1 (MQ=255)
                                                                                                                                  cGATCGTAAATCGTACATTCCGGGTACAGCGTTTGCAGGGCACTAATTAATGCCGCGCGCTGATATTCTGCCCCAGCACTCAGCAGTTGCAGCACCAGAAAATTACCGAAACGATCGATAGTAATACCCGGCAGGCCATCAGATtctc               <  6:269992/148‑1 (MQ=255)
                                                                                                                                  cGATCGTAAATCGAACATTCCGGGTACAGCGTTTGCAGGGCACTAATTAATGCCGCGCGCTGATATTCTGCCCCAGCACTCAGCAGTTGCAGCACCAGAAAATTACCGAAACGATCGATAGTAATACCCGGCAGGCCATCAGATTCTcc              >  5:97472/1‑149 (MQ=255)
                                                                                                                                   gATCGTAAATCGAACATTCCGGGTACAGCGTTTGCAGGGCACTAATTAATGCCGAGCGCTGATATTCTGCCCCAGCACTCAGCAGTTGCAGCACCAGAAAATTACCGAAACGATCGATAGTAATACCCGGCAGGCCATCAGATTCTCCg             >  1:93781/1‑149 (MQ=255)
                                                                                                                                     tCGTAAATCGAACATTCCGGGTACAGCGTTTGCAGGGCACTAATTAATGCCGCGCGCTGATATTCTGCCCCAGCACTCAGCAGTTGCAGCACCAGAAAATTACCGAAACGATCGATAGTAATACCCGGCAGGCCATCAGATTCTCCGGc           <  7:228374/149‑1 (MQ=255)
                                                                                                                                     tCGTAAATCGAACATTCCGGGTACAGCGTTTGCAGGGCACTAATTAATGCCGCGCGCTGATATTCTGCCCCAGCACTCAGCAGTTGCAGCACCAGAAAATTACCGAAACGATCGATAGTAATACCCGGCAGGCCATCAGATTCTCCGGc           >  2:197636/1‑149 (MQ=255)
                                                                                                                                       gTAAATCGTACATTCCGGGTACAGCGTTTGCAGGGCACTAATTAATGCCGCGCGCTGATATTCTGCCCCAGCACTCAGCAGTTGCAGCACCAGAAAATTACCGAAACGATCGATAGTAATACCCGGCAGGCCATCAGATTCTCCGGCGa         <  8:187298/149‑1 (MQ=255)
                                                                                                                                       gTAAATCGAACATTCCGGGTACAGCGTTTGCAGGGCACTAATTAATGCCGCGCGCTGATATTCTGCCCCAGCACTCAGCAGTTGCAGCACCAGAAAATTACCGAAACGATCGATAGTAATACCCGGCAGGCCATCAGATTCTCCGGCGa         <  1:38734/149‑1 (MQ=255)
                                                                                                                                        tAAATCGAACATTCCGGGTACAGCGTTTGCAGGGCACTAATTAATGCCGCGCGCTGATATTCTGCCCCAGCACTCAGCAGTTGCAGCACCAGAAAATTACCGAAACGATCGATAGTAATACCCGGCAGGCCATCAGATTCTCCGGCGAt        <  7:386293/149‑1 (MQ=255)
                                                                                                                                         aaaTCGAACATTCCGGGTACAGCGTTTGCAGGGCACTAATTAATGCCGCGCGCTGATATTCTGCCCCAGCACTCAGCAGTTGCAGCACCAGAAAATTACCGAAACGGTCGATAGTAATACCCGGCAGGCCATCAGATTCTCCGGCGAt        >  3:35089/1‑148 (MQ=255)
                                                                                                                                              gAACATTCCGGGTACAGCGTTTGCAGGGCACTAATTAATGTCGCGCGCTGATATTCTGCCCCAGCACTCAGCAGTTGCAGCACCAGAAAATTACCGAAACGATCGATAGTAATACCCGGCAGGCCATCAGATTCTCCGGCGATTAAACg  <  4:372703/149‑1 (MQ=255)
                                                                                                                                              gAACATTCCGGGTACAGCGTTTGCAGGGCACTAATTAATGCCGCGCGCTGATATTCTGCCCCAGCACTCAGCAGTTGCAGCACCAGAAAATTACCGAAACGATCGATAGTAATACCCGGCAGGGCATCAGATTCTCCGGCGATTAAAc   >  1:336316/1‑148 (MQ=255)
                                                                                                                                              gAACATTCCGGGTACAGCGTTTGCAGGGCACTAATTAATGCCGCGCGCTGATATTCTGCCCCAGCACTCAGCAGTTGCAGCACCAGAAAATTACCGAAACGATCGATAGTAATACCCGGCAGGCCATCAGATTCTCCGGCGATTAAAc   <  2:192424/148‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
TCCGTGCTGAATATCCACCAGCAGTTTCATTCCGTGTTCTTCAATCGGCAGCAGGGCAGGTGGCAACTCGCCGGTGACGGGGCCCTGGGTCAGCTCCATTCCTTCTTTTTTACGTACCGCGACGTCGCTGCGATCGTAAATCGAACATTCCGGGTACAGCGTTTGCAGGGCACTAATTAATGCCGCGCGCTGATATTCTGCCCCAGCACTCAGCAGTTGCAGCACCAGAAAATTACCGAAACGATCGATAGTAATACCCGGCAGGCCATCAGATTCTCCGGCGATTAAACG  >  NC_000913/1029373‑1029663

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: ‑