Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,423,367 C→T pseudogene (23/300 nt) rzpR → Rac prophage; putative endopeptidase RzpR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,423,3670CT100.0% 736.9 / NA 215pseudogene (23/300 nt)rzpRRac prophage; putative endopeptidase RzpR
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (100/115);  total (100/115)

Only 100 of 220 total aligned reads displayed.

CATCGCCTACAAAGAACAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTGATGATGT  >  NC_000913/1423269‑1423436
                                                                                                  |                                                                     
cATCGCCTACAAAGAACAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTgtgc                                                                     >  1:2067437/1‑101 (MQ=255)
 aTCGCCTACAAAGAACAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGTGCt                                                                    <  1:2067048/101‑1 (MQ=255)
 aTCGCCTACAAAGAACAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGTGCt                                                                    >  1:2067439/1‑101 (MQ=255)
  tCGCCTACAAAGAACAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGTGCTc                                                                   <  2:2066841/101‑1 (MQ=255)
   cGCCTACAAAGAACAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGTGCTCg                                                                  >  2:2067442/1‑101 (MQ=255)
    gCCTACAAAGAACAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGTGCTCGa                                                                 >  2:2067445/1‑101 (MQ=255)
     ccTACAAAGAACAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGTGCTCGAt                                                                >  1:2067446/1‑101 (MQ=255)
     ccTACAAAGAACAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGTGCTCGAt                                                                <  1:2067148/101‑1 (MQ=255)
     ccTACAAAGAACAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGTGCTCGAt                                                                <  2:2067023/101‑1 (MQ=255)
      cTACAAAGAACAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGTGCTCGATg                                                               >  1:2067448/1‑101 (MQ=255)
       tACAAAGAACAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGTGCTCGATGc                                                              <  1:2066856/101‑1 (MQ=255)
        aCAAAGAACAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGTGCTCGATGCa                                                             <  2:2067105/101‑1 (MQ=255)
          aaaGAACAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGTGCTCGATGCaaa                                                           <  2:2066920/101‑1 (MQ=255)
          aaaGAACAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGTGCTCGATGCaaa                                                           >  1:2067453/1‑101 (MQ=255)
          aaaGAACAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGTGCTCGATGCaaa                                                           >  1:2067451/1‑101 (MQ=255)
           aaGAACAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGTGCTCGATGCaaaa                                                          >  2:2067454/1‑101 (MQ=255)
            aGAACAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGTGCTCGATGCAAAAt                                                         >  2:2067456/1‑101 (MQ=255)
             gAACAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGTGCTCGATGCAAAATa                                                        >  1:2067459/1‑101 (MQ=255)
              aaCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGTGCTCGATGCAAAATac                                                       <  2:2066983/101‑1 (MQ=255)
              aaCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGTGCTCGATGCAAAATac                                                       <  2:2067149/101‑1 (MQ=255)
              aaCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGTGCTCGATGCAAAATac                                                       <  1:2066701/101‑1 (MQ=255)
              aaCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGTGCTCGATGCAAAATaa                                                       >  2:2067460/1‑100 (MQ=255)
               aCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGTGCTCGATGCAAAATaca                                                      >  1:2067463/1‑101 (MQ=255)
                   cgcgACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTgtgc                                                                     >  1:2067467/1‑82 (MQ=255)
                   cgcgACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGTGCTCGATGCAAAATACACGaa                                                  <  1:2066802/101‑1 (MQ=255)
                     cgACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGTGCTCGATGCAAAATACACGAAgg                                                >  1:2067469/1‑101 (MQ=255)
                      gACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGTGCTCGATGCAAAATACACGAAGGa                                               <  2:2067218/101‑1 (MQ=255)
                      gACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGTGCTCGATGCAAAATACACGAAGGa                                               <  2:2067094/101‑1 (MQ=255)
                      gACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGTGCTCGATGCAAAATACACGAAGGa                                               >  1:2067472/1‑101 (MQ=255)
                      gACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGTGCTCGATGCAAAATACACGAAGGa                                               <  1:2067185/101‑1 (MQ=255)
                      gACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGTGCTCGATGCAAAATACACGAAGGa                                               <  1:2067130/101‑1 (MQ=255)
                      gACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGTGCTAGATGCAAAATACACGAAGGa                                               <  2:2066950/101‑1 (MQ=255)
                       aCAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGTGCTCGATGCAAAATACACGAAGGAg                                              <  2:2067191/101‑1 (MQ=255)
                        cAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGTGCTCGATGCAAAATACACGAAGGAGt                                             >  2:2067476/1‑101 (MQ=255)
                        cAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGTGCTCGATGCAAAATACACGAAGGAGt                                             <  1:2067279/101‑1 (MQ=255)
                        cAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGTGCTCGATGCAAAATACACGAAGGAGt                                             >  1:2067475/1‑101 (MQ=255)
                         aaaaaTGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGTGCTCGATGCAAAATACACGAAGGAGtt                                            <  1:2067089/101‑1 (MQ=255)
                          aaaaTGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGTGCTCGATGCAAAATACACGAAGGAGTTa                                           >  2:2067479/1‑101 (MQ=255)
                          aaaaTGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGTGCTCGATGCAAAATACACGAAGGAGTTa                                           <  2:2067055/101‑1 (MQ=255)
                          aaaaTGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGTGCTCGATGCAAAATACACGAAGGAGTTa                                           <  2:2067269/101‑1 (MQ=255)
                          aaaaTGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGTGCTCGATGCAAAATACACGAAGGAGTTa                                           <  1:2067225/101‑1 (MQ=255)
                            aaTGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGTGCTCGATGCAAAATACACGAAGGAGTTAGc                                         <  1:2067114/101‑1 (MQ=255)
                             aTGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGTGCTCGATGCAAAATACACGAAGGAGTTAGCt                                        >  1:2067481/1‑101 (MQ=255)
                             aTGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGTGCTCGATGCAAAATACACGAAGGAGTTAGCt                                        <  2:2067070/101‑1 (MQ=25)
                              tGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGTGCTCGATGCAAAATACACGAAGGAGTTAGCTg                                       <  2:2066992/101‑1 (MQ=25)
                              tGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGTGCTCGATGCAAAATACACGAAGGAGTTAGCTg                                       <  2:2067030/101‑1 (MQ=25)
                              tGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGTGCTCGATGCAAAATACACGAAGGAGTTAGCTg                                       <  1:2066977/101‑1 (MQ=25)
                                ccAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGTGCTCGATGCAAAATACACGAAGGAGTTAGCTGAt                                     <  2:2067007/101‑1 (MQ=25)
                                 cAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGTGCTCGATGCAAAATACACGAAGGAGTTAGCTGAt                                     <  2:2067032/100‑1 (MQ=37)
                                 cAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGTGCTCGATGCAAAATACACGAAGGAGTTAGCTGATg                                    >  1:2067484/1‑101 (MQ=38)
                                 cAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGTGCTCGATGCAAAATACACGAAGGAGTTAGCTGATg                                    >  2:2067485/1‑101 (MQ=25)
                                  agagAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGTGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGc                                   >  1:2067487/1‑101 (MQ=25)
                                     gaACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGTGCTCGATGCAAAATACAGGAAGGAGTTAGCTGATGCGaa                                <  1:2066973/101‑1 (MQ=255)
                                      aaCTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGTGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGaaa                               >  1:2067490/1‑101 (MQ=38)
                                      aaCTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGTGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGaaa                               <  2:2067084/101‑1 (MQ=25)
                                       aCTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGTGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAg                              >  2:2067491/1‑101 (MQ=25)
                                       aCTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGTGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAg                              <  1:2067182/101‑1 (MQ=25)
                                        cTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGTGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGc                             <  2:2066995/101‑1 (MQ=25)
                                         tGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGTGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCt                            <  1:2067062/101‑1 (MQ=25)
                                          gAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGTGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTg                           >  1:2067494/1‑101 (MQ=38)
                                           aaGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGTGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGa                          >  2:2067495/1‑101 (MQ=38)
                                            aGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGTGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGaa                         <  1:2067086/101‑1 (MQ=25)
                                            aGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGTGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGaa                         <  1:2067020/101‑1 (MQ=25)
                                             gCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGTGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGaaa                        <  1:2066905/101‑1 (MQ=25)
                                              cTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGTGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGaaaa                       <  1:2067402/101‑1 (MQ=25)
                                               tGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGTGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAat                      <  1:2067195/101‑1 (MQ=25)
                                               tGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGTGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAat                      <  2:2066871/101‑1 (MQ=25)
                                               tGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGTGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAat                      >  2:2067501/1‑101 (MQ=25)
                                               tGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGTGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAat                      <  1:2066949/101‑1 (MQ=25)
                                               tGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGTGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAat                      >  1:2067500/1‑101 (MQ=25)
                                                ggCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGTGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAatg                     >  2:2067504/1‑101 (MQ=25)
                                                ggCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGTGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAatg                     <  1:2067011/101‑1 (MQ=18)
                                                 gCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGTGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAatga                    <  1:2067050/101‑1 (MQ=17)
                                                 gCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGTGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAatga                    >  2:2067507/1‑101 (MQ=17)
                                                 gCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGTGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAatga                    >  1:2067506/1‑101 (MQ=18)
                                                   gAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGTGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCg                                  <  1:2067324/85‑1 (MQ=18)
                                                     aCGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGTGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTg                           <  2:2067112/90‑1 (MQ=18)
                                                     aCGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGTGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGct                <  1:2067204/101‑1 (MQ=17)
                                                     aCGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGTGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGct                <  2:2067133/101‑1 (MQ=17)
                                                     aCGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGTGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGct                >  1:2067510/1‑101 (MQ=18)
                                                     aCGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGTGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGct                >  2:2067511/1‑101 (MQ=17)
                                                       gcgGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGTGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGctct              >  1:2067516/1‑101 (MQ=18)
                                                       gcgGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGTGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGctct              >  1:2067514/1‑101 (MQ=18)
                                                        cgGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGTGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTg             >  2:2067517/1‑101 (MQ=18)
                                                        cgGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGTGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTg             <  1:2066775/101‑1 (MQ=17)
                                                        cgGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGTGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTg             <  2:2066897/101‑1 (MQ=17)
                                                        cgGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGTGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTg             >  1:2067521/1‑101 (MQ=17)
                                                         ggCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGTGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGc            <  1:2067143/101‑1 (MQ=17)
                                                           cAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGTGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGt          <  2:2067051/101‑1 (MQ=17)
                                                           cAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGTGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGt          >  2:2067524/1‑101 (MQ=17)
                                                            aaTTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGTGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGtg         >  1:2067526/1‑101 (MQ=17)
                                                            aaTTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGTGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGtg         >  2:2067528/1‑101 (MQ=17)
                                                             aTTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGTGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGtga        >  1:2067529/1‑101 (MQ=17)
                                                               tACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGTGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGtgatg      <  2:2067166/101‑1 (MQ=17)
                                                               tACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGTGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGtgatg      >  1:2067531/1‑101 (MQ=17)
                                                                aCTGAGATGCAGATGCGTCAGCGTGATGTTGCTGTGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGc            >  2:2067534/1‑94 (MQ=18)
                                                                 cTGAGATGCAGATGCGTCAGCGTGATGTTGCTGTGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTgatgat    >  2:2067536/1‑101 (MQ=17)
                                                                  tGAGATGCAGATGCGTCAGCGTGATGTTGCTGTGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTGatgatg   >  2:2067539/1‑101 (MQ=17)
                                                                   gagaTGCAGATGCGTCAGCGTGATGTTGCTGTGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTGATGATGt  <  2:2067009/101‑1 (MQ=17)
                                                                   gagaTGCAGATGCGTCAGCGTGATGTTGCTGTGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTGATGATGt  <  1:2067108/101‑1 (MQ=17)
                                                                                                  |                                                                     
CATCGCCTACAAAGAACAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTGATGATGT  >  NC_000913/1423269‑1423436

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