Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 435,380 A→C K82Q (AAA→CAA)  nusB → transcription antitermination protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913435,3800AC100.0% 239.9 / NA 64K82Q (AAA→CAA) nusBtranscription antitermination protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (39/25);  total (39/25)

GTTCCTGGCTGAACAGGATGTAAAAGACGTTGACGTCCTGTACTTCCGTGAGCTGCTGGCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATA  >  NC_000913/435235‑435521
                                                                                                                                                 |                                                                                                                                             
gTTCCTGGCTGAACAGGATGTAAAAGACGTTGACGTCCTGTACTTCCGTGAGCTGCTGGCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAACAAGc                                                                                                                                           >  2:12026/1‑150 (MQ=255)
gTTCCTGGCTGAACAGGATGTAAAAGACGTTGACGTCCTGTACTTCCGTGAGCTGCTGGCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAACAAGCa                                                                                                                                          <  2:120257/151‑1 (MQ=255)
  tCCTGGCTGAACAGGATGTAAAAGACGTTGACGTCCTGTACTTCCGTGAGCTGCTGGCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAACAAGCAg                                                                                                                                         <  2:224172/150‑1 (MQ=255)
     tGGCTGAACAGGATGTAAAAGACGTTGACGTCCTGTACTTCCGTGAGCTGCTGGCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAACAAGCAGTACt                                                                                                                                     >  1:598088/1‑151 (MQ=255)
     tGGCTGAACAGGATGTAAAAGACGTTGACGTCCTGTACTTCCGTGAGCTGCTGGCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAACAAGCAGTACt                                                                                                                                     >  2:162872/1‑151 (MQ=255)
     tGGCTGAACAGGATGTAAAAGACGTTGACGTCCTGTACTTCCGTGAGCTGCTGGCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAACAAGCAGTACt                                                                                                                                     >  1:401620/1‑151 (MQ=255)
        cTGAACAGGATGTAAAAGACGTTGACGTCCTGTACTTCCGTGAGCTGCTGGCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAACAAGCAGTACTgcg                                                                                                                                  >  1:832957/1‑151 (MQ=255)
             cAGGATGTAAAAGACGTTGACGTCCTGTACTTCCGTGAGCTGCTGGCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAACAAGCAGTACTGCGCATTg                                                                                                                             >  2:37399/1‑151 (MQ=255)
              aGGATGTAAAAGACGTTGACGACCTGTACTTCCGTGAGCTGCTGGCCGGGGTGGCGATTAATACCGCACGCCTCCGCGGTGTAATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAACAAGCAGTACTGCGCATTgc                                                                                                                            <  2:365792/151‑1 (MQ=255)
               ggATGTAAAAGACGTTGACGTCCTGTCCTTCCGTGAGCTGCTGGCCGGGGTGGCGACTAATACCGCATCCCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAACAAGCAGTACTGCGCATTgcg                                                                                                                           <  2:107340/151‑1 (MQ=255)
               ggATGTAAAAGACGTTGACGTCCTGTACTTCCGTGAGCTGCTGGCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAACAAGCAGTACTGCGCATTgcg                                                                                                                           >  1:854828/1‑151 (MQ=255)
               ggATGTAAAAGACGTTGACGTCCTGTACTTCCGTGAGCTGCTGGCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAACAAGCAGTACTGCGCATTgcg                                                                                                                           >  2:723903/1‑151 (MQ=255)
                  tGTAAAAGACGTTGACGTCCTGTACTTCCGTGAGCTGCTGGCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAACAAGCAGTACTGCGCATTGCGCTg                                                                                                                        >  1:897107/1‑151 (MQ=255)
                     aaaaGACGTTGACGTCCTGTACTTCCGTGAGCTGCTGGCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAACAAGCAGTACTGCGCATTGCGCTGTa                                                                                                                      <  1:5053/150‑1 (MQ=255)
                      aaaGACGTTGACGTCCTGTACTTCCGTGAGCTGCTGGCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAACAAGCAGTACTGCGCATTGCGCTGTACg                                                                                                                    >  1:823147/1‑151 (MQ=255)
                       aaGACGTTGACGTCCTGTACTTCCGTGAGCTGCTGGCCGGGGTGGCGACTAATCCCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAACAAGCAGTACTGCGCATTGCGCTGTACGa                                                                                                                   <  1:434647/151‑1 (MQ=255)
                          aCGTTGACGTCCTGTACTTCCGTGAGCTGCTGGCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAACAAGCAGTACTGCGCATTGCGCTGTACGAACt                                                                                                                >  1:145543/1‑151 (MQ=255)
                             ttGACGTCCTGTACTTCCGTGAGCTGCTGGCCGGGGTGGCGACTAATACCGCATCCCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAACAAGCAGTACTGCGCATTGCGCTGTACGAACt                                                                                                                <  1:94660/148‑1 (MQ=255)
                             ttGACGTCCTGTACTTCCGTGAGCTGCTGGCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAACAAGCAGTACTGCGCATTGCGCTGTACGAACTGTc                                                                                                             >  2:449447/1‑151 (MQ=255)
                               gACGTCCTGTACTTCCGTGAGCTGCTGGCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAACAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCt                                                                                                            >  1:179657/1‑150 (MQ=255)
                                 cGTCCTGTACTTCCGTGAGCTGCTGGCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAACAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTaa                                                                                                          >  2:670858/1‑150 (MQ=255)
                                 cGTCCTGTACTTCCGTGAGCTGCTGGCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAACAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTaa                                                                                                          >  1:693184/1‑150 (MQ=255)
                                  gTCCTGTACTTCCGTGAGCTGCTGGCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAACAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTaaa                                                                                                         >  1:861697/1‑150 (MQ=255)
                                  gTCCTGTACTTCCGTGAGCTGCTGGCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAACAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAAc                                                                                                        >  1:771946/1‑151 (MQ=255)
                                      tGTACTTCCGTGAGCTGCTGGCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATCCCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAACAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACg                                                                                                       <  2:377931/148‑1 (MQ=255)
                                      tGTACTTCCGTGAGCTGCTGGCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAACAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGt                                                                                                      >  1:11348/1‑149 (MQ=255)
                                       gTACTTCCGTGAGCTGCTGGCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAACAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTa                                                                                                     <  2:359363/149‑1 (MQ=255)
                                          cTTCCGTGAGCTGCTGGCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTATAACAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGAt                                                                                                >  1:870215/1‑151 (MQ=255)
                                          cTTCCGTGAGCTGCTGGCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAACAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGAt                                                                                                >  2:28706/1‑151 (MQ=255)
                                           ttCCGTGAGCTGCTGGCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAACAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGAtg                                                                                               >  1:609239/1‑151 (MQ=255)
                                           ttCCGTGAGCTGCTGGCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAACAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGAtg                                                                                               >  2:643930/1‑151 (MQ=255)
                                               gTGAGCTGCTGGCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAACAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGcc                                                                                           >  2:169418/1‑151 (MQ=255)
                                               gTGAGCTGCTGGCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAACAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGcc                                                                                           >  1:247487/1‑151 (MQ=255)
                                                    ctgctgGCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAACAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATAc                                                                                       >  1:595442/1‑150 (MQ=255)
                                                           ccGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAACAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTgg                                                                                >  1:129414/1‑150 (MQ=255)
                                                             ggggTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAACAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGcc                                                                              >  2:518492/1‑150 (MQ=255)
                                                             ggggTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAACAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCa                                                                             <  2:355817/151‑1 (MQ=255)
                                                                  ggCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAACAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTaa                                                                         <  2:332528/150‑1 (MQ=255)
                                                                   gCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAACAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACg                                                                       <  2:612340/151‑1 (MQ=255)
                                                                     gACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAACAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGa                                                                      <  2:35319/150‑1 (MQ=255)
                                                                        tAATACCGCAAACCCCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAACAAGCAGTACTGCTCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGaa                                                                     <  1:610706/148‑1 (MQ=255)
                                                                         aaTACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAACAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCg                                                                  >  1:802430/1‑150 (MQ=255)
                                                                            aCCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAACAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCg                                                              >  2:819961/1‑151 (MQ=255)
                                                                            aCCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAACAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCg                                                              >  2:391305/1‑151 (MQ=255)
                                                                            aCCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAACAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCg                                                              >  1:306/1‑151 (MQ=255)
                                                                             ccGCATCCCTCGACGGACTGAGGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAACAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCg                                                              <  2:398881/150‑1 (MQ=255)
                                                                                cATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAACAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACt                                                          >  1:523315/1‑151 (MQ=255)
                                                                                 aTACCTCGACGGACTGATGAAGCCATCCCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAACAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTg                                                         <  1:321340/151‑1 (MQ=255)
                                                                                  tACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAACAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTgg                                                        >  2:239109/1‑151 (MQ=255)
                                                                                    ccTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAACAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCg                                                      <  2:302979/151‑1 (MQ=255)
                                                                                           ggACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAACAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGt                                               >  2:807649/1‑151 (MQ=255)
                                                                                                gatgaGGCCATGCCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAACAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGgc                                          <  2:911955/151‑1 (MQ=255)
                                                                                                      gCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAACAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGAGAAATCGTTCGGCGCa                                       >  2:47735/1‑148 (MQ=255)
                                                                                                                  ccGCCTGCTGGAAGAACTGGGACAGGTAGAACAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATa                          >  1:547560/1‑149 (MQ=255)
                                                                                                                          tgGAAGAACTGGGACAGGTAGAACAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCa                 >  2:26162/1‑150 (MQ=255)
                                                                                                                           ggAAGAACTGGGACAGGTAGAACAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGt                   <  1:219204/147‑1 (MQ=255)
                                                                                                                           ggAAGAACTGGGACAGGTAGAACAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAAc               <  1:902764/151‑1 (MQ=255)
                                                                                                                                aaCTGGGACAGGTAGAACAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGt          <  1:284056/151‑1 (MQ=255)
                                                                                                                                  cTGGGACAGGTAGAACAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTAc        >  2:432030/1‑151 (MQ=255)
                                                                                                                                      gACAGGTAGAACAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGa    <  2:148023/151‑1 (MQ=255)
                                                                                                                                       aCAGGTAGAACAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGAt   >  1:218431/1‑151 (MQ=255)
                                                                                                                                       aCAGGTAGAACAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGAt   <  1:171617/151‑1 (MQ=255)
                                                                                                                                       aCAGGTAGAACAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGAt   <  1:713653/151‑1 (MQ=255)
                                                                                                                                        cAGGTAGAACAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATa  <  2:282718/151‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                             
GTTCCTGGCTGAACAGGATGTAAAAGACGTTGACGTCCTGTACTTCCGTGAGCTGCTGGCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATA  >  NC_000913/435235‑435521

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