Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,484,843 T→C V353A (GTC→GCC)  yjgP → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,484,8430TC97.0% 497.1 / 1.6 166V353A (GTC→GCC) yjgPhypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (2/3);  new base C (70/91);  total (72/94)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 175 total aligned reads displayed.

GGTAAGCTGGACCCGACGCTGTGGATGTGGACCGTTAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGTCCGCCGCCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCTATATCGGTAAAACTATTTT  >  NC_000913/4484752‑4484932
                                                                                           |                                                                                         
ggTAAGCTGGACCCGACGCTGTGGATGTGGACCGTTAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGTc                                                                                          >  2:923460/1‑93 (MQ=255)
   aaGCTGGACCCGACGCTGTGGATGTGGACCGTTAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCccgc                                                                                       >  1:3180298/1‑93 (MQ=255)
     gCTGGACCCGACGCTGTGGATGTGGACCGTTAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCccgccc                                                                                     >  2:1695551/1‑92 (MQ=255)
     gCTGGACCCGACGCTGTGGATGTGGACCGCTAACCTGATTTATCTGGCTTTAGCGATTGGTCTCAACCCACGGGACACCGCGCCGGCccgcca                                                                                     >  1:4224560/1‑92 (MQ=255)
        gatcccGACGCTGTGGATGTGGACCGTTAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCt                                                                                  <  2:4260879/90‑1 (MQ=255)
         gACCCGACGCTGTGGATGTGGACCGTTAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTg                                                                                 <  1:1770973/93‑1 (MQ=255)
           cccGACGCTGTGGATGTGGACCGTTAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTgcg                                                                               <  1:4548976/93‑1 (MQ=255)
            ccGACGCTGTGGATGTGGACCGTTAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTgcgc                                                                              <  2:3505417/93‑1 (MQ=255)
              gACGCTGTGGATGTGGACCGTTAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTgcgcgc                                                                            >  1:1760005/1‑93 (MQ=255)
                cGCTGTGGATGTGGACCGTTAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCa                                                                          >  1:4080397/1‑93 (MQ=255)
                 gCTGTGGATGTGGACCGTTAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAg                                                                         <  2:2428034/93‑1 (MQ=255)
                 gCTGTGGATGTGGACCGTTAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGAAACCGTGCCGGCCCGCCACCTGCGCGCCAg                                                                         >  1:3902499/1‑93 (MQ=255)
                   tgtggatgtggaCCGTTAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGtt                                                                       <  2:1623719/93‑1 (MQ=255)
                   tgtggatgtggaCCGTTAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGtt                                                                       >  1:4283089/1‑93 (MQ=255)
                     tggatgtggaCCGTTAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGtttt                                                                     >  1:4250158/1‑93 (MQ=255)
                       gatgttgaCTGTTACCCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCGCCTGCGCGCCAGTTTTTc                                                                   <  2:3301124/93‑1 (MQ=255)
                          gtggaCCGTTAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTcgcg                                                                <  1:4467049/93‑1 (MQ=255)
                           tggaCCGTTAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGt                                                               >  1:1257686/1‑93 (MQ=255)
                            ggaCCGTTAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTa                                                              >  1:3161689/1‑93 (MQ=255)
                            ggaCCGTTAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTa                                                              <  2:812721/93‑1 (MQ=255)
                            ggaCCGTTAACCTGATTTATCTGACTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTa                                                              <  2:2579454/93‑1 (MQ=255)
                             gaCCGTTAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTaa                                                             >  2:4178620/1‑93 (MQ=255)
                              aCCGTTAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTaaa                                                            <  1:3585498/93‑1 (MQ=255)
                              aCCGTTAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTaaa                                                            >  1:2358944/1‑93 (MQ=255)
                               ccGTTAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAg                                                           <  1:297649/93‑1 (MQ=255)
                               ccGTTAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAg                                                           >  1:709082/1‑93 (MQ=255)
                               ccGTTAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCACGCCACCTGCGCGCCAGTTTTTCGCGTAAAg                                                           >  1:7157/1‑93 (MQ=255)
                                 gTTAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGa                                                         <  1:1958086/93‑1 (MQ=255)
                                 gTTAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGa                                                         <  2:2644445/93‑1 (MQ=255)
                                 gTTAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGa                                                         >  2:2538973/1‑93 (MQ=255)
                                     acatGATTTATCTGGCTATAGCGATTGTTCTCACCCTTTGGGACACCGTGCCGGCACGCCACCTGCGAGCCAGTTTTTCGCGTAAAGGAGCgg                                                     <  2:4223528/90‑1 (MQ=255)
                                     aCCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCgg                                                     <  2:592895/93‑1 (MQ=255)
                                      ccTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGgt                                                    >  2:4602229/1‑93 (MQ=255)
                                      ccTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACATGCGCGCCAGTTTTTCGCGTAAAGGAGCGgt                                                    >  1:1864882/1‑93 (MQ=255)
                                           tttATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGAt                                               <  1:4352573/93‑1 (MQ=255)
                                            ttATCTGGCTTTAGCGATTGTTCTCATCCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATg                                              <  1:2145580/93‑1 (MQ=255)
                                              aTCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCa                                            >  1:2057319/1‑93 (MQ=255)
                                               tCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCaa                                           >  1:1004960/1‑93 (MQ=255)
                                               tCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCaa                                           <  1:3279459/93‑1 (MQ=255)
                                                 tGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAAcc                                         >  1:3006994/1‑93 (MQ=255)
                                                  ggCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCt                                        <  2:1113313/93‑1 (MQ=255)
                                                   gCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCtt                                       <  2:1522953/93‑1 (MQ=255)
                                                   gCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCtt                                       <  1:3741948/93‑1 (MQ=255)
                                                   gCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCtt                                       <  1:3362677/93‑1 (MQ=255)
                                                    cTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCttt                                      <  1:3776639/93‑1 (MQ=255)
                                                    cTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCttt                                      >  1:829325/1‑93 (MQ=255)
                                                    cTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCttt                                      <  1:1975539/93‑1 (MQ=255)
                                                     tttAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCtttt                                     <  2:2771140/93‑1 (MQ=255)
                                                     tttAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCtttt                                     <  2:1857107/93‑1 (MQ=255)
                                                     tttAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCtttt                                     <  1:861984/93‑1 (MQ=255)
                                                     tttAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCtttt                                     <  1:731312/93‑1 (MQ=255)
                                                       tAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTgg                                   <  2:2003402/93‑1 (MQ=255)
                                                        aGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGc                                  <  2:3450395/93‑1 (MQ=255)
                                                           gATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTa                               <  1:535503/93‑1 (MQ=255)
                                                           gATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTa                               <  2:397788/93‑1 (MQ=255)
                                                            aTTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTAc                              <  1:4000963/93‑1 (MQ=255)
                                                             ttGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCACCCTTTTGGCGTACt                             >  2:832106/1‑93 (MQ=255)
                                                             ttGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACt                             >  1:2608996/1‑93 (MQ=255)
                                                              tGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCCCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACtt                            >  2:2343091/1‑93 (MQ=255)
                                                               tttCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTg                           <  2:1524575/92‑1 (MQ=255)
                                                               gTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTg                           >  2:3738458/1‑93 (MQ=255)
                                                                ttCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGa                          <  2:727681/93‑1 (MQ=255)
                                                                 tctcAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGAc                         <  2:3041289/93‑1 (MQ=255)
                                                                 tctcAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGAc                         >  1:1031800/1‑93 (MQ=255)
                                                                 tctcAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGAc                         <  1:3928105/93‑1 (MQ=255)
                                                                  ctcAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGAcc                        <  1:4661925/93‑1 (MQ=255)
                                                                  ctcAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGAcc                        >  1:3949671/1‑93 (MQ=255)
                                                                   tcAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCg                       >  1:2700222/1‑93 (MQ=255)
                                                                   tcAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCg                       >  1:1449933/1‑93 (MQ=255)
                                                                     aaCCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCt                     <  1:2378830/93‑1 (MQ=255)
                                                                     aaCCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCt                     <  2:3912871/93‑1 (MQ=255)
                                                                     aaCCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCt                     <  1:3086658/93‑1 (MQ=255)
                                                                     aaCCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCt                     >  2:2179123/1‑93 (MQ=255)
                                                                      aCCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCta                    <  1:1441679/93‑1 (MQ=255)
                                                                      aCCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCta                    <  1:2519297/93‑1 (MQ=255)
                                                                      aCCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCta                    <  2:1291850/93‑1 (MQ=255)
                                                                      aCCTTTGGGACACCGTGCCGGCCCGCCACATGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCta                    <  2:817522/93‑1 (MQ=255)
                                                                       ccTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTTACCGCtat                   <  2:246631/93‑1 (MQ=255)
                                                                        cTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCtata                  >  1:951304/1‑93 (MQ=255)
                                                                         tttGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCTatat                 <  2:4158253/93‑1 (MQ=255)
                                                                          ttGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCTATATc                <  2:3537522/93‑1 (MQ=255)
                                                                           tGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCTATATCg               <  2:2304136/93‑1 (MQ=255)
                                                                             ggACACCGTGCCGGCCCGCCGCCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCTATATCGGt             <  1:2793741/93‑1 (MQ=255)
                                                                             ggACACCGTGCCGGCCCGCCCCCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGACCTTGCCCGCTCAATCgta             >  1:1205773/1‑91 (MQ=255)
                                                                             ggACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCTATATCGGt             <  1:4130297/93‑1 (MQ=255)
                                                                             ggACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTGTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCTATATCGGt             <  1:3743453/93‑1 (MQ=255)
                                                                              gACACCGTGCCGGCCCGCCGCCTGCGAGCCAGTTAGTCTCGAAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCTATATCGGTa            <  2:3020722/93‑1 (MQ=255)
                                                                               acacCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCTATATCGGTaa           >  1:3007470/1‑93 (MQ=255)
                                                                                cacCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCTATATCGGTaaa          >  2:3615634/1‑93 (MQ=255)
                                                                                 acCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCTATATCGGTaaaa         >  1:2827698/1‑93 (MQ=255)
                                                                                 acCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCTATATCGGTaaaa         >  1:4387282/1‑93 (MQ=255)
                                                                                 acCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCTATATCGGTaaaa         <  1:1378771/93‑1 (MQ=255)
                                                                                 acCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTGCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCTATATCGGTaaaa         <  1:4569115/93‑1 (MQ=255)
                                                                                    gTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCTATATCGGTAAAACTa      <  2:407875/93‑1 (MQ=255)
                                                                                    gTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCTATATCGGTAAAACTa      >  1:27636/1‑93 (MQ=255)
                                                                                    gTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCTATATCGGTAAAACTa      <  2:4542236/93‑1 (MQ=255)
                                                                                    gTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCTATATCGGTAAAACTa      >  2:1825550/1‑93 (MQ=255)
                                                                                      gCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCTATATCGGTAAAACTAtt    <  2:3838843/93‑1 (MQ=255)
                                                                                       gcGGCCCGCCACCTGCTCGACAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCTATATCGGTAAAACTAttt   <  1:1646864/92‑1 (MQ=255)
                                                                                        cGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCTATATCGGTAACACTAtttc  >  1:4415566/1‑92 (MQ=255)
                                                                                           |                                                                                         
GGTAAGCTGGACCCGACGCTGTGGATGTGGACCGTTAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGTCCGCCGCCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCTATATCGGTAAAACTATTTT  >  NC_000913/4484752‑4484932

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