Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,879,763 G→A L4F (CTT→TTT)  ygbT ← multifunctional endonuclease Cas1, CRISPR adaptation protein; DNA repair enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,879,7630GA100.0% 226.1 / NA 63L4F (CTT→TTT) ygbTmultifunctional endonuclease Cas1, CRISPR adaptation protein; DNA repair enzyme
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (30/33);  total (30/33)

ATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTACAAGATCTATTAACGCTGGCGCGTCGTTGATGGTGAGCACACCTTCAAAGCA  >  NC_000913/2879615‑2879904
                                                                                                                                                    |                                                                                                                                             
atgatgCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAAcc                                                                                                                                             <  1:589118/151‑1 (MQ=255)
  gatgCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAACCa                                                                                                                                            <  1:806370/150‑1 (MQ=255)
      cAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAACCAGGTc                                                                                                                                        >  1:434912/1‑150 (MQ=255)
         gCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAACCAGGTCAtt                                                                                                                                     >  1:995341/1‑150 (MQ=255)
             ccGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAACCAGGTCATTttatt                                                                                                                                <  2:656548/151‑1 (MQ=255)
                 gCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAACCAGGTCATTTTATTacac                                                                                                                            >  2:209974/1‑151 (MQ=255)
                 gCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAACCAGGTCATTTTATTacac                                                                                                                            <  2:1038853/151‑1 (MQ=255)
                  ccAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAACCAGGTCATTTTATTACAcc                                                                                                                           >  1:1184575/1‑151 (MQ=255)
                      cAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAACCAGGTCATTTTATTACACCTCaa                                                                                                                       <  1:260310/151‑1 (MQ=255)
                       aGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAACCAGGTCATTTTATTACACCTCAAt                                                                                                                      <  2:1007877/151‑1 (MQ=255)
                       aGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAACCAGGTCATTTTATTACACCTCAAt                                                                                                                      <  2:544368/151‑1 (MQ=255)
                         gAATATGAGTGCGGATCGCTGTCTTGTCGATAAGTACAAACGCGCGATCTATTACATCGATGTGAGCAGATTGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAACCAGGTCATTTTATTACACCTCAATca                                                                                                                    <  1:462222/151‑1 (MQ=255)
                         gAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAACCAGGTCATTTTATTACACCTCAATc                                                                                                                     >  2:987934/1‑150 (MQ=255)
                          aaTATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAACCAGGTCATTTTATTACACCTCAATcac                                                                                                                   <  1:1148023/151‑1 (MQ=255)
                          aaTATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAACCAGGTCATTTTATTACACCTCAATc                                                                                                                     <  1:63552/149‑1 (MQ=255)
                           atatGTGTGCCGTTCTCAATCTTGTCGATAAGTACAAACTTGCCATCTATTTCATCGATCTGCCCATATTGCAGAAGGATCATGGCGTCGCGATCTTTGAGTGGAATGGGATTAAGGTGAAACCAGGTCATTTTCTTACCCCTCAATcc                                                                                                                    <  1:180996/149‑2 (MQ=255)
                           atatGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAACCAGGTCATTTTATTACACCTCAATcaca                                                                                                                  <  2:823886/151‑1 (MQ=255)
                            tatGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAACCAGGTCATTTTATTACACCTCAATcac                                                                                                                   >  2:913202/1‑149 (MQ=255)
                              tGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAACCAGGTCATTTTATTACACCTCAATCACAGTg                                                                                                               >  1:633342/1‑151 (MQ=255)
                                aGTGCAGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAACCAGGTCATTTTATTACACCTCAATCACAGTg                                                                                                               <  1:190278/149‑1 (MQ=255)
                                        tCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAACCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCaa                                                                                                        >  1:1151817/1‑148 (MQ=255)
                                            tGTCTTGTCGATAAGTACAAACGGGCCATCTATTACATCGATCTGCCCGTATTGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAACCAGGTCATTTTATTACACCTCAATCACAGCGGAGCCAAAGATAg                                                                                                  >  1:289451/1‑150 (MQ=255)
                                            tGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAACCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGc                                                                                                 >  1:638161/1‑151 (MQ=255)
                                             gTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAACCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGc                                                                                                 <  2:453002/150‑1 (MQ=255)
                                             gTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAACCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGc                                                                                                 <  1:422628/150‑1 (MQ=255)
                                                 tGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAACCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGc                                                                                             <  2:809118/150‑1 (MQ=255)
                                                  gTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAACCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGc                                                                                             >  2:204981/1‑149 (MQ=255)
                                                        aaGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAACCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATccc                                                                                     <  1:293673/151‑1 (MQ=255)
                                                         aGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAACCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATccc                                                                                     >  2:180113/1‑150 (MQ=255)
                                                         aGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAACCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATccc                                                                                     >  2:1082989/1‑150 (MQ=255)
                                                            aCAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAACCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATc                                                                                  <  2:430489/150‑1 (MQ=255)
                                                            aCAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAACCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCg                                                                                 <  1:241920/151‑1 (MQ=255)
                                                             cAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAACCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCg                                                                                 >  1:1069698/1‑150 (MQ=255)
                                                              aaaCGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAACCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGAt                                                                               >  2:941605/1‑151 (MQ=255)
                                                                  gcgcCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAACCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAg                                                                           <  2:269677/151‑1 (MQ=255)
                                                                   cgcCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAACCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGc                                                                          <  1:181367/151‑1 (MQ=255)
                                                                     ccATCTATTACCTCGATCTGCCCATATTGCAGAAAGATCATGGAGACGCGCTCTTTGAGTGGAATGGGATTAAGGGGACACCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTg                                                                        <  1:42140/151‑1 (MQ=255)
                                                                     ccATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAACCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTg                                                                        <  1:397581/151‑1 (MQ=255)
                                                                     ccATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAACCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTg                                                                        >  1:15502/1‑151 (MQ=255)
                                                                       aTCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAACCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGc                                                                      >  2:296556/1‑151 (MQ=255)
                                                                             tACATCGATCTGCCCATATTGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAACCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAAt                                                                 >  2:76979/1‑150 (MQ=255)
                                                                              aCATCGATCTGCCCATATTGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAACCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATAc                                                               >  1:1093536/1‑151 (MQ=255)
                                                                               cATCGATCTGCCCATATTGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAACCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATAcc                                                              >  2:651453/1‑151 (MQ=255)
                                                                                aTCGATCTGCCCATATTGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAACCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCt                                                             >  2:234628/1‑151 (MQ=255)
                                                                                     tCTGCCCATATTGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAACCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTg                                                        >  1:762616/1‑151 (MQ=255)
                                                                                        gCCCATATTGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAACCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTg                                                           <  2:297278/145‑1 (MQ=255)
                                                                                           cATATTGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAACCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTACa                                                    >  1:155292/1‑149 (MQ=255)
                                                                                                      aaGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAACCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTACAAGATCTATTAACg                                       >  1:965237/1‑151 (MQ=255)
                                                                                                        gATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAACCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTACAAGATCTATTAACGCt                                     <  1:296547/151‑1 (MQ=255)
                                                                                                         aTCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAACCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCAAATCCCTTGCTGTACAAGATCTATTAACGCtt                                    <  1:522450/151‑2 (MQ=255)
                                                                                                           cATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAACCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTACAAGATCTATTAACGCTGgc                                  <  1:155155/151‑1 (MQ=255)
                                                                                                           cATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAACCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTACAAGATCTATTAACGCTGgc                                  >  2:381069/1‑151 (MQ=255)
                                                                                                                agaCGCGATCTTTGAGTGGAATGGGATTAAGGGGAAACCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTACAAGATCTATTAACGCTGGCGcgtc                             >  2:1157238/1‑151 (MQ=255)
                                                                                                                 gaCGCGATCTTTAAGTGGAATGGGATTAAGGGGAAACCAGGTCATTTTATTCCACCTCAATCGCAGTGGAGCCAAAGATAGCAGGCCACATCCCATCGATTTAGCTGGGCCAATCCCTTGCTGTACAAGATCTATTAGCGCTGGcgcg                               <  1:465712/148‑1 (MQ=255)
                                                                                                                   cgcgATCTTTGAGTGGAATGGGATTAAGGGGAAACCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTACAAGATCTATTAACGCTGGCGCGTCGtt                          <  2:988999/151‑1 (MQ=255)
                                                                                                                     cgATCTTTGAGTGGAATGGGATTAAGGGGAAACCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTACAAGATCTATTAACGCTGGCGCGTCGTTg                         >  2:688017/1‑150 (MQ=255)
                                                                                                                      gATCTTTGAGTGGAATGGGATTAAGGGGAAACCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTACAAGATCTATTAACGCTGGCGCGTCGTTGAt                       >  2:615244/1‑151 (MQ=255)
                                                                                                                           ttGAGTGGAATGGGATTAAGGGGAAACCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTACACGATCTATTAACGCTGGCGCGTCGTTGATgg                     >  2:126589/1‑148 (MQ=255)
                                                                                                                           ttGAGTGGAATGGGATTAAGGGGAAACCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTACAAGATCTATTAACGCTGGCGCGTCGTTGATg                      <  2:259903/147‑1 (MQ=255)
                                                                                                                           ttGAGTGGAATGGGATTAAGGGGAAACCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTACAAGATCTATTAACGCTGGCGCGTCGTTGATg                      <  2:262250/147‑1 (MQ=255)
                                                                                                                              aGTGGAATGGGATTAAGGGGAAACCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTACAAGATCTATTAACGCTGGCGCGTCGTTGATGGTGAGca               <  2:764/151‑1 (MQ=255)
                                                                                                                              aGTGGAATGGGATTAAGGGGAAACCAGGTCATTTAATTACACCTCAATCACAGTAGAGCCAAAGATAGCAAGTCACATACCATCGATTTAGCTGGTCCAATACCTTGTTGTACAAGATCTATTTACGCTGGCGTGTCGTTGATGGTGAGc                >  1:601102/1‑150 (MQ=255)
                                                                                                                                            aaGGGGAAACCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTACAAGATCTATTAACGCTGGCGCGTCGTTGATGGTGAGCACACCTTCAAAGCa  <  2:632278/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                             
ATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTACAAGATCTATTAACGCTGGCGCGTCGTTGATGGTGAGCACACCTTCAAAGCA  >  NC_000913/2879615‑2879904

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