Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,406,402 T→G K80N (AAA→AAC lrhA ← transcriptional repressor of flagellar, motility and chemotaxis genes

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,406,4020TG96.8% 740.6 / 10.9 217K80N (AAA→AAClrhAtranscriptional repressor of flagellar, motility and chemotaxis genes
Reads supporting (aligned to +/- strand):  ref base T (2/5);  new base G (104/106);  total (106/111)
Fisher's exact test for biased strand distribution p-value = 4.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.78e-01

Only 100 of 217 total aligned reads displayed.

GTATCGGCAGATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAACAGTTCTTTCCCAACGAGTTGTTCCAG  >  NC_000913/2406307‑2406497
                                                                                               |                                                                                               
   tCGGCAGATTCAGCTGAAGCACCGATAGTTGATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATGTTCCTgg                                                                                          <  1:3281830/100‑1 (MQ=255)
    cGGCAGATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGGACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATGTTCCTggn                                                                                         <  1:4450215/100‑2 (MQ=255)
    cGGCAGATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATGTTCCTGGc                                                                                         <  1:2303511/100‑1 (MQ=255)
     ggCAGATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATGTTCCTGGCg                                                                                        <  1:3293683/100‑1 (MQ=255)
     ggCAGATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATGTTCCTGGCg                                                                                        >  1:2946767/1‑100 (MQ=255)
     ggCAGATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATGTTCCTGGCg                                                                                        >  1:1434286/1‑100 (MQ=255)
         gATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATGTTCCTGGCGTAAc                                                                                    >  2:800123/1‑100 (MQ=255)
          aTTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATGTTCCTGGCGTAAcc                                                                                   >  1:728449/1‑100 (MQ=255)
             cATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAGCATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATGTTCCTGGCGTAACCaag                                                                                <  1:1962133/100‑1 (MQ=255)
              aTCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATGTTCCTGGCGTAACCaaga                                                                               >  2:2299479/1‑100 (MQ=255)
              aTCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATGTTCCTGGCGTAACCaaga                                                                               >  1:3394815/1‑100 (MQ=255)
              aTCTGAAGCACCGATAGTTAATACGCCCTCAAGATTACTGAACATTAATGATGCGCAGGCCTCATCATTAAAACGCAGGATGTTCCTGGCGTAACCaaga                                                                               <  2:2037304/100‑1 (MQ=255)
                cTGAAGCACCGATAGTTAATACGCCCTGAAGATTATTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATGTTCCTGGCGTAACCaagaag                                                                             <  2:5425996/100‑1 (MQ=255)
                cTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATGTTCCTGGCGTAACCaagaag                                                                             >  2:5049704/1‑100 (MQ=255)
                cTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATGTTCCTGGCGTAACCaagaag                                                                             <  1:1937191/100‑1 (MQ=255)
                 tGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATGTTCCTGGCGTAACCAAGAAGt                                                                            >  2:1226679/1‑100 (MQ=255)
                  gAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATGTTCCTGGCGTAACCAAGAAGtt                                                                           >  2:3257879/1‑100 (MQ=255)
                   aaGCAGCGATAGTTAATACGCCCTGAAGATTACTGAAGAATAAGGATGAGAAGGCCTCATCATTAAAACGCAGGATGTTCCTGGCGTAACCAAGAAGTTg                                                                          <  2:3971619/100‑1 (MQ=255)
                   aaGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAACCGCAGGATGTTCCTGGCGTAACCAAGAAGTTg                                                                          <  2:4953770/100‑1 (MQ=255)
                   aaGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATGTTCCTGGCGTAACCAAGAAGTTg                                                                          <  2:4010506/100‑1 (MQ=255)
                    aGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATGTTCCTGGCGTAACCAAGAAGTTGa                                                                         >  1:563010/1‑100 (MQ=255)
                    aGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATGTTCCTGGCGTAACCAAGAAGTTGa                                                                         >  2:844544/1‑100 (MQ=255)
                    aGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATGTTCCTGGCGTAACCAAGAAGTTGa                                                                         <  1:4964859/100‑1 (MQ=255)
                      cACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATGTTCCTGGCGTAACCAAGAAGTTGAAt                                                                       <  2:4830603/100‑1 (MQ=255)
                       aCCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATGTTCCTGGCGTAACCAAGAAGTTGAATg                                                                      <  2:4272707/100‑1 (MQ=255)
                       aCCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATGTTCCTGGCGTAACCAAGAAGTTGAATg                                                                      >  1:1853997/1‑100 (MQ=255)
                        ccGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATGTTCCTGGCGTAACCAAGAAGTTGAAtag                                                                     >  2:4572601/1‑98 (MQ=255)
                        ccGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATGTTCCTGGCGTAACCAAGAAGTTGAATGc                                                                     >  2:3894021/1‑100 (MQ=255)
                        ccGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATGTTCCTGGCGTAACCAAGAAGTTGAATGc                                                                     <  1:2110267/100‑1 (MQ=255)
                        ccGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATGTTCCTGGCGTAACCAAGAAGTTGAATGc                                                                     <  1:4028503/100‑1 (MQ=255)
                        ccGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCACCATTAAAACGCAGGATGTTCCTGGCGTAACCAAGAAGTTGAATGc                                                                     >  1:623275/1‑100 (MQ=255)
                          gATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATGTTCCTGGCGTAACCAAGAAGTTGAATGCCa                                                                   >  2:4752193/1‑100 (MQ=255)
                          gATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATGTTCCTGGCGTAACCAAGAAGTTGAATGCCa                                                                   <  1:1033523/100‑1 (MQ=255)
                          gATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATGTTCCTGGCGTAACCAAGAAGTTGAATGCCa                                                                   <  1:3918289/100‑1 (MQ=255)
                            tAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATGTTCCTGGCGTAACCAAGAAGTTGAATGCCATg                                                                 <  2:4809881/100‑1 (MQ=255)
                             aGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATGTTCCTGGCGTAACCAAGAAGTTGAATGCCATGt                                                                >  2:2466767/1‑100 (MQ=255)
                             aGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATGTTCCTGGCGTAACCAAGAAGTTGAATGCCATGt                                                                >  2:1429159/1‑100 (MQ=255)
                              gTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATGTTCCTGGCGTAACCAAGAAGTTGAATGCCATGtt                                                               >  1:2992735/1‑100 (MQ=255)
                                tAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATGTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCa                                                             >  2:2601712/1‑100 (MQ=255)
                                 gatctgCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATGTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAg                                                            <  2:196647/95‑1 (MQ=255)
                                 aaTACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAg                                                            >  1:1693254/1‑100 (MQ=255)
                                 aaTACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATGTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAg                                                            <  1:1000461/100‑1 (MQ=255)
                                  aTACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATGTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGt                                                           <  1:4262952/100‑1 (MQ=255)
                                  aTACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATGTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGt                                                           <  2:4084869/100‑1 (MQ=255)
                                   tACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATGTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGtt                                                          <  2:3146846/100‑1 (MQ=255)
                                   tACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATGTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGtt                                                          >  1:2515898/1‑100 (MQ=255)
                                   tACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATGTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGtt                                                          >  1:652087/1‑100 (MQ=255)
                                    aCGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATGTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTa                                                         <  1:940340/100‑1 (MQ=255)
                                      gCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATGTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAgagatc                                                       >  1:198556/1‑95 (MQ=255)
                                      gCCCTGAAGAATACTGAACATTAATGATGAGCAGGACTCATCATTATAACGCAGGTTGTTCCGGGCGTAACCAAGAAGTTGATTGCCATGTTTAGTTAAc                                                       <  2:3609127/100‑1 (MQ=255)
                                       cccTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATGTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACa                                                      <  2:3482005/100‑1 (MQ=255)
                                        ccTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATGTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAg                                                     >  2:5390854/1‑100 (MQ=255)
                                        ccTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATGTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAg                                                     >  2:139149/1‑100 (MQ=255)
                                         cTGAAGATTACTGAACATTCATGATGAGCAGGCCTCATCATTCAAACGCAGGATGTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGt                                                    >  1:4745129/1‑100 (MQ=255)
                                            aaGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATGTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTttg                                                 <  1:810604/100‑1 (MQ=255)
                                              gATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATGTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTttgtt                                               <  1:2245495/100‑1 (MQ=255)
                                               aTTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTttgttg                                              <  1:3873103/100‑1 (MQ=255)
                                               aTTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATGTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTttgttg                                              <  1:4639038/100‑1 (MQ=255)
                                                ttACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATGTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGc                                             <  2:3363293/100‑1 (MQ=255)
                                                  aCTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATGTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAAAAGTTTGTTGCGa                                           >  2:1449010/1‑100 (MQ=255)
                                                     gAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATGTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCg                                        >  1:5420952/1‑100 (MQ=255)
                                                     gAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATGTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCg                                        >  1:5347045/1‑100 (MQ=255)
                                                       aCATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATGTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTg                                      <  1:4008487/100‑1 (MQ=255)
                                                        ttttAATGATGAGGCCTGCTCATCATTAAAACGCAGGATGTTCCTGGCGTAACCAAGAAGCTGAATGCCATGATCAGTTAACAGTTTGTTGCGACCGTGa                                     >  2:2990655/3‑100 (MQ=37)
                                                        cATTAATGATGAGCAGGCCTCATCATTCAAACGAAGGATGTTCCTGGCGTAACCTATGCGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGa                                     >  2:5484161/1‑100 (MQ=255)
                                                         aTTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATGTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGAc                                    >  2:1274563/1‑100 (MQ=255)
                                                            aaTGATGAGCAGGCCTCATCATTAAAACGCAGGATGTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAg                                 <  1:504446/100‑1 (MQ=255)
                                                             atgatgAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGc                                <  1:280947/100‑1 (MQ=255)
                                                             atgatgAGCAGGCCTCATCATTAAAACGCAGGATGTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGc                                <  1:2790054/100‑1 (MQ=255)
                                                              tgatgaGCAGGCCTCATCATTAAAACGCAGGATGTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCg                               <  1:3288427/100‑1 (MQ=255)
                                                               gatgaGCAGGCCTCAACATTAAAACGCAGGCTGTTCCTGGCGTAAGCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGa                              >  1:4842537/1‑100 (MQ=255)
                                                                atgaGCAGGCCTCATCATTAAAACGCAGGATGTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGAGACCGTGACGAGCGaa                             >  1:4577671/1‑100 (MQ=255)
                                                                 tgaGCAGGCGGCATCATTAAAACGCAGGATGTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAAc                            <  2:1677964/100‑1 (MQ=255)
                                                                 tgaGCAGGCCTCATCATTAAAACGCAGGATGTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAAc                            <  1:855767/100‑1 (MQ=255)
                                                                   aGCAGGCCTCATCATTAAAACGCAGGATGTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAACAg                          <  1:1386546/100‑1 (MQ=255)
                                                                     cAGGCCTCATCATTAAAACGCAGGATGTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAACAGtt                        >  1:3264804/1‑100 (MQ=255)
                                                                        gCCTCATCATTAAAACGCAGGATGTTCCTTGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGTGAACAGTTCtt                     >  2:1534719/1‑100 (MQ=255)
                                                                         ccTCATCATTAAAACGCAGGATGTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAACAGTTCttt                    <  2:1166598/100‑1 (MQ=255)
                                                                         ccTCATCATTAAAACGCAGGATGTTCCTGGCGTAACCAAGAAGTTGAATGCCATGATCAGTTAACAGGTTGTTGCGACCGTGACGAGCGAACAGTTCttt                    >  1:5106655/1‑100 (MQ=255)
                                                                          cccatcatTAAAACGCAGGATGTTCCTGGCGTAATCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAACAGTTCTTTc                   <  2:3274859/98‑1 (MQ=255)
                                                                             atcatTAAAACGCAGGATGTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAACAGTTCTTTCCCa                <  1:3181353/100‑1 (MQ=255)
                                                                             atcatTAAAACGCAGGATGTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAACAGTTCTTTCCCa                >  1:2223019/1‑100 (MQ=255)
                                                                              tcatTAAAACGCAGGATGTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAACAGTTCTTTCCCaa               >  1:2895688/1‑100 (MQ=255)
                                                                                atTAAAACGCAGGATGTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAACAGTTCTTTCCCAACg             >  1:3667576/1‑100 (MQ=255)
                                                                                atTAAAACGCAGGATGTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAACAGTTCTTTCCCAACg             <  2:3800230/100‑1 (MQ=255)
                                                                                atTAAAACGCAGGATGTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAACAGTTCTTTCCCAACg             <  2:3081328/100‑1 (MQ=255)
                                                                                atTAAAACGCAGGATGTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGGGACGAGCGAACGGTTCTTTCCCAACg             >  1:2527431/1‑100 (MQ=255)
                                                                                atTAAAACGCAGGATGTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCACCCGTGACGAGCGAACACTTCTTTCCCAACg             >  2:4306991/1‑100 (MQ=255)
                                                                                  tAAAACGCAGGATGTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAACAGTTCTTTCCCAACGAg           <  1:1170326/100‑1 (MQ=255)
                                                                                   aaaaCGCAGGATGTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAACAGTTCTTTCCCAACGAgt          <  2:3691647/100‑1 (MQ=255)
                                                                                   aaaaCGCAGGATGTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAACAGTTCTTTCCCAACGAgt          >  1:523380/1‑100 (MQ=255)
                                                                                   aaaaCGCAGGATGTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAACAGTTCTTTCCCAACGAgt          <  1:1291621/100‑1 (MQ=255)
                                                                                     aaCGCAGGATGTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAACAGTTCTTTCCCAACGAgttg        >  2:2909964/1‑100 (MQ=255)
                                                                                      aCGCAGGATGTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAACAGTTCTTTCCCAACGAgttgt       <  2:1604606/100‑1 (MQ=255)
                                                                                      aCGCAGGATGTTCCTGGCGTAACCAAGAAGCTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAACAGTTCTTTCCCAACGAgttgt       <  1:830481/100‑1 (MQ=255)
                                                                                         cAGGATGTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAACAGTTCTTTCCCAACGAGTTGTTcc    <  2:186096/100‑1 (MQ=255)
                                                                                          aGGATGTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAACAGTTCTTTCCCAACGAGTTGTTCCa   <  1:2461387/100‑1 (MQ=255)
                                                                                          aGGATGTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAACAGTTCTTTCCCAACGAGTTGTTCCa   >  2:4768765/1‑100 (MQ=255)
                                                                                           ggATGTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAACAGTTCTTTCCCAACGAGTTGTTCCAg  >  1:308749/1‑100 (MQ=255)
                                                                                           ggATGTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAACAGTTCTTTCCCAACGAGTTGTTCCAg  >  1:3805551/1‑100 (MQ=255)
                                                                                               |                                                                                               
GTATCGGCAGATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAACAGTTCTTTCCCAACGAGTTGTTCCAG  >  NC_000913/2406307‑2406497

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