Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,450,334 G→A A64T (GCT→ACT)  ytfT → putative transport system permease protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,450,3340GA97.9% 494.3 / ‑2.0 147A64T (GCT→ACT) ytfTputative transport system permease protein
Reads supporting (aligned to +/- strand):  ref base G (1/2);  new base A (83/59);  total (84/63)
Fisher's exact test for biased strand distribution p-value = 5.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.32e-01

Only 100 of 152 total aligned reads displayed.

GTCGATAGCCTGGTAGCCCCGCATTTCTGGCAGGTGGTGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAGCTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGATCTCTCCGTAGGGGCGGTGATGGCTATCGC  >  NC_000913/4450244‑4450422
                                                                                          |                                                                                        
gTCGATAGCCTGGTAGCCACGCATTTCTGGGAGGTGGTGCTTCAGGGTGGTCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAGCt                                                                                        <  1:1951284/93‑1 (MQ=255)
   gATAGCCTGGTAGCCCCGCATTTCTGGCAGGTGGTGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTccc                                                                                     >  2:3891355/1‑93 (MQ=255)
    aTAGCCTGGTAGCCCCGCATTTCTGGCAGGTGGTGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCg                                                                                    <  1:4118871/93‑1 (MQ=255)
     tAGCCTGGTAGCCCCGCATTTCTGGCAGGTGGTGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGt                                                                                   >  2:3857767/1‑93 (MQ=255)
       gCCTGGTAGCCCCGCATTTCTGGCAGGAGGTGCTTCAGGATGGGCGTTTGTTCGGTAGCACCATAGACATTCTTAACCGTGCAACTCCCGTTg                                                                                 >  1:408224/1‑93 (MQ=255)
        ccTGGTAGCCCCGCATTTCTGGCAGGTGGTGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGc                                                                                <  2:3429206/93‑1 (MQ=255)
         cTGGTAGCCCCGCATTTCTGGCAGGTGGTGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCg                                                                               >  1:2057369/1‑93 (MQ=255)
         cTGGTAGCCCCGCATTTCTGGCAGGTGGTGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTATTAACCGTGCAACTCCCGTTGCg                                                                               >  2:2796940/1‑93 (MQ=255)
         cTGGTAGCCCCGCATTTCTGGCAGGTAGTGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCg                                                                               <  1:3333401/93‑1 (MQ=255)
            gTAGCCCCGCATTTCTGGCAGGTGGTGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTa                                                                            <  2:1645111/93‑1 (MQ=255)
            gTAGCCCCGCATTTCTGGCAGGTGGTGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTa                                                                            <  2:3222643/93‑1 (MQ=255)
            gTAGCCCCGCATTTCTGGCAGGTGGTGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTa                                                                            >  1:4378813/1‑93 (MQ=255)
            gTAGCCCCGCATTTCTGGCAGGTGGTGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTa                                                                            <  1:4505180/93‑1 (MQ=255)
             tAGCCCCGCATTTCTGGCAGGTGGTGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTAc                                                                           >  1:2253993/1‑93 (MQ=255)
               gCCCCGCATTTCTGGCAGGTGGTGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTg                                                                         >  2:2433613/1‑93 (MQ=255)
                ccccGCATTTCTGGCAGGTGGTGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTgg                                                                        <  2:2266042/93‑1 (MQ=255)
                ccccGCATTTCTGGCAGGTGGTGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTgg                                                                        >  2:572978/1‑93 (MQ=255)
                   cGCATTTCTGGCAGGTGGTGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGa                                                                     >  2:2902741/1‑93 (MQ=255)
                   cGCATTTCTGGCAGGTGGTGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGa                                                                     >  2:3271120/1‑93 (MQ=255)
                    gCATTTCTGGCAGGTGGTGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGAt                                                                    >  1:474314/1‑93 (MQ=255)
                     cATTTCTGGCAGGTGGTGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGAtt                                                                   <  1:686261/93‑1 (MQ=255)
                      aTTTCTGGCAGGTGGTGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTg                                                                  >  2:1254441/1‑93 (MQ=255)
                      aTTTCTGGCAGGTGGTGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTg                                                                  <  2:267190/93‑1 (MQ=255)
                      aTTTCTGGCAGGTGGTGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTg                                                                  >  2:1043165/1‑93 (MQ=255)
                      aTTTCTGGCAGGTGGTGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTg                                                                  <  1:2250586/93‑1 (MQ=255)
                        ttCTGGCAGGTGGTGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGc                                                                <  1:3698931/93‑1 (MQ=255)
                         tCTGGCAGGTGGTGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCa                                                               <  1:3777286/93‑1 (MQ=255)
                         tCTGGCAGGTGGTGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCa                                                               >  1:2925105/1‑93 (MQ=255)
                           tGGCAGGTGGTGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATg                                                             >  1:3950734/1‑93 (MQ=255)
                           tGGCAGGTGGTGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATg                                                             >  2:1775231/1‑93 (MQ=255)
                             gCAGGTGGTGCTTCAGGATGGGCGTTTGTGCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGAc                                                           <  1:965047/93‑1 (MQ=255)
                               aGGTGGTGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGc                                                         >  1:93601/1‑93 (MQ=255)
                                ggtggtGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCt                                                        <  2:3513228/93‑1 (MQ=255)
                                ggtggtGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCt                                                        <  2:3513229/93‑1 (MQ=255)
                                ggtggtGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCt                                                        <  2:3699094/93‑1 (MQ=255)
                                ggtggtGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCt                                                        <  2:4232763/93‑1 (MQ=255)
                                 gtggtgCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTg                                                       <  2:4031119/93‑1 (MQ=255)
                                   ggtgCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGt                                                     >  1:2003498/1‑93 (MQ=255)
                                    gtgCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTg                                                    <  2:2937886/93‑1 (MQ=255)
                                      gCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGAt                                                  >  2:3205115/1‑93 (MQ=255)
                                      gCTACAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGAt                                                  <  2:2979277/93‑1 (MQ=255)
                                        ttCAGGAGGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCg                                                <  1:2852158/93‑1 (MQ=255)
                                         tCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTTCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGc                                               >  1:522146/1‑93 (MQ=255)
                                           aGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAGCTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCa                                             >  2:4692083/1‑93 (MQ=255)
                                           aGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCa                                             <  1:1775520/93‑1 (MQ=255)
                                             gATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCAcc                                           >  2:3130613/1‑93 (MQ=255)
                                                gggTGTGTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGTTTGGCATGACGCTGGTGATCGCCACCGGt                                        <  1:1212841/93‑1 (MQ=255)
                                                gggCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCggg                                        >  1:1646158/1‑92 (MQ=255)
                                                gggCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTAGCGATTGGCATGACGCTGGTGATCGCCCCCggg                                        >  1:1646157/1‑92 (MQ=255)
                                                 ggCGTTTGATAGGTAGCCCCATAGACATTCTAAACCGTACAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTg                                       <  1:4215999/93‑1 (MQ=255)
                                                  gCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTgg                                      >  2:4111654/1‑93 (MQ=255)
                                                   cGTTTGTTCTGTAGCCCCATAGACATTCTTAACCGTCCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTggg                                     <  2:2213625/93‑1 (MQ=255)
                                                     tttGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGAt                                   >  1:922759/1‑93 (MQ=255)
                                                      ttGTTCGGTAGCCCCATAGACATTCTTACCCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGAtt                                  <  2:134775/93‑1 (MQ=255)
                                                       tGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTg                                 <  1:621989/93‑1 (MQ=255)
                                                       tGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTg                                 <  2:17794/93‑1 (MQ=255)
                                                         ttCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCCGGTGATCGCCACCGGTGGGATTGAt                               >  2:222664/1‑93 (MQ=255)
                                                          tCGCTAGCCCCATAGACATTCTTAACCGTGCATCTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGAtc                              <  2:2062169/93‑1 (MQ=255)
                                                           cGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGAtct                             <  2:3697683/93‑1 (MQ=255)
                                                           cGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGAtct                             >  1:2756863/1‑93 (MQ=255)
                                                            ggTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGAtctc                            >  2:2697479/1‑93 (MQ=255)
                                                              tAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGAtctctc                          <  2:1875866/93‑1 (MQ=255)
                                                               aGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGATCTCTcc                         >  1:3732907/1‑93 (MQ=255)
                                                                gCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGATCTCTCCg                        >  2:396250/1‑93 (MQ=255)
                                                                 ccccATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCTACCGGTGGGATTGATCTCTCCGt                       <  2:2251781/93‑1 (MQ=255)
                                                                 ccccATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGATCTCTCCga                       >  2:254036/1‑92 (MQ=255)
                                                                  cccATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGATCTCTCCGTa                      <  1:3077857/93‑1 (MQ=255)
                                                                  cccATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGATCTCTCCGTa                      <  1:4395014/93‑1 (MQ=255)
                                                                  cccATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGATCTCTCCGTa                      <  1:3077858/93‑1 (MQ=255)
                                                                   ccATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGATCTCTCCGTAg                     >  1:1169057/1‑93 (MQ=255)
                                                                    cATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGATCTCTCCGTAgg                    >  1:4075690/1‑93 (MQ=255)
                                                                     aTAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGATCTCTCCGTAggg                   >  2:2798030/1‑93 (MQ=255)
                                                                     aTAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGATCTCTCCGTAggg                   >  2:4445759/1‑93 (MQ=255)
                                                                     aTAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGATCTCTCCGTAggg                   <  2:2468406/93‑1 (MQ=255)
                                                                      tAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGATCTCTCCGTAgggg                  >  2:2929873/1‑93 (MQ=255)
                                                                      tAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGATCTCTCCGTAgggg                  >  1:1802797/1‑93 (MQ=255)
                                                                       aGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGATCTCTCCGTAGGGGc                 <  1:3320821/93‑1 (MQ=255)
                                                                        gACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGATCTCTCCGTAGGGGCg                >  1:4191933/1‑93 (MQ=255)
                                                                        gACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGATCTCTCCGTAGGGGCg                >  1:4289359/1‑93 (MQ=255)
                                                                         aCATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGATCTCTCCGTAGGGGCgg               <  1:3834384/93‑1 (MQ=255)
                                                                           aTTCTTAACCGTGCAACTTCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGATCTCTCCGTAGGGGCGGTg             >  2:1879963/1‑93 (MQ=255)
                                                                           aTTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGATCTCTCCGTAGGGGCGGTg             >  2:1139184/1‑93 (MQ=255)
                                                                           aTTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGATCTCTCCGTAGGGGAGGTg             >  2:2515631/1‑93 (MQ=255)
                                                                            ttCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGATCTCTCCGTAGGGGCGGTGa            <  2:1645321/93‑1 (MQ=255)
                                                                            ttCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGATCTCTCCGTAGGGGCGGTGa            <  2:3012036/93‑1 (MQ=255)
                                                                            ttCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGATCTCTCCGTAGGGGCGGTGa            >  1:2408987/1‑93 (MQ=255)
                                                                            ttCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGATCTCTCCGTAGGGGAGGCGa            >  2:485312/1‑93 (MQ=255)
                                                                             tCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGATCTCTCCGTAGGGGCGGTGAt           >  2:3917611/1‑93 (MQ=255)
                                                                             tCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGATCTCTCCGTAGGGGCGGTGAt           >  2:1120211/1‑93 (MQ=255)
                                                                             tCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGATCTCTCCGTAGGGGCGGTGAt           >  2:1524290/1‑93 (MQ=255)
                                                                             tCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGATCTCTCCGTAGGGGCGGTGAt           <  1:4170504/93‑1 (MQ=255)
                                                                             tCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGATCTCTCCGTAGGGGCGGTGAt           <  1:251539/93‑1 (MQ=255)
                                                                               ttAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGATCTCTCCGTAGGGGCGGTGATgg         >  1:2285072/1‑93 (MQ=255)
                                                                               ttAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGATCTCTCCGTAGGGGCGGTGATgg         >  1:4503201/1‑93 (MQ=255)
                                                                                   ccGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGATCTCTCCGTAGGGGCGGTGATGGCTag     >  2:2074202/1‑92 (MQ=255)
                                                                                   ccGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGATCTCTCCGTAGGGGCGGTGATGGCTAt     <  2:2132668/93‑1 (MQ=255)
                                                                                   ccGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGATCTCTCCGTAGGGGCGGTGATGGCTAt     <  2:2344670/93‑1 (MQ=255)
                                                                                    cGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGATCTCTCCGTAGGGGCGGTGATGGCTATc    >  1:2416113/1‑93 (MQ=255)
                                                                                     gTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGATCTCTCCGTAGGGGCGGTGATGGCTATCg   <  1:4426255/93‑1 (MQ=255)
                                                                                     gTGCAACTCCCGTTGCGCTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGATCTCTCCGTAGGGGCGGTGATGGCTATCg   >  2:3125795/1‑93 (MQ=255)
                                                                                          |                                                                                        
GTCGATAGCCTGGTAGCCCCGCATTTCTGGCAGGTGGTGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAGCTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGATCTCTCCGTAGGGGCGGTGATGGCTATCGC  >  NC_000913/4450244‑4450422

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