Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,138,489 A→G T297A (ACG→GCG)  flgK → hook‑filament junction protein 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,138,4890AG88.3% 727.1 / 74.6 283T297A (ACG→GCG) flgKhook‑filament junction protein 1
Reads supporting (aligned to +/- strand):  ref base A (17/16);  new base G (121/129);  total (138/145)
Fisher's exact test for biased strand distribution p-value = 8.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.89e-01

Only 100 of 286 total aligned reads displayed.

ATATTGAGATCCCGGAGAAATTACTGAATACCGGGTCGCTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATACGCTTGGACAACTGGCGCTGGCATTTGCCGAGGCTTTCAACACCCAACACAAAGCCGGATTTGATGCTAACGGCGATGCCGGTGAAGATTT  >  NC_000913/1138400‑1138580
                                                                                         |                                                                                           
atatTGAGATCCCGGAGAAATTACTGAATACCGGGTCGCTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATgcgc                                                                                          >  2:1887688/1‑93 (MQ=255)
 tatTGAGATCCCGGAGAAATTACTGAATACCGGGTCGCTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATGCGCt                                                                                         <  1:4017872/93‑1 (MQ=255)
 tatTGAGATCCCGGAGAAATTACTGAATACCGGGTCGCTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATACGCt                                                                                         >  1:4335029/1‑93 (MQ=255)
  atTGAGATCCCGGAGAAATTACTGAATACCGGGTCGCTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATGCGCtt                                                                                        <  2:388399/93‑1 (MQ=255)
     gagaTCCCGGAGAAATTACTGAATACCGGGTCGCTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATGCGCTTGGa                                                                                     <  2:1513715/93‑1 (MQ=255)
       gaTCCCGGAGAAATTACTGAATACCGGGTCGCTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATGCGCTTGGACa                                                                                   >  2:3631178/1‑93 (MQ=255)
        aTCCCGGAGAAATTACTGAATACCGGGTCGCTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATGCGCTTGGACaa                                                                                  >  2:1241357/1‑93 (MQ=255)
         tCCCGGAGAAATTACTGAATACCGGGTCGCTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATGCGCTTGGACAAc                                                                                 >  1:498625/1‑93 (MQ=255)
           ccGGAGAAATTACTGAATACCGGGTCGCTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATGCGCTTGGACAACTg                                                                               <  2:3123980/93‑1 (MQ=255)
           ccGGAGAAATTACTGAATACCGGGTCGCTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATGCGCTTGGACAACTg                                                                               <  1:4469988/93‑1 (MQ=255)
           ccGGAGAAATTACTGAATACCGGGTCGCTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATGCGCTTGGACAACTg                                                                               <  2:2863673/93‑1 (MQ=255)
           ccGGAGAAATTACTGAATACCGGGTCGCTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATACGCTTGGACAACTg                                                                               >  2:2609988/1‑93 (MQ=255)
            cGGAGAAATTACTGAATACCGGGTCGCTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATGCGCTTGGACAACTgg                                                                              <  2:6885194/93‑1 (MQ=255)
             ggAGAAATTACTGAATACCGGGTCGCTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATGCGCTTGGACAACTGgc                                                                             <  2:6258605/93‑1 (MQ=255)
             ggAGAAATTACTGAATACCGGGTCGCTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATGCGCTTGGACAACTGgc                                                                             <  2:4279397/93‑1 (MQ=255)
             ggAGAAATTACTGAATACCGGGTCGCTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATGCGCTTGGACAACTGgc                                                                             >  1:3574304/1‑93 (MQ=255)
             ggAGAAATTACTGAATACCGGGTCGCTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATGCGCTTGGACAACTGgc                                                                             >  1:745141/1‑93 (MQ=255)
               agaAATTACTGAATACCGGGTCGCTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATGCGCTTGGACAACTGgcgc                                                                           <  1:5895438/93‑1 (MQ=255)
               agaAATTACTGAATACCGGGTCGCTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATGCGCTTGGACAACTGgcgc                                                                           <  2:4517407/93‑1 (MQ=255)
               agaAATTACTGAATACCGGGTCGCTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATGCGCTTGGACAACTGgcgc                                                                           >  2:5362744/1‑93 (MQ=255)
                gaAATTACTGAATACCGGGTCGCTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATGCGCTTGGACAACTGGCGCt                                                                          >  1:3904699/1‑93 (MQ=255)
                 aaaTTACTGAATACCGGGTCGCTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATGCGCTTGGACAACTGGCGCTg                                                                         <  2:139093/93‑1 (MQ=255)
                  aaTTACTGAATACCGGGTCGCTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATGCGCTTGGACAACTGGCGCTgg                                                                        <  1:4555203/93‑1 (MQ=255)
                  aaTTACTGAATACCGGGTCGCTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATGCGCTTGGACAACTGGCGCTgg                                                                        >  2:5721793/1‑93 (MQ=255)
                   aTTACTGAATACCGGGTCGCTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATGCGCTTGGACAACTGGCGCTGGc                                                                       >  2:3705064/1‑93 (MQ=255)
                     tACTGAATACCGGGTCGCTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATGCGCTTGGACAACTGGCGCTGGCAt                                                                     <  2:6074848/93‑1 (MQ=255)
                      aCTGAATACCGGGTCGCTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATGCGCTTGGACAACTGGCGCTGGCAtt                                                                    >  2:6697244/1‑93 (MQ=255)
                      aCTGAATACCGGGTCGCTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATACGCTTGGACAACTGGCGCTGGCAtt                                                                    <  2:7256169/93‑1 (MQ=255)
                       cTGAATACCGGGTCGCTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATGCGCTTGGACAACTGGCGCTGGCAttt                                                                   <  1:5787368/93‑1 (MQ=255)
                       cTGAATACCGGGTCGCTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATGCGCTTGGACAACTGGCGCTGGCAttt                                                                   <  1:5647443/93‑1 (MQ=255)
                         gAATACCGGGTCGCTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATGCGCTTGGACAACTGGCGCTGGCATTTGc                                                                 <  1:3465126/93‑1 (MQ=255)
                           aTACCGGGTCGCTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATGCGCTTGGACAACTGGCGCTGGCATTTGCCg                                                               >  1:177404/1‑93 (MQ=255)
                           aTACCGGGTCGCTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATGCGCTTGGACAACTGGCGCTGGCATTTGCCg                                                               >  2:6552108/1‑93 (MQ=255)
                            tACCGGGTCGCTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATGCGCTTGGACAACTGGCGCTGGCATTTGCCGa                                                              <  2:2457844/93‑1 (MQ=255)
                            tACCGGGTCGCTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATGCGCTTGGACAACTGGCGCTGGCATTTGCCGa                                                              >  1:650850/1‑93 (MQ=255)
                             aCCGGGTCGCTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATGCGCTTGGACAACTGGCGCTGGCATTTGCCGAg                                                             >  2:3006097/1‑93 (MQ=255)
                                gggTCGCTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATGCGCTTGGACAACTGGCGCTGGCATTTGCCGAGGCt                                                          >  2:5322992/1‑93 (MQ=255)
                                 ggTCGCTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATGCGCTTGGACAACTGGCGCTGGCATTTGCCGAGGCtt                                                         <  2:1654703/93‑1 (MQ=255)
                                 ggTCGCTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATACGCTTGGACAACTGGCGCTGGCATTTGCCGAGGCtt                                                         >  1:2280838/1‑93 (MQ=255)
                                  gTCGCTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATGCGCTTGGACAACTGGCGCTGGCATTTGCCGAGGCttt                                                        <  1:931333/93‑1 (MQ=255)
                                   tCGCTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATGCGCTTGGACAACTGGCGCTGGCATTTGCCGAGGCTTTc                                                       >  1:6672761/1‑93 (MQ=255)
                                     gCTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATGCGCTTGGACAACTGGCGCTGGCATTTGCCGAGGCTTTCaa                                                     <  2:4032202/93‑1 (MQ=255)
                                     gCTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATGCGCTTGGACAACTGGCGCTAGCATTTGCCGAGGCTTTCaa                                                     <  2:5603140/93‑1 (MQ=255)
                                        gggCGGCATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATGCGCTTGGACAACTGGCGCTGGCATTTGCCGAGGCTTTCAacac                                                  >  2:5635298/1‑93 (MQ=255)
                                         ggcggcATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATACGCTTGGACAACTGGCGCTGGCATTTGCCGAGGCTTTCAACAcc                                                 >  1:2205527/1‑93 (MQ=255)
                                            ggcATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATGCGCTTGGACAACTGGCGCTGGCATTTGCCGAGGCTTTCAACACCCaa                                              >  2:2049864/1‑93 (MQ=255)
                                            ggcATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATGCGCTTGGACAACTGGCGCTGGCATTTGCCGAGGCTTTCAACACCCaa                                              >  2:7092382/1‑93 (MQ=255)
                                            ggcATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATGCGCTTGGACAACTGGCGCTGGCATTTGCCGAGGCTTTCAACACCCaa                                              >  2:5009410/1‑93 (MQ=255)
                                            ggcATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATGCGCTTGGACAACTGGCGCTGGCATTTGCCGAGGCTTTCAACACCCaa                                              >  2:836750/1‑93 (MQ=255)
                                             gcATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATGCGCTTGGACAACTGGCGCTGGCATTTGCCGAGGCTTTCAACACCCAac                                             >  2:5900710/1‑93 (MQ=255)
                                             gcATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATGCGCTTGGACAACTGGCGCTGGCATTTGCCGAGGCTTTCAACACCCAac                                             <  2:4208354/93‑1 (MQ=255)
                                             gcATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATGCGCTTGGACAACTGGCGCTGGCATTTGCCGAGGCTTTCAACACCCAac                                             <  2:1306414/93‑1 (MQ=255)
                                              cATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATGCGCTTGGACAACTGGCGCTGGCATTTGCCGAGGCTTTCAACACCCAaca                                            <  2:1956802/93‑1 (MQ=255)
                                               aTTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATGCGCTTGGACAACTGGCGCTGGCATTTGCCGAGGCTTTCAACACCCAacac                                           >  1:1024625/1‑93 (MQ=255)
                                               aTTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATGCGCTTGGACAACTGGCGCTGGCATTTGCCGAGGCTTTCAACACCCAacac                                           <  1:7372463/93‑1 (MQ=255)
                                               aTTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATGCGCTTGGACAACTGGCGCTGGCATTTGCCGAGGCTTTCAACACCCAacac                                           >  1:4900350/1‑93 (MQ=255)
                                               aTTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATGCGCTTGGACAACTGGCGCTGGCATTTGCCGAGGCTTTCAACACCCAacac                                           >  2:6641491/1‑93 (MQ=255)
                                                ttCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATGCGCTTGGACAACTGGCGCTGGCATTTGCCGAGGCTTTCAACACCCAAcgca                                          >  1:4088180/1‑93 (MQ=255)
                                                  cTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATGCGCTTGGACAACTGGCGCTGGCATTTGCCGAGGCTTTCAACACCCAACACaaa                                        <  1:5493650/93‑1 (MQ=255)
                                                  cTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATGCGCTTGGACAACTGGCGCTGGCATTTGCCGAGGCTTTCAACACCCAACACaaa                                        <  1:4597626/93‑1 (MQ=255)
                                                  cTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATGCGCTTGGACAACTGGCGCTGGCATTTGCCGAGGCTTTCAACACCCAACACaaa                                        <  1:1947215/93‑1 (MQ=255)
                                                   tGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATGCGCTTGGACAACTGGCGCTGGCATTTGCCGAGGCTTTCAACACCCAACACAAAg                                       >  1:5327398/1‑93 (MQ=255)
                                                    gACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATGCGCTTGGACAACTGGCGCTGGCATTTGCCGAGGCTTTCAACACCCAACACAAAGc                                      >  1:6237742/1‑93 (MQ=255)
                                                    gACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATGCGCTTGGACAACTGGCGCTGGCATTTGCCGAGGCTTTCAACACCCAACACAAAGc                                      <  2:3295776/93‑1 (MQ=255)
                                                     aCATTCCGTTCTCAGGATCTGGACCAGACGCGTAATGCGCTTGGACAACTGGCGCTGGCATTTGCCGAGGCTTTCAACACCCAACACAAAGcc                                     <  1:2723529/93‑1 (MQ=255)
                                                       aTTCCGTTCTCAGGATCTGGACCAGACGCGTAATGCGCTTGGACAACTGGCGCTGGCATTTGCCGAGGCTTTCAACACCCAACACAAAGCCgg                                   <  2:3259309/93‑1 (MQ=255)
                                                        ttCCGTTCTCAGGATCTGGACCAGACGCGTAATGCGCTTGGACAACTGGCGCTGGCATTTGCCGAGGCTTTCAACACCCAACACAAAGCCGGa                                  >  2:6004077/1‑93 (MQ=255)
                                                        ttCCGTTCTCAGGATCTGGACCAGACGCGTAATGCGCTTGGACAACTGGCGCTGGCATTTGCCGAGGCTTTCAACACCCAACACAAAGCCGGa                                  >  2:1751523/1‑93 (MQ=255)
                                                        ttCCGTTCTCAGGATCTGGACCAGACGCGTAATGCGCTTGGACAACTGGCGCTGGCATTTGCCGAGGCTTTCAACACCCAACACAAAGCCGGa                                  >  2:4397543/1‑93 (MQ=255)
                                                         tCCGTTCTCAGGATCTGGACCAGACGCGTAATGCGCTTGGACAACTGGCGCTGGCATTTGCCGAGGCTTTCAACACCCAACACAAAGCCGGAt                                 <  1:3156344/93‑1 (MQ=255)
                                                         tCCGTTCTCAGGATCTGGACCAGACGCGTAATACGCTTGGACAACTGGCGCTGGCATTTGCCGAGGCTTTCAACACCCAACACAAAGCCGGAt                                 <  1:4538862/93‑1 (MQ=255)
                                                           cGTTCTCAGGATCTGGACCAGACGCGTAATGCGCTTGGACAACTGGCGCTGGCATTTGCCGAGGCTTTCAACACCCAACACAAAGCCGGGttt                               >  1:6971105/1‑93 (MQ=255)
                                                           cGTTCTCAGGATCTGGACCAGACGCGTAATGCGCTTGGACAACTGGCGCTGGCATTTGCCGAGGCTTTCAACACCCAACACAAAGCCGGAttt                               <  2:1825184/93‑1 (MQ=255)
                                                           cGTTCTCAGGATCTGGACCAGACGCGTAATGCGCTTGGACAACTGGCGCTGGCATTTGCCGAGGCTTTCAACACCCAACACAAAGCCGGAttt                               >  1:3024612/1‑93 (MQ=255)
                                                             ttCTCAGGATCTGGACCAGACGCGTAATGCGCTTGGACAACTGGCGCTGGCATTTGCCGAGGCTTTCAACACCCAACACAAAGCCGGATTTGa                             >  1:4685063/1‑93 (MQ=255)
                                                             ttCTCAGGATCTGGACCAGACGCGTAATGCGCTTGGACAACTGGCGCTGGCATTTGCCGAGGCTTTCAACACCCAACACAAAGCCGGATTTGa                             >  2:512768/1‑93 (MQ=255)
                                                              tctcAGGATCTGGACCAGACGCGTAATGCGCTTGGACAACTGGCGCTGGCATTTGCCGAGGCTTTCAACACCCAACACAAAGCCGGATTTGAt                            <  1:3768168/93‑1 (MQ=255)
                                                              tctcAGGATCTGGACCAGACGCGTAATGCGCTTGGACAACTGGCGCTGGCATTTGCCGAGGCTTTCAACACCCAACACAAAGCCGGATTTGAt                            >  2:4734799/1‑93 (MQ=255)
                                                                tcAGGATCTGGACCAGACGCGTAATACGCTTGGACAACTGGCGCTGGCATTTGCCGAGGCTTTCAACACCCAACACAAAGCCGGATTTGATGc                          <  1:1314637/93‑1 (MQ=255)
                                                                  aGGATCTGGACCAGACGCGTAATGCGCTTGGACAACTGGCGCTGGCATTTGCCGAGGCTTTCAACACCCAACACAAAGCCGGATTTGATGCTa                        <  1:2883372/93‑1 (MQ=255)
                                                                   ggATCTGGACCAGACGCGTAATGCGCTTGGACAACTGGCGCTGGCATTTGCCGAGGCTTTCAACACCCAACACAAAGCCGGATTTGATGCTaa                       >  2:641746/1‑93 (MQ=255)
                                                                    gATCTGGACCAGACGCGTAATGCGCTTGGACAACTGGCGCTGGCATTTGCCGAGGCTTTCAACACCCAACACAAAGCCGGATTTGATGCTAAc                      >  1:2533965/1‑93 (MQ=255)
                                                                    gATCTGGACCAGACGCGTAATGCGCTTGGACAACTGGCGCTGGCATTTGCCGAGGCTTTCAACACCCAACACAAAGCCGGATTTGATGCTAAc                      <  2:5472271/93‑1 (MQ=255)
                                                                    gATCTGGACCAGACGCGTAATGCGCTTGGACAACTGGCGCTGGCATTTGCCGAGGCTTTCAACACCCAACACAAAGCCGGATTTGATGCTAAc                      <  1:6789911/93‑1 (MQ=255)
                                                                     aTCTGGACCAGACGCGTAATGCGCTTGGACAACTGGCGCTGGCATTTGCCGAGGCTTTCAACACCCAACACAAAGCCGGATTTGATGCTAACg                     <  2:110598/93‑1 (MQ=255)
                                                                      tCTGGACCAGACGCGTAATGCGCTTGGACAACTGGCGCTGGCATTTGCCGAGGCTTTCAACACCCAACACAAAGCCGGATTTGATGCTAACgg                    <  2:2265653/93‑1 (MQ=255)
                                                                        tGGACCAGACGCGTAATGCGCTTGGACAACTGGCGCTGGCATTTGCCGAGGCTTTCAACACCCAACACAAAGCCGGATTTGATGCTAACGGCg                  >  1:1484685/1‑93 (MQ=255)
                                                                         ggACCAGACGCGTAATGCGCTTGGACAACTGGCGCTGGCATTTGCCGAGGCTTTCAACACCCAACACAAAGCCGGATTTGATGCTAACGGCGa                 <  1:5383006/93‑1 (MQ=255)
                                                                         ggACCAGACGCGTAATACGCTTGGACAACTGGCGCTGGCATTTGCCGAGGCTTTCAACACCCAACACAAAGCCGGATTTGATGCTAACGGCGa                 <  1:703202/93‑1 (MQ=255)
                                                                            ccAGACGCGTAATGCGCTTGGACAACTGGCGCTGGCATTTGCCGAGGCTTTCAACACCCAACACAAAGCCGGATTTGATGCTAACGGCGATGc              <  2:4890443/93‑1 (MQ=255)
                                                                            ccAGACGCGTAATACGCTTGGACAACTGGCGCTGGCATTTGCCGAGGCTTTCAACACCCAACACAAAGCCGGATTTGATGCTAACGGCGATGc              <  2:7529059/93‑1 (MQ=255)
                                                                             cAGACGCGTAATGCGCTTGGACAACTGGCGCTGGCATTTGCCGAGGCTTTCAACACCCAACACAAAGCCGGATTTGATGCTAACGGCGATGcc             <  1:5131020/93‑1 (MQ=255)
                                                                              aGACGCGTAATGCGCTTGGACAACTGGCGCTGGCATTTGCCGAGGCTTTCAACACCCAACACAAAGCCGGATTTGATGCTAACGGCGATGCCg            <  2:2827804/93‑1 (MQ=255)
                                                                                  gcgTAATGCGCTTGGACAACTGGCGCTGGCATTTGCCGAGGCTTTCAACACCCAACACAAAGCCGGATTTGATGCTAACGGCGATGCCGGTGa        <  1:4198264/93‑1 (MQ=255)
                                                                                  gcgTAATACGCTTGGACAACTGGCGCTGGCATTTGCCGAGGCTTTCAACACCCAACACAAAGCCGGATTTGATGCTAACGGCGATGCCGGTGa        >  1:2060069/1‑93 (MQ=255)
                                                                                   cgTAATGCGCTTGGACAACTGGCGCTGGCATTTGCCGAGGCTTTCAACACCCAACACAAAGCCGGATTTGATGCTAACGGCGATGCCGGTGaa       <  2:5813826/93‑1 (MQ=255)
                                                                                   cgTAATGCGCTTGGACAACTGGCGCTGGCATTTGCCGAGGCTTTCAACACCCAACACAAAGCCGGATTTGATGCTAACGGCGATGCCGGTGaa       <  1:2429377/93‑1 (MQ=255)
                                                                                     tAATGCGCTTGGACAACTGGCGCTGGCATTTGCCGAGGCTTTCAACACCCAACACAAAGCCGGATTTGATGCTAACGGCGATGCCGGTGAAGa     >  1:2143459/1‑93 (MQ=255)
                                                                                      aaTGCGCTTGGACAACTGGCGCTGGCATTTGCCGAGGCTTTCAACACCCAACACAAAGCCGGATTTGATGCTAACGGCGATGCCGGTGAAGAt    <  1:5557698/93‑1 (MQ=255)
                                                                                       aTACGCTTGGACAACTGGCGCTGGCATTTGCCGAGGCTTTCAACACCCAACACAAAGCCGGATTTGATGCTAACGGCGATGCCGGTGAAGAtt   >  1:1906914/1‑93 (MQ=255)
                                                                                         |                                                                                           
ATATTGAGATCCCGGAGAAATTACTGAATACCGGGTCGCTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATACGCTTGGACAACTGGCGCTGGCATTTGCCGAGGCTTTCAACACCCAACACAAAGCCGGATTTGATGCTAACGGCGATGCCGGTGAAGATTT  >  NC_000913/1138400‑1138580

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