Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,018,530 (C)6→7 coding (145/1128 nt) fliK → flagellar hook‑length control protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,018,5241.C99.2% 1084.6 / ‑0.5 267coding (139/1128 nt)fliKflagellar hook‑length control protein
Reads supporting (aligned to +/- strand):  ref base . (0/2);  new base C (143/122);  total (143/124)
Fisher's exact test for biased strand distribution p-value = 2.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

Only 100 of 296 total aligned reads displayed.

CACCACCACATTGCCTGGCGGCAAAGCCAGCGATGCTGCACAAGATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGG‑‑CCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTGTTTCCGACGCGCAACAAGCTAATTTACTG  >  NC_000913/2018424‑2018619
                                                                                                     |                                                                                                
caccaccacATTGCCTGGCGGCAAAGCCAGCGATGCTGCACAAGATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCgg                                                                                                   >  2:3379478/1‑101 (MQ=255)
   caccacATTGCCTGGCGGCAAAGCCAGCGATGCTGCACAAGATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGG‑‑ccc                                                                                              <  2:3378925/101‑1 (MQ=255)
     ccacATTGCCTGGCGGCAAAGCCAGCGATGCTGCACAAGATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGG‑‑ccccc                                                                                            <  2:3378783/101‑1 (MQ=255)
      cacaTTGCCTGGCGGCAAAGCCAGCGATGCTGCACAAGATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGG‑‑cccccc                                                                                           >  1:3379481/1‑101 (MQ=255)
        caTTGCCTGGCGGCAAAGCCAGCGATGCTGCACAAGATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGG‑‑ccccccca                                                                                         >  2:3379482/1‑99 (MQ=255)
        caTTGCCTGGCGGCAAAGCCAGCGATGCTGCACAAGATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGG‑‑ccccccca                                                                                         <  1:3379085/101‑3 (MQ=255)
        caTTGCCTGGCGGCAAAGCCAGCGATGCTGCACAAGATTTTCTCGCGTTGTTGAGCGAAGCATTAGAAGGCGAGACAACTACCGACAAAGCgg                                                                                                   <  2:3379147/93‑1 (MQ=255)
         aTTGCCTGGCGGCAAAGCCAGCGATGCTGCACAAGATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGG‑‑cccccccag                                                                                        <  2:3378939/101‑4 (MQ=255)
          ttGCCTGGCGGCAAAGCCAGCGATGCTGCACAAGATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGG‑‑cccccccagt                                                                                       <  2:3379054/101‑5 (MQ=255)
           tGCCTGGCGGCAAAGCCAGCGATGCTGCACAAGATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGG‑‑cccccccagtt                                                                                      >  2:3379486/1‑96 (MQ=255)
            gCCTGGCGGCAAAGCCAGCGATGCTGCACAAGATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGG‑‑ccc                                                                                              >  1:3379488/1‑92 (MQ=255)
              cTGGCGGCAAAGCCAGCGATGCTGCACAAGATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCCCCCCCAGTTGCt                                                                                    >  1:3379491/1‑101 (MQ=255)
               tGGCGGCAAAGCCAGCGATGCTGCACAAGATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCCCCCCCAGTTGCtg                                                                                   >  2:3379493/1‑101 (MQ=255)
                   ggcAAAGCCAGCGATGCTGCACAAGATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCCCCCCCAGTTGCtggtgg                                                                               >  2:3379500/1‑101 (MQ=255)
                   ggcAAAGCCAGCGATGCTGCACAAGATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCCCCCCCAGTTGCtggtgg                                                                               >  1:3379498/1‑101 (MQ=255)
                     cAAAGCCAGCGATGCTGCACAAGATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCCCCCCCAGTTGCTGGTGGCa                                                                             >  1:3379503/1‑101 (MQ=255)
                      aaaGCCAGCGATGCTGCACAAGATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCCCCCCCAGTTGCTGGTGGCaa                                                                            >  2:3379507/1‑101 (MQ=255)
                       aaGCCAGCGATGCTGCACAAGATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCCCCCCCAGTTGCTGGTGGCAAc                                                                           >  1:3379509/1‑101 (MQ=255)
                           cAGCGATGCTGCACAAGATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCCCCCCCAGTTGCTGGTGGCAACAGAt                                                                       >  1:3379514/1‑101 (MQ=255)
                           cAGCGATGCTGCACAAGATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCCCCCCCAGTTGCTGGTGGCAACAGAt                                                                       >  2:3379517/1‑101 (MQ=255)
                            aGCGATGCTGCACAAGATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCCCCCCCAGTTGCTGGTGGCAACAGATa                                                                      >  2:3379522/1‑101 (MQ=255)
                            aGCGATGCTGCACAAGATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCCCCCCCAGTTGCTGGTGGCAACAGATa                                                                      >  1:3379520/1‑101 (MQ=255)
                              cGATGCTGCACAAGATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCCCCCCCAGTTGCTGGTGGCAACAGATAAg                                                                    <  2:3378949/101‑1 (MQ=255)
                               gATGCTGCACAAGATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCCCCCCCAGTTGCTGGTGGCAACAGATAAGc                                                                   <  1:3378477/101‑1 (MQ=255)
                                aTGCTGCACAAGATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCCCCCCCAGTTGCTGGTGGCAACAGATAAGcc                                                                  >  1:3379525/1‑101 (MQ=255)
                                aTGCTGCACAAGATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCCCCCCCAGTTGCTGGTGGCAACAGATAAGcc                                                                  >  2:3379527/1‑101 (MQ=255)
                                 tgctgcACAAGATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCCCCCCCAGTTGCTGGTGGCAACAGATAAGccc                                                                 <  2:3379133/101‑1 (MQ=255)
                                  gctgcACAAGATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCa                                                                >  1:3379532/1‑101 (MQ=255)
                                  gctgcACAAGATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCa                                                                >  2:3379533/1‑101 (MQ=255)
                                   ctgcACAAGATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCAc                                                               <  2:3379023/101‑1 (MQ=255)
                                   ctgcACAAGATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCAc                                                               <  2:3379151/101‑1 (MQ=255)
                                   ctgcACAAGATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCAc                                                               <  1:3379036/101‑1 (MQ=255)
                                   ctgcACAAGATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCAc                                                               <  1:3379077/101‑1 (MQ=255)
                                     gcACAAGATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGa                                                             <  1:3379166/101‑1 (MQ=255)
                                        caAGATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACaa                                                          >  2:3379538/1‑101 (MQ=255)
                                        caAGATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACaa                                                          <  2:3379087/101‑1 (MQ=255)
                                          aGATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACaaaa                                                        <  2:3379177/101‑1 (MQ=255)
                                          aGATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACaaaa                                                        >  1:3379539/1‑101 (MQ=255)
                                           gATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAg                                                       <  1:3379155/101‑1 (MQ=255)
                                           gATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAg                                                       >  2:3379542/1‑101 (MQ=255)
                                           gATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAg                                                       <  2:3379012/101‑1 (MQ=255)
                                            aTTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAgg                                                      <  2:3379281/101‑1 (MQ=255)
                                             ttttCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGc                                                     <  1:3379064/101‑1 (MQ=255)
                                             ttttCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGc                                                     >  2:3379545/1‑101 (MQ=255)
                                             ttttCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGc                                                     <  1:3379067/101‑1 (MQ=255)
                                               ttCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCAc                                                               >  2:3379550/1‑89 (MQ=255)
                                               ttCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGa                                                   <  1:3379108/101‑1 (MQ=255)
                                               ttCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGa                                                   <  2:3379193/101‑1 (MQ=255)
                                                tctcGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAg                                                  >  1:3379551/1‑101 (MQ=255)
                                                 ctcGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAg                                                  >  2:3379554/1‑100 (MQ=255)
                                                  tcGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGcc                                                <  1:3378915/101‑1 (MQ=255)
                                                   cgcgTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCg                                               >  2:3379557/1‑101 (MQ=255)
                                                     cgTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCt                                             >  2:3379561/1‑101 (MQ=255)
                                                      gttgttGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGcc                                                <  2:3379231/97‑1 (MQ=255)
                                                      gttgttGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTg                                            >  1:3379568/1‑101 (MQ=255)
                                                      gttgttGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTg                                            >  1:3379565/1‑101 (MQ=255)
                                                        tgttgAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGAt                                          <  1:3379025/101‑1 (MQ=255)
                                                        tgttgAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGAt                                          >  2:3379570/1‑101 (MQ=255)
                                                        tgttgAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGG‑‑ccccccca                                                                                         >  1:3379572/1‑51 (MQ=255)
                                                           tgAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAg                                       <  1:3379171/101‑1 (MQ=255)
                                                           tgAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAg                                       <  2:3379097/101‑1 (MQ=255)
                                                           tgAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAg                                       >  1:3379575/1‑101 (MQ=255)
                                                            gAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGc                                      >  2:3379578/1‑101 (MQ=255)
                                                               cGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCCCCCCCAGTTGCTGGTGGCAACAGATaa                                                                     >  1:3379581/1‑67 (MQ=255)
                                                               cGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGat                                   <  1:3379245/101‑1 (MQ=255)
                                                                 aaGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGatat                                 >  2:3379586/1‑101 (MQ=255)
                                                                 aaGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGatat                                 >  1:3379587/1‑101 (MQ=255)
                                                                   gCATTAGCAGGCGAGACAACTACCGACAAAGCGGCCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTg                               <  1:3379072/101‑1 (MQ=255)
                                                                   gCATTAGCAGGCGAGACAACTACCGACAAAGCGGCCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTg                               >  1:3379592/1‑101 (MQ=255)
                                                                    cATTAGCAGGCGAGACAACTACCGACAAAGCGGCCCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGAt                                          <  1:3379056/90‑1 (MQ=255)
                                                                    cATTAGCAGGCGAGACAACTACCGACAAAGCGGCCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTGt                              >  2:3379593/1‑101 (MQ=255)
                                                                    cATTAGCAGGCGAGACAACTACCGACAAAGCGGCCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTGt                              <  1:3378834/101‑1 (MQ=255)
                                                                        aGCAGGCGAGACAACTACCGACAAAGCGGCCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTGTTTcc                          <  2:3379118/101‑1 (MQ=255)
                                                                        aGCAGGCGAGACAACTACCGACAAAGCGGCCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTGTTTcc                          >  1:3379597/1‑101 (MQ=255)
                                                                         gCAGGCGAGACAACTACCGACAAAGCGGCCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTg                                            <  2:3379103/82‑1 (MQ=255)
                                                                         gCAGGCGAGACAACTACCGACAAAGCGGCCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTGTTTCCg                         <  1:3378739/101‑1 (MQ=255)
                                                                         gCAGGCGAGACAACTACCGACAAAGCGGCCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTGTTTCCg                         >  2:3379600/1‑101 (MQ=255)
                                                                           aGGCGAGACAACTACCGACAAAGCGGCCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTGTTTCCGAc                       <  1:3378950/101‑1 (MQ=255)
                                                                            ggCGAGACAACTACCGACAAAGCGGCCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTGTTTCCGAcg                      >  1:3379603/1‑101 (MQ=255)
                                                                             gCGAGACAACTACCGACAAAGCGGCCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTGTTTCCGAcgc                     <  2:3378861/101‑1 (MQ=255)
                                                                                ataCAACTACCGACAAAGCGGCCCCCCCAGTTGCTGGTGGCAACAGAt                                                                       <  1:3379000/46‑1 (MQ=38)
                                                                                agaCAACTACCGACAAAGCGGCCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTGTTTCCGACGCGca                  >  1:3379606/1‑101 (MQ=255)
                                                                                  aCAACTACCGACAAAGCGGCCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTGTTTcc                          <  2:3378981/91‑1 (MQ=255)
                                                                                  aCAACTACCGACAAAGCGGCCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTGTTTCCGACGCGcaac                >  2:3379609/1‑101 (MQ=255)
                                                                                   cAACTACCGACAAAGCGGCCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTGTTTCCGACGCGcaaca               >  1:3379611/1‑101 (MQ=255)
                                                                                    aaCTACCGACAAAGCGGCCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGatat                                 >  2:3379616/1‑82 (MQ=255)
                                                                                    aaCTACCGACAAAGCGGCCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTGTTTCCGACGCGcaacaa              <  1:3379214/101‑1 (MQ=255)
                                                                                     aCTACCGACAAAGCGGCCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTGTTTCCGACGCGCAACAAg             >  1:3379617/1‑101 (MQ=255)
                                                                                      cTACCGACAAAGCGGCCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCg                                                    >  1:3379622/1‑61 (MQ=255)
                                                                                       tACCGACAAAGCGGCCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTGTTTCCGACGCGCAACAAGCt           >  2:3379624/1‑101 (MQ=255)
                                                                                       tACCGACAAAGCGGCCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTGTTTCCGACGCGCAACAAGCt           <  1:3378978/101‑1 (MQ=255)
                                                                                        aCCGACAAAGCGGCCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCAATATTGTTTCCGACGCGCAACAAGCTa          >  1:3379626/1‑101 (MQ=255)
                                                                                          cGACAAAGCGGCCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTGTTTCCGACGCGCAACAAGCTAAt        <  1:3379016/101‑1 (MQ=255)
                                                                                           gACAAAGCGGCCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGc                                                 >  1:3379630/1‑59 (MQ=255)
                                                                                            aCAAAGCGGCCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTGTTTCCGACGCGCAACAAGCt           <  1:3379234/96‑1 (MQ=255)
                                                                                             cAAAGCGGCCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTGTTTCCGACGCGCAACAAGCTAATTTa     <  1:3379334/101‑1 (MQ=255)
                                                                                              aaaGCGGCCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTGTTTCCGACGCGCAACAAGCTAATTTa     <  2:3378984/100‑1 (MQ=255)
                                                                                              aaaGCGGCCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTGTTTCCGACGCGCAACAAGCTAATTTAc    <  2:3379032/101‑1 (MQ=255)
                                                                                               aaGCGGCCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTGTTTCCGACGCGCAACAAGCTAATTTACt   >  1:3379633/1‑101 (MQ=255)
                                                                                               aaGCGGCCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTGTTTCCGACGCGCAACAAGCTAATTTACt   >  2:3379634/1‑101 (MQ=255)
                                                                                                     |                                                                                                
CACCACCACATTGCCTGGCGGCAAAGCCAGCGATGCTGCACAAGATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGG‑‑CCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTGTTTCCGACGCGCAACAAGCTAATTTACTG  >  NC_000913/2018424‑2018619

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