Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 668,691 C→A G9* (GGA→TGA)  rlmH ← 23S rRNA m(3)Psi1915 pseudouridine methyltransferase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913668,6910CA100.0% 250.7 / NA 65G9* (GGA→TGA) rlmH23S rRNA m(3)Psi1915 pseudouridine methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (29/36);  total (29/36)

CACCCTCTTTGTCGAGTATGCGCTTGATGTCCGCATTCTTGCCGCGTTTTCCGGCCGGAATTTCAATCAGCTCGAAGGGCATATCTTTCGGAAAACGACGCAGGTACTCGGTAAAACCGGTTTGTACCCAGTCCGGCATTTTCGTTCCCACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGCGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCCCAAATCCACGACAATCCAGTCGGCGCTGTTTTCACCTTCTAC  >  NC_000913/668544‑668832
                                                                                                                                                   |                                                                                                                                             
cACCCTCTTTGTCGAGTATGCGCTTGATGTCCGCATTCTTGCCGCGTTTTCCGGCCGGAATTTCAATCAGCTCGAAGGGCATATCTTTCGGAAAACGACGCAGGTACTCGGTAAAACCGGTTTGTACCCAGTCCGGCATTTTCGTTCacac                                                                                                                                            <  1:690185/151‑1 (MQ=255)
       tttGTCGAGTATGCGCTTGATGTCCGCATTCTTGCCGCGTTTTCCGGCCGGAATTTCAATCAGCTCGAAGGGCATATCTTTCGGAAAACGACGCAGGTACTCGGTAAAACCGGTTTGTACCCAGTCCGGCATTTTCGTTCACACGGCGAc                                                                                                                                      >  1:314522/1‑150 (MQ=255)
         tGTCGAGTATGCGCTTGATGTCCGCATTCTTGCCGCGTTTTCCGGCCGGAATTTCAATCAGCTCGAAGGGCATATCTTTCGGAAAACGACGCAGGTACTCGGTAAAACCGGTTTGTACCCAGTCCGGCATTTTCGTTCACACGGCGACa                                                                                                                                     >  1:728419/1‑149 (MQ=255)
         agTCGAGTATGCGCTTGATGTCCGCATTCTTGCCGCGTTTTCCGGCCGGAATTTCAATCAGCTCGAAGGGCATATCTTTCGGAAAACGACGCAGGTACTCGGTAAAACCGGTTTGTACCCAGTCCGGCATTTTCGTTCACACGGCGACa                                                                                                                                     >  1:561913/2‑149 (MQ=255)
            cGAGTATGCGCTTGATGTCCGCATTCTTGCCGCGTTTTCCGGCCGGAATTTCAATCAGCTCGAAGGGCATATCTTTCGGAAAACGACGCAGGTACTCGGTAAAACCGGTTTGTACCCAGTCCGGCATTTTCGTTCACACGGCGACAAGTTg                                                                                                                                >  2:887860/1‑151 (MQ=255)
              aGTATGCGCTTGATGTCCGCATTCTTGCCGCGTTTTCCGGCCGGAATTTCAATCAGCTCGAAGGGCATATCTTTCGGAAAACGACGCAGGTACTCGGTAAAACCGGTTTGTACCCAGTCCGGCATTTTCGTTCACACGGCGACAAGTTGc                                                                                                                               >  1:961112/1‑150 (MQ=255)
                   gcgcTTGATGTCCGCATTCTTGCCGCGTTTTCCGGCCGGAATTTCAATCAGCTCGAAGGGCATATCTTTCGGAAAACGACGCAGGTACTCGGTAAAACCGGTTTGTACCCAGTCCGGCATTTTCGTTCACACGGCGACAAGTTGCAGCtt                                                                                                                          <  1:1150146/150‑1 (MQ=255)
                   gcgcTTGATGTCCGCATTCTTGCCGCGTTTTCCGGCCGGAATTTCAATCAGCTCGAAGGGCATATCTTTCGGAAAACGACGCAGGTACTCGGTAAAACCGGTTTGTACCCAGTCCGGCATTTTCGTTCACACGGCGACAAGTTGCAGCTTc                                                                                                                         >  2:112134/1‑151 (MQ=255)
                       ttGATGTCCGCATTCTTGCCGCGTTTTCCGGCCGGAATTTCAATCAGCTCGAAGGGCATATCTTTCGGAAAACGACGCAGGTACTCGGTAAAACCGGTTTGTACCCAGTCCGGCATTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGc                                                                                                                     >  1:560090/1‑151 (MQ=255)
                        tGATGTCCGCATTCTTGCCGCGTTTTCCGGCCGGAATTTCAATCAGCTCGAAGGGCATATCTTTCGGAAAACGACGCAGGTACTCGGTAAAACCGGTTTGTACCCAGTCCGGCATTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGc                                                                                                                     >  1:612081/1‑150 (MQ=255)
                         gATGTCCGCATTCTTGCCGCGTTTTCCGGCCGGAATTTCAATCAGCTCGAAGGGCATATCTTTCGGAAAACGACGCAGGTACTCGGTAAAACCGGTTTGTACCCAGTCCGGCATTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGCAt                                                                                                                   <  2:973299/151‑1 (MQ=255)
                            gTCCGCATTCTTGCCGCGTTTTCCGGCCGGAATTTCAATCAGCTCGAAGGGCATATCTTTCGGAAAACGACGCAGGTACTCGGTAAAACCGGTTTGTACCCAGTCCGGCATTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGCATTa                                                                                                                 >  2:26183/1‑150 (MQ=255)
                             tCCGCATTCTTGCCGCGTTTTCCGGCCGGAATTTCAATCAGCTCGAAGGGCATATCTTTCGGAAAACGACGCAGGTACTCGGTAAAACCGGTTTGTACCCAGTCCGGCATTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGCATTAAc                                                                                                               >  1:1133712/1‑151 (MQ=255)
                             tCCGCATTCTTGCCGCGTTTTCCGGCCGGAATTTCAATCAGCTCGAAGGGCATATCTTTCGGAAAACGACGCAGGTACTCGGTAAAACCGGTTTGTACCCAGTCCGGCATTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGCATTAAc                                                                                                               >  1:594366/1‑151 (MQ=255)
                                 cATTCTTGCCGCGTTTTCCGGCCGGAATTTCAATCAGCTCGAAGGGCATATCTTTCGGAAAACGACGCAGGTACTCGGTAAAACCGGTTTGTACCCAGTCCGGCATTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGCATTAACTCCa                                                                                                           <  1:374716/151‑1 (MQ=255)
                                  aTTCTTGCCGCGTTTTCCGGCCGGAATTTCAATCAGCTCGAAGGGCATATCTTTCGGAAAACGACGCAGGTACTCGGTAAAACCGGTTTGTACCCAGTCCGGCATTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGCATTAACTc                                                                                                             >  1:1100778/1‑148 (MQ=255)
                                    tCTTGCCGCGTTTTCCGGCCGGAATTTCAATCAGCTCGAAGGGCATATCTTTCGGAAAACGACGCAGGTACTCGGTAAAACCGGTTTGTACCCAGTCCGGCATTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGCATTAACTCCaga                                                                                                         <  1:325191/150‑1 (MQ=255)
                                       ttccGCGTTTTCCGGCCGTAATTTCAATCAGCTCGAAGGGCATATCTTTCGGAAAACGACGCAGGTACTCGGTAAAACCGGTTTGCACCCAGTCCGGCATTTTCGTTCACACGGCGACAAGTTGCAGCTTCTCGCATTATCTCCAGAGtt                                                                                                      <  2:540844/148‑1 (MQ=255)
                                       tGCCGCGTTTTCCGGCCGGAATTTCAATCAGCTCGAAGGGCATATCTTTCGGAAAACGACGCAGGTACTCGGTAAAACCGGTTTGTACCCAGTCCGGCATTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGCATTAACTCCagag                                                                                                        <  1:418039/148‑1 (MQ=255)
                                            cgTTTTCCGGCCGGAATTTCAATCAGCTCGAAGGGCATATCTTTCGGAAAACGACGCAGGTACTCGGTAAAACCGGTTTGTACCCAGTCCGGCATTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCa                                                                                                <  2:1002597/151‑1 (MQ=255)
                                                   cggccggAATTTCAATCAGCTCGAAGGGCATATCTTTCGGAAAACGACGCAGGTACTCGGTAAAACCGGTTTGTACCCAGTCCGGCATTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCAt                                                                                          >  2:25297/1‑150 (MQ=255)
                                                   cggccggAATTTCAATCAGCTCGAAGGGCATATCTTTCGGAAAACGACGCAGGTACTCGGTAAAACCGGTTTGTACCCAGTCCGGCATTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATa                                                                                         <  2:497134/151‑1 (MQ=255)
                                                     gccggAATTTCAATCAGCTCGAAGGGCATATCTTTCGGAAAACGACGCAGGTACTCGGTAAAACCGGTTTGTACCCAGTCCGGCATTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATAc                                                                                        >  2:1135337/1‑150 (MQ=255)
                                                     gccggAATTTCAATCAGCTCGAAGGGCATATCTTTCGGAAAACGACGCAGGTACTCGGTAAAACCGGTTTGTACCCAGTCCGGCATTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATAc                                                                                        >  2:859171/1‑150 (MQ=255)
                                                      ccggAATTTCAATCAGCTCGAAGGGCATATCTTTCGGAAAACGACGCAGGTACTCGGTAAAACCGGTTTGTACCCAGTCCGGCATTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAg                                                                                      <  2:742700/151‑1 (MQ=255)
                                                       cggAATTTCAATCAGCTCGAAGGGCATATCTTTCGGAAAACGACGCAGGTACTCGGTAAAACCGGTTTGTACCCAGTCCGGCATTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAgg                                                                                     <  1:668526/151‑1 (MQ=255)
                                                       cggAATTTCAATCAGCTCGAAGGGCATATCTTTCGGAAAACGACGCAGGTACTCGGTAAAACCGGTTTGTACCCAGTCCGGCATTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAgg                                                                                     <  2:801779/151‑1 (MQ=255)
                                                                aaTCAGCTCGAAGGGCATATCTTTCGGAAAACGACGCAGGTACTCGGTAAAACCGGTTTGTACCCAGTCCGGCATTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGCGACGGctc                                                                            >  1:94367/1‑151 (MQ=255)
                                                                         gAAGGGCATATCTTTCGGAAAACGACGCAGGTACTCGGTAAAACCGGTTTGTACCCAGTCCGGCATTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGCGACGGCTCTCTTCCTGc                                                                   >  1:893332/1‑151 (MQ=255)
                                                                               cATATCTTTCGGAAAACGACGCAGGTACTCGGTAAAACCGGTTTGTACCCAGTCCGGCATTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGCGACGGCTCTCTTCCTGCATGACa                                                             <  1:112160/151‑1 (MQ=255)
                                                                                      ttCGGAAAACGACGCAGGTACTCGGTAAAACCGGTTTGTACCCAGTCCGGCATTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGCGACGGCTCTCTTCCTGCATGACATGGAct                                                       <  1:906832/150‑2 (MQ=255)
                                                                                      ttCGGAAAACGACGCAGGTACTCGGTAAAACCGGTTTGTACCCAGTCCGGCATTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGCGACGGCTCTCTTCCTGCATGACATGGACa                                                       >  2:746313/1‑150 (MQ=255)
                                                                                      ttCGGAAAACGACGCAGGTACTCGGTAAAACCGGTTTGTACCCAGTCCGGCATTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGCGACGGCTCTCTTCCTGCATGACATGGACa                                                       <  1:880024/150‑1 (MQ=255)
                                                                                      ttCGGAAAACGACGCAGGTACTCGGTAAAACCGGTTTGTACCCAGTCCGGCATTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGCGACGGCTCTCTTCCTGCATGACATGGACa                                                       >  2:916020/1‑150 (MQ=255)
                                                                                       tCGGAAAACGACGCAGGTACTCGGTAAAACCGGTTTGTACCCAGTCCGGCATTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGCGACGGCTCTCTTCCTGCATGACATGGAc                                                        <  2:1104925/148‑1 (MQ=255)
                                                                                       tCGGAAAACGACGCAGGTACTCGGTAAAACCGGTTTGTACCCAGTCCGGCATTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGCGACGGCTCTCTTCCTGCATGACATGGACaa                                                      <  1:375653/150‑1 (MQ=255)
                                                                                        cGGAAAACGACGCAGGTACTCGGTAAAACCGGTTTGTACCCAGTCCGGCATTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGCGACGGCTCTCTTCCTGCATGACATGGACAAt                                                     <  2:729586/150‑1 (MQ=255)
                                                                                        cGGAAAACGACGCAGGTACTCGGTAAAACCGGTTTGTACCCAGTCCGGCATTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGCGACGGCTCTCTTCCTGCATGACATGGACAAt                                                     <  1:505834/150‑1 (MQ=255)
                                                                                        cGGAAAACGACGCAGGTACTCGGTAAAACCGGTTTGTACCCAGTCCGGCATTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGCGACGGCTCTCTTCCTGCATGACATGGACAATc                                                    <  2:409897/151‑1 (MQ=255)
                                                                                         ggAAAACGACGCAGGTACTCGGTAAAACCGGTTTGTACCCAGTCCGGCATTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGCGACGGCTCTCTTCCTGCATGACATGGACAATca                                                   <  1:543716/151‑1 (MQ=255)
                                                                                                     aGGTACTCGGTAAAACCGGTTTGTACCCAGTCCGGCATTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGCGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCCCaaa                                        <  1:468685/150‑1 (MQ=255)
                                                                                                         aCTCGGTAAAACCGGTTTGTACCCAGTCCGGCATTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGCGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCCCAAATCCAc                                   <  1:859445/151‑1 (MQ=255)
                                                                                                          cTCGGTAAAACCGGTTTGTACCCAGTCCGGCATTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGCGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCCCAAATCCAc                                   <  2:484646/150‑1 (MQ=255)
                                                                                                            cGGTAAAACCGGTTTGTACCCAGTCCGGCATTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGCGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCCCAAATCCACGAc                                <  1:939295/151‑1 (MQ=255)
                                                                                                             ggTAAAACCGGTTTGTACCCAGTCCGGCATTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGCGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCCCAAATCCACGACa                               <  1:100002/151‑1 (MQ=255)
                                                                                                              gTAAAACCGGTTTGTACCCAGTCCGGCATTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGCGACGGCTCTCTTCCTGCATGACATGGCCAATCACATCGCCCAAATCCACGACaa                              <  2:624850/151‑1 (MQ=255)
                                                                                                              gTAAAACCGGTTTGTACCCAGTCCGGCATTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGCGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCCCAAATCCACGACaa                              <  2:1081506/151‑1 (MQ=255)
                                                                                                               tAAAACCGGTTTGTACCCAGTCCGGCATTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGCGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCCCAAATCCACGACAAt                             >  1:471049/1‑151 (MQ=255)
                                                                                                                aaaaCCGGTTTGTACCCAGTCCGGCATTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGCGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCCCAAATCCACGACAATc                            >  1:241448/1‑151 (MQ=255)
                                                                                                                  aaCCGGTTTGTACCCAGTCCGGCATTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGCGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCCCAAATCCACGaaatcgga                          <  2:1123506/151‑8 (MQ=255)
                                                                                                                  aaCCGGTTTGTACCCAGTCCGGCATTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGCGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCCCAAATCCACGACAATcc                           >  1:924223/1‑150 (MQ=255)
                                                                                                                   aCCGGTTTGTACCCAGTCCGGCATTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGCGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCCCAAATCCACGACAATcc                           >  2:522830/1‑149 (MQ=255)
                                                                                                                    ccGGTTTGTACCCAGTCCGGCATTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGCGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCCCAAATCCACGACAATCCAGt                        >  2:591765/1‑151 (MQ=255)
                                                                                                                            tACCCAGTCCGGCATTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGCGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCCCAAATCCACGACAATCCAGTCGGCGCTg                >  2:769198/1‑151 (MQ=255)
                                                                                                                             aCCCAGTCCGGCATTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGCGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCCCAAATCCACGACAATCCAGTCGGCGCTg                <  1:746541/150‑1 (MQ=255)
                                                                                                                              cccAGTCCGGCATTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGCGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCCCAAATCCACGACAATCCAGTCGGCGCTGtt              <  2:115342/151‑1 (MQ=255)
                                                                                                                               ccAGTCCGGCATTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGCGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCCCAAATCCACGACAATCCAGTCGGCGCTGttt             <  2:198158/151‑1 (MQ=255)
                                                                                                                                cAGTCCGGCATTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGCGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCCCAAATCCACGACAATCCAGTCGGCGCTGtttt            <  1:1147894/151‑1 (MQ=255)
                                                                                                                                 aGTCCGGCATTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGCGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCCCAAATCCACGACAATCCAGTCGGCGCTGTTTTc           >  2:367224/1‑151 (MQ=255)
                                                                                                                                    ccGGCATTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGCGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCCCAAATCCACGACAATCCAGTCGGCGCTGTTTTCAc         <  1:888149/150‑1 (MQ=255)
                                                                                                                                     cGGCATTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGCGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCCCAAATCCACGACAATCCAGTCGGCGCTGTTTTCACCt       <  2:728201/151‑1 (MQ=255)
                                                                                                                                     cGGCATTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGCGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCCCAAATCCACGACAATCCAGTCGGCGCTGTTTTCACCt       <  1:26190/151‑1 (MQ=255)
                                                                                                                                     cGGCATTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGCGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCCCAAATCCACGACAATCCAGTCGGCGCTGTTTTCACCt       >  2:942972/1‑151 (MQ=255)
                                                                                                                                       gCATTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGCGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCCCAAATCCACGACAATCCAGTCGGCGCTGTTTTCACCTTc     >  1:595892/1‑151 (MQ=255)
                                                                                                                                         aTTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGCGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCCCAAATCCACGACAATCCAGTCGGCGCTGTTTTCACCTTCt    >  1:144916/1‑150 (MQ=255)
                                                                                                                                           tttCGTTCACACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGCGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCCCAAATCCACGACAATCCAGTCGGCGCTGTTTTCACCTTCTAc  <  1:148754/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                             
CACCCTCTTTGTCGAGTATGCGCTTGATGTCCGCATTCTTGCCGCGTTTTCCGGCCGGAATTTCAATCAGCTCGAAGGGCATATCTTTCGGAAAACGACGCAGGTACTCGGTAAAACCGGTTTGTACCCAGTCCGGCATTTTCGTTCCCACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGCGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCCCAAATCCACGACAATCCAGTCGGCGCTGTTTTCACCTTCTAC  >  NC_000913/668544‑668832

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: