Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 344,915 T→G I247R (ATA→AGA)  yahL → uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913344,9150TG93.0% 378.9 / 18.1 130I247R (ATA→AGA) yahLuncharacterized protein
Reads supporting (aligned to +/- strand):  ref base T (5/4);  new base G (61/58);  total (67/63)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.26e-01

Only 100 of 133 total aligned reads displayed.

TATCTACCCTGTTTTTCCTGAATCACTGGTGATGTTATTGCTCTCTGTTCTGGACGTAAAAACATACAGAGATGATGCGCGATTAAATGAATGGATATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAGTTCAGAAAGGGTTCCACCTGCAA  >  NC_000913/344820‑345014
                                                                                               |                                                                                                   
     aCCCTGTTTTTCCTGAATCACTGGTGATGTTATTGCTCTCTGTTCTGGACGTAAAAACATACAGAGATGATGCGCGATTAAATGAATGGAGATCTAATag                                                                                            <  1:419688/100‑1 (MQ=255)
       ccTGTTTTTCCTGAATCACTGGTGATGTTATTGCTCTCTGTTCTGGACGTAAAAACATACAGAGATGATGCGCGATTAAATGAATGGAGATCTAATagag                                                                                          >  1:1640737/1‑100 (MQ=255)
        cTGTTTTTCGTGAATCACTGGTGATGTTGTTGCGCTCTGTTGTGGGCGTAAAAAAATATAGACATGTTGCTCTATTAGATGGATGGAGATTTAATAGAGt                                                                                         >  2:587893/1‑100 (MQ=255)
          gTTTTTCCTGAATCAGTGGTGATGTTATTGCTCTCTGTTCAGGACCTAAAAACATACAGAGATGATGCGCGAATAAATGAATGGAGATCTAATAGAGTGc                                                                                       >  1:3543231/1‑100 (MQ=255)
           tttttCCTGAATCACTGGTGATGTTATTGCTCTCTGTTCTGGACGTAAAAACATACAGAGATGATGCGCGATTAAATGAATGGAGATCTAATAGAGTGCa                                                                                      <  1:3226415/100‑1 (MQ=255)
             tttCCTGAATCACTGGTGATGTTATTGCTCTCTGTTCTGGACGTAAAAACATACAGAGATGATGCGCGATTAAATGAATGGAGATCTAATAGAGTGCAgg                                                                                    >  1:1890646/1‑100 (MQ=255)
              ttCCTGAATCACTGGTGATGTTATTGCTCTCTGTTCTGGACGTAAAAACATACAGAGATGATGCGCGATTAAATGAATGGAGATCTAATAGAGTGCAGGa                                                                                   >  2:3396/1‑100 (MQ=255)
                ccTGAATCACTGGTGATGTTATTGCTCTCTGTTCTGGACGTAAAAACATACAGAGATGATGCGCGATTAAATGAATGGATATCTAATAGAGTGCAGGAAt                                                                                 >  1:384234/1‑100 (MQ=255)
                ccTGAATCACTGGTGATGTTATTGCTCTCTGTTCTGGACGTAAAAACATACAGAGATGATGCGCGATTAAATGAATGGAGATCTAATAGAGTGCAGGAAt                                                                                 <  1:3293787/100‑1 (MQ=255)
                 cTGAATCACTGGTGATGTTATTGCTCTCTGTTCTGGACGTAAAAACATACAGAGATGATGCGCGATTAAATGAATGGAGATCTAATAGAGTGCAGGAAtt                                                                                >  1:718365/1‑100 (MQ=255)
                    aagcACTGGTGATGTTATTGCTCTCTGTTCTGGACGTAAAAACATACAGAGATGATGCGCGATTAAATGAATGGAGATCTAGTAGAGTGCAGGCATTAgg                                                                             <  1:3404249/97‑1 (MQ=255)
                    aaTCACTGGTGATGTTATTGCTCTCTGTTCTGGACGTAAAAACATACAGAGATGATGCGCGATTAAATGAATGGCGATCTAATAGAGTGCAGGAATTAgg                                                                             <  2:1727717/100‑1 (MQ=255)
                      tCACTGGTGATGTTATTGCTCTCTGTTCTGGACGTAAAAACATACAGAGATGATGCGCGATTAAATGAATGGAGATCTAATAGAGTGCAGGAATTAGGTg                                                                           <  1:2613136/100‑1 (MQ=255)
                        aCTGGTGATGTTATTGCTCTCTGTTCTGGACGTAAAAACATACAGAGATGATGCGCGATTAAATGAATGGAGATCTAATAGAGTGCAGGAATTAGGTGAc                                                                         <  2:3470287/100‑1 (MQ=255)
                        aCTAGTGATGTTATTGCTCTCTGTTCTGGACGTAAAAACATACAGAGATGATGCGCGATTAAATGAATGGAGATCTAATAGAGTGCAGGAATTAGGTGAc                                                                         <  2:1254546/100‑1 (MQ=255)
                         cTGGTGATGTTATTGCTCTCTGTTCTGGACGTAAAAACATACAGAGATGATGCGCGATTAAATGAATGGAGATCTAATAGAGTGCAGGAATTAGGTGACa                                                                        <  1:2098228/100‑1 (MQ=255)
                           ggTGATGTTATTGCTCTCTGTTCTGGACGTAAAAACATACAGAGATGATGCGCGATTAAATGAATGGATATCTAATAGAGTGCAGGAATTAGGTGACAGa                                                                      >  1:3721581/1‑100 (MQ=255)
                           ggTGATGTTATTGCTCTCTGTTCTGGACGTAAAAACATACAGAGATGATGCGCGATTAAATGAATGGAGATCTAATAGAGTGCAGGAATTAGGTGACAGa                                                                      >  1:522974/1‑100 (MQ=255)
                           ggTGATGTTATTGCTCTCTGTTCTGGACGTAAAAACATACAGAGATGATGCGCGATTAAATGAATGGAGATCTAATAGAGTGCAGGAATTAGGTGACAGa                                                                      >  1:3521419/1‑100 (MQ=255)
                           ggTGATGTTATTGCTCTCTGTTCTGGACGTAAAAACATACAGAGATGATGAGCGATTAAATGAATGGAGATCTAATAGAGTGCAGGAATTAGGTGACAGa                                                                      >  2:3602485/1‑100 (MQ=255)
                             tGATGTTATTGCTCTCTGTTCTGGACGTAAAAACATACAGAGATGATGCGCGATTAAATGAATGGAGATCTAATAGAGTGCAGGAATTAGGTGACAGata                                                                    >  2:2354135/1‑100 (MQ=255)
                              gATGTTATTGCTCTCTGTTCTGGACGTAAAAACATACAGAGATGATGCGCGATTAAATGAATGGAGATCTAATAGAGTGCAGGAATTAGGTGACAGatat                                                                   >  2:2100934/1‑100 (MQ=255)
                               aTGTTATTGCTCTCTGTTCTGGACGTAAAAACATACAGAGATGATGCGCGATTAAATGAATGGAGATCTAATAGAGTGCAGGAATTAGGTGACAGAtatt                                                                  >  2:3647185/1‑100 (MQ=255)
                               aTGTTATTGCTCTCTGTTCTGGACGTAAAAACATACAGAGATGATGCGCGATTAAATGAATGGAGATCTAATAGAGTGCAGGAATTAGGTGACAGAtatt                                                                  <  2:3047720/100‑1 (MQ=255)
                                 gTTATTGCTCTCTGTTCTGGACGTAAAAACATACAGAGATGATGCGCGATTAAATGAATGGAGATCTAATAGAGTGCAGGAATTAGGTGACAGAtattat                                                                <  2:2975787/100‑1 (MQ=255)
                                  ttATTGCTCTCTGTTCTGGACGTAAAAACATACAGAGATGATGCGCGATTAAATGAATGGAGATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATc                                                               >  1:1291767/1‑100 (MQ=255)
                                   tATTGCTCTCTGTTCTGGACGTAAAAACATACAGAGATGATGCGCGATTAAATGAATGGAGATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATcc                                                              >  1:2410326/1‑100 (MQ=255)
                                   tATTGCTCTCTGGTCTTGACGTAAAAACATACGGAGATGATGCGCGATTAAATGAATGGAGATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATcc                                                              <  2:2401147/100‑1 (MQ=255)
                                       gCTCTCTGTTCTGGACGTAAAAACATACAGAGATGATGCGCGATTAAATGAATGGAGATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTa                                                          >  2:102138/1‑100 (MQ=255)
                                       gCTCTCTGTTCTGGACGTAAAAACATACAGAGATGATGCGCGATTAAATGAATGGAGATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTa                                                          <  2:3159842/100‑1 (MQ=255)
                                       gCTCTCTGTTCTGGACGTAAAAACATACAGAGATGATGCGCGATTAAATGAATGGAGATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTa                                                          <  1:3695745/100‑1 (MQ=255)
                                       gCTCTCTGTTCTGGACGTAAAAACATACAGAGATGATGCGCGATTAAATGAATGGAGATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTa                                                          <  2:3218685/100‑1 (MQ=255)
                                         tctctGTTCTGGACGTAAAAACATACAGAGATGATGCGCGATTAAATGAATGGAGATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTaaa                                                        >  1:2809669/1‑100 (MQ=255)
                                         tctctGTTCTGGACGTAAAAACATACAGAGATGATGCGCGATTAAATGAATGGAGATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTaaa                                                        >  1:26410/1‑100 (MQ=255)
                                         tctatGTCCTGGACGTAAAACCAAACAGAGATGATGCGCGATTAAATGAATGGAAATCTAATAGAGTGCAGGAATTATTTGACAGATATTATGCTGTaaa                                                        >  2:2083190/1‑100 (MQ=255)
                                          ctctGTTCTGGACGTAAAAACATACAGAGATGATGCGCGATTAAATGAATGGAGATCTAATAGAGTGCAGGAATTAGGTGACAGATATGATCCTGTAAAt                                                       >  2:1220240/1‑100 (MQ=255)
                                            ctGTTCTGGACGTAAAAACATACAGAGATGATGCGCGATTAAATGAATGGAGATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATaa                                                     >  2:3101254/1‑100 (MQ=255)
                                              gTTCTGGACGTAAAAACATACAGAGATGATGCGCGATTAAATGAATGGAGATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGc                                                   <  2:1235188/100‑1 (MQ=255)
                                              gTTCTGGAAGTAAAAACATACAGAGATGATGCGCGATTAAATGAATGGAGATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGc                                                   <  2:3287252/100‑1 (MQ=255)
                                                 cTGGACGTAAAAACATACAGAGATGATGCGCGATTAAATGAATGGAGATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACg                                                >  1:1202339/1‑100 (MQ=255)
                                                 cTGGACGTAAAAACATACAGAGATGATGCGCGATTAAATGAATGGAGATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACg                                                >  1:1632740/1‑100 (MQ=255)
                                                   ggACGTAAAAACATACAGAGATGATGCGCGATTAAATGAATGGAGATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGtt                                              >  1:391077/1‑100 (MQ=255)
                                                   ggACGTAAAAACATACAGAGATGATGCGCGATTAAATGAATGGAGATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGtt                                              >  2:673018/1‑100 (MQ=255)
                                                   ggACGTAAAAACATACAGAGATGATGCGCGATTAAATGAATGGAGATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGtt                                              >  2:2289657/1‑100 (MQ=255)
                                                   ggACGTAAAAACATACAGAGATGATGCGCGATTAAATGAATGGAGATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGtt                                              >  2:2499227/1‑100 (MQ=255)
                                                   ggACGTAAAAACATACAGAGATGATGCGCGATTAAATGAATGGAGATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGtt                                              >  1:1587192/1‑100 (MQ=255)
                                                   ggACGTAAAAACATACAGAGATGATGCGCGATTAAATGAATGGAGATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGtt                                              >  1:2686908/1‑100 (MQ=255)
                                                    gACGTAAAAACATACAGAGATGATGCGCGATTAAATGAATGGAGATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTa                                             <  2:2191752/100‑1 (MQ=255)
                                                       gTAAAAACATACAGAGATGATGCGCGATTAAATGAATGGAGATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTaaaa                                          <  1:3348716/100‑1 (MQ=255)
                                                        tAAAAACATACAGAGATGATGCGCTATTAAATGAATGGAGATCTAATAGCGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAAt                                         <  1:3690972/100‑1 (MQ=255)
                                                        tAAAAACATACAGAGATGATGCGCGATTAAATGAATGGAGATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAAt                                         >  2:2465567/1‑100 (MQ=255)
                                                        tAAAAACATACAGAGATGATGCGCGATTAAATGAATGGAGATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAAt                                         >  1:3970352/1‑100 (MQ=255)
                                                        tAAAAACATACAGAGATGATGCGCGATTAAATGAATGGAGATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAAt                                         <  2:3622668/100‑1 (MQ=255)
                                                        tAAAAACATACAGAGATGATGCGCGATTAAATGAATGGAGATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAAt                                         >  1:2235911/1‑100 (MQ=255)
                                                           aaaCATACAGAGATGATGCGCGATTAAATGAATGGAGATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCg                                      >  2:2367517/1‑100 (MQ=255)
                                                            aaCATACAGAGATGATGCGCGATTAAATGAATGGAGATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGc                                     <  2:1627894/100‑1 (MQ=255)
                                                            aaCATACAGAGATGATGCGCGATTAAATGAATGGAGATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGc                                     >  2:2927386/1‑100 (MQ=255)
                                                            aaCATACAGAGATGATGCGCGATTAAATGAATGGAGATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGc                                     <  1:796604/100‑1 (MQ=255)
                                                             aCATACAGAGATGATGCGCGATTAAATGAATGGATATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCt                                    <  2:1851954/100‑1 (MQ=255)
                                                             aCATACAGAGATGATGCGCGATTAAATGAATGGAGATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCt                                    <  2:237546/100‑1 (MQ=255)
                                                               aTACAGAGATGATGCGCGATTAAATGAATGGATATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTac                                  <  1:666934/100‑1 (MQ=255)
                                                               aTACAGAGATGATGCGCGATTAAATGAATGGAGATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTac                                  >  1:2839223/1‑100 (MQ=255)
                                                               aTACAGAGATGATGCGCGATTAAATGAATGGAGATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTac                                  <  2:2735565/100‑1 (MQ=255)
                                                               aTACAGAGATGATGCGCGATTAAATGAATGGAGATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTac                                  <  2:503846/100‑1 (MQ=255)
                                                                tACAGAGATGATGCGCGATTAAATGAATGGAGATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTaca                                 >  1:784090/1‑100 (MQ=255)
                                                                 aCAGAGATGATGCGCGATTAAATGAATGGAGATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTacac                                >  2:1532064/1‑100 (MQ=255)
                                                                     agaTGATGCGCGATTAAATGAATGGAGATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCtt                            <  2:488243/100‑1 (MQ=255)
                                                                      gatgatGCGCGATTAAATGAATGGAGATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCttt                           <  1:812047/100‑1 (MQ=255)
                                                                       atgatgCGCGATTAAATGAATGGAGATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTa                          <  1:2735187/100‑1 (MQ=255)
                                                                       atgatgCGCGATTAAATGAATGGAGATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTa                          <  1:2645416/100‑1 (MQ=255)
                                                                        tgatgCGCGATTAAATGAATGGAGATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAg                         >  2:719327/1‑100 (MQ=255)
                                                                        tgatgCGCGATTAAATGAATGGAGATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAg                         >  2:3234509/1‑100 (MQ=255)
                                                                        tgatgCGCGATTAAATGAATGGAGATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAg                         >  1:1728324/1‑100 (MQ=255)
                                                                        tgatgCGCGATTAAATGAATGGAGATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAg                         <  1:1272487/100‑1 (MQ=255)
                                                                        tgatgCGCGATTAAATGAATGGAGATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAg                         <  1:3882174/100‑1 (MQ=255)
                                                                         gatgCGCGATTAAATGAATGGAGATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAGt                        <  1:3099688/100‑1 (MQ=255)
                                                                           tgCGTGATTAAATGAATGCAGCTCTAATAGAGTGAAGGAATTAGGTGACTGATATTATCCTGTAAATACGCAAGTTAAAATCCGCTCCACCCTTTAGTTc                      <  2:3441694/100‑1 (MQ=255)
                                                                           tgCGCGATTAAATGAATGGAGATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAGTTc                      <  2:165706/100‑1 (MQ=255)
                                                                            ncgcgATTAAATGAATGGAGATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAGTTCa                     >  1:3234567/2‑100 (MQ=255)
                                                                               cgATTAAATGAATGGAGATCTAATAGAGTGCAGGAATTATGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAGTTCAGaa                  <  2:1114445/100‑1 (MQ=255)
                                                                                   tAAATGAATGGATATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAGTTCAGAAAggg              <  2:3816343/100‑1 (MQ=255)
                                                                                    aaaTGAATGGAGATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAGTTCAGAAAGGGt             <  1:2039578/100‑1 (MQ=255)
                                                                                    aaaTGAATGGAGATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAGTTCAGAAAGGGt             >  1:2546626/1‑100 (MQ=255)
                                                                                     aatgaatgGAGATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAGTTCAGAAAGGGtt            <  2:686008/100‑1 (MQ=255)
                                                                                     aatgaatgGAGATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAGTTCAGAAAGGGtt            >  2:2968769/1‑100 (MQ=255)
                                                                                     aatgaatgAAGATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAGTTCAGAAAGGGtt            >  2:3829630/1‑100 (MQ=255)
                                                                                      atgaatgGAGATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAGTTCAGAAAGGGTTc           <  1:124179/100‑1 (MQ=255)
                                                                                       tgaatgGTGTTCTAATAGAGTGCAGGGATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACTCTCTTTAGTTCAGAAAGGGTTcc          <  1:1890325/100‑1 (MQ=255)
                                                                                       tgaatgGAGATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAGTTCAGAAAGGGTTcc          <  1:1079418/100‑1 (MQ=255)
                                                                                       tgaatgGAGATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAGTTCAGAAAGGGTTcc          <  1:1234641/100‑1 (MQ=255)
                                                                                       tgaatgGAGATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAGTTCAGAAAGGGTTcc          <  2:2143657/100‑1 (MQ=255)
                                                                                       tgaatgGAGATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAGTTCAGAAAGGGTTcc          <  1:467353/100‑1 (MQ=255)
                                                                                         aatgGAGATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAGTTCAGAAAGGGTTCCAc        <  2:164787/100‑1 (MQ=255)
                                                                                         aatgGAGATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAGTTCAGAAAGGGTTCCAc        >  2:2704265/1‑100 (MQ=255)
                                                                                          atgGAGATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAGTTCAGAAAGGGTTCCAcc       <  1:3366788/100‑1 (MQ=255)
                                                                                          atgGAGATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAGTTCAGAAAGGGTTCCAcc       >  2:2058252/1‑100 (MQ=255)
                                                                                          atgGAAATCTAAGAGAGTGCAGGAATTAAGTGAAAGATATTATCCTGTAAATAAGCACGTTAAAATACGCTACACTCTTTAGATCAGAAAGGGTTCCAcc       <  2:1779300/100‑1 (MQ=255)
                                                                                           tgGAGATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAGTTCAGAAAGGGTTCCACCt      <  1:1460573/100‑1 (MQ=255)
                                                                                              atatCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAGTTCAGAAAGGGTTCCACCTGCa   <  2:3666967/100‑1 (MQ=255)
                                                                                               tatCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAGTTCAGAAAGGGTTCCACCTGCaa  >  2:2608841/1‑100 (MQ=255)
                                                                                               |                                                                                                   
TATCTACCCTGTTTTTCCTGAATCACTGGTGATGTTATTGCTCTCTGTTCTGGACGTAAAAACATACAGAGATGATGCGCGATTAAATGAATGGATATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAGTTCAGAAAGGGTTCCACCTGCAA  >  NC_000913/344820‑345014

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: