Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,008,222 A→C I156L (ATC→CTC)  rhtC → threonine efflux pump

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,008,2220AC99.5% 640.7 / NA 183I156L (ATC→CTC) rhtCthreonine efflux pump
Reads supporting (aligned to +/- strand):  ref base A (0/1);  major base C (101/80);  minor base T (1/0);  total (102/81)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 187 total aligned reads displayed.

GCTAATCCGAAAGCGATTATCTACTTTGGCTCGGTGTTCTCATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGATCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTGGTTATCAACGTCTGGCGAAGTGGATTG  >  NC_000913/4008126‑4008318
                                                                                                |                                                                                                
 cTAATCCGAAAGCGATTAACTACTTTGGCTCGGTGTTCTCATTGTTTGTCGGTGATAACATTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCAt                                                                                              >  2:3576675/1‑100 (MQ=255)
  tAATCCGAAAGCGATTATCTACTTTGGCTCGGTGTTCTCATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCAtt                                                                                             >  2:3393275/1‑100 (MQ=255)
   aaTCCGAAAGCGATTATCTACTTTGGCTCGGTGTTCTCATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTg                                                                                            >  1:2579366/1‑100 (MQ=255)
     tCCGAAAGCGATTATCTACTTTGGCTCGGTGTTCTCATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTc                                                                                          >  1:2505965/1‑100 (MQ=255)
       cGAAAGCGATTATCTACTTTGGCTCGGTGTTCTCATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGa                                                                                        >  2:1136491/1‑100 (MQ=255)
       cGAAAGCGATTATCTACTTTGGCTCGGTGTTCTCATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGa                                                                                        >  1:4242881/1‑100 (MQ=255)
         aaaGCGATTATCTACTTTGGCTCGGTGTTCTCATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGaaa                                                                                      <  1:4935975/100‑1 (MQ=255)
          aaGCGATTATCTACTTTGGCTCGGTGTTCTCATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAAc                                                                                     >  2:3682945/1‑100 (MQ=255)
           aGCGATTATCTACTTTGGCTCGGTGTTCTCATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACg                                                                                    <  1:3113407/100‑1 (MQ=255)
            gCGATTATCTACTTTGGCTCGGTGTTCTCATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGc                                                                                   >  2:4795892/1‑100 (MQ=255)
             cGATTATCTACTTTGGCTCGGTGTTCTCATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCt                                                                                  >  2:2962313/1‑100 (MQ=255)
              gATTATCTACTTTGGCTCGGTGTTCTCATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCtn                                                                                 <  1:2164053/100‑2 (MQ=255)
                ttATCTACTTTGGCTCGGTGTTCTCATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTggt                                                                               <  2:2736080/100‑2 (MQ=255)
                  aTCTACTTTGGCTCGGTGTTCTCATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGt                                                                             >  2:4185069/1‑100 (MQ=255)
                  aTCTACTTTGGCTCGGTGTTCTCATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGt                                                                             >  1:4572100/1‑100 (MQ=255)
                    cTACTTTGGCTCGGTGTTCTCATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTgg                                                                           <  2:4479861/100‑1 (MQ=255)
                     tACTTTGGCTCGGTGTTCTCATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGt                                                                          <  2:3211731/100‑1 (MQ=255)
                     tACTTTGGCTCGGTGTTCTCATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGt                                                                          <  1:840421/100‑1 (MQ=255)
                     tACTTTGGCTCGGTGTTCTCATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGt                                                                          >  2:3479949/1‑100 (MQ=255)
                      aCTTTGGCTCGGTGTTCTCATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGtt                                                                         <  1:5360380/100‑1 (MQ=255)
                        tttGGCTCGGTGTTCTCATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTa                                                                       >  1:5481402/1‑100 (MQ=255)
                          tGGCTCGGTGTTCTCATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTAcc                                                                     <  1:3784804/100‑1 (MQ=255)
                          tGGCTCGGTGTTCACATTGTTTGTCGGTGATAACGTTGGCACTCCCGCGAGCTGGGGCATCTTTGCGATGCTCATTGTCGAAGCGCTGGCGGGGTTTAcc                                                                     >  2:1341819/1‑100 (MQ=255)
                           ggCTCGTTGTTCTCATTGTTTGTCGGTGATAACGTTGGTACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACcg                                                                    <  2:4103376/100‑1 (MQ=255)
                            gCTCGGTGTTCTCATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACcgt                                                                   >  1:5526727/1‑100 (MQ=255)
                               cGGTGTTCTCATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACcgtcgt                                                                <  2:663661/100‑1 (MQ=255)
                                ggTGTTCTCATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGtt                                                               <  1:3481798/100‑1 (MQ=255)
                                ggTGTTCTCATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCCGGCGTGGTTTACCgttgtt                                                               >  2:3141888/1‑100 (MQ=255)
                                 gttttCTCATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTg                                                              <  2:246244/97‑1 (MQ=255)
                                 gtgtTCTCCGTGTTTGTCGGTGATAGCGGTGGCACTACCGCGGGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTg                                                              <  2:2425309/100‑1 (MQ=255)
                                 gtgtTCTCATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTg                                                              >  1:3667026/1‑100 (MQ=255)
                                 gtgtTCTCATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTg                                                              <  1:1479839/100‑1 (MQ=255)
                                  tgtTCTCATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGc                                                             <  2:2013980/100‑1 (MQ=255)
                                   gtTCTCATTGTTTGTCGGTGATAACGTTTGCACAACCGCGCGCTGGGTAATTTTTGCGCTGCTCACTGTCGAAACGCTGGCGTGGTTTTCCGTCGTTGcc                                                            >  2:1507125/1‑100 (MQ=255)
                                     tctcATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAg                                                          >  1:2426980/1‑100 (MQ=255)
                                      ctcATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGc                                                         >  1:4664690/1‑100 (MQ=255)
                                      ctcATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGc                                                         >  2:1393423/1‑100 (MQ=255)
                                       tcATTGTTTGTCGGTGATAACGTTGGCACTACCTCGCGCTGGGGCATTTTTGAGCTGCTCATTGTCTTAACCCTGGCGTGCTTTACCCTCGTtgcctgcc                                                        >  2:5057077/1‑100 (MQ=255)
                                       tcATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGcc                                                        >  1:3652276/1‑100 (MQ=255)
                                       tcATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGcc                                                        >  2:3276464/1‑100 (MQ=255)
                                       tcATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGcc                                                        <  1:2473113/100‑1 (MQ=255)
                                       tcATTGTTTGTCGGTGATAACGTGGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGcc                                                        <  2:4326828/100‑1 (MQ=255)
                                        cATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCt                                                       >  1:332663/1‑100 (MQ=255)
                                        cATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCt                                                       >  1:330266/1‑100 (MQ=255)
                                        cATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCt                                                       >  2:5104228/1‑100 (MQ=255)
                                        cATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCt                                                       >  1:1242084/1‑100 (MQ=255)
                                         aTTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTg                                                      >  1:1180408/1‑100 (MQ=255)
                                         aTTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTg                                                      >  1:3173348/1‑100 (MQ=255)
                                         aTTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTg                                                      <  2:4855228/100‑1 (MQ=255)
                                           tgtttgtCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGtt                                                    <  1:4399822/100‑1 (MQ=255)
                                           tgtttgtCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGtt                                                    >  1:4222668/1‑100 (MQ=255)
                                           tgtttgtCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGtt                                                    <  2:2108472/100‑1 (MQ=255)
                                           tgtttgtCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGtt                                                    >  2:1639588/1‑100 (MQ=255)
                                           tgtttgtCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGtt                                                    >  2:1749385/1‑100 (MQ=255)
                                            gtttgtCGGTGATAACGTTGGCGCTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGttt                                                   >  2:9317/1‑100 (MQ=255)
                                              ttgtCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGc                                                 >  2:2632454/1‑100 (MQ=255)
                                              ttgtCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGATCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGc                                                 <  2:1593822/100‑1 (MQ=255)
                                                gtCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGccc                                               <  1:5175257/100‑1 (MQ=255)
                                                  cGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTg                                             >  2:5457085/1‑100 (MQ=255)
                                                  cGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTg                                             >  2:2763935/1‑100 (MQ=255)
                                                   ggTGATATCGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGCCGTTGCCAGCCTGTTTGCCCTGc                                            >  2:2115902/1‑100 (MQ=255)
                                                          aCGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAAt                                     <  1:2706567/100‑1 (MQ=255)
                                                           cGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATg                                    >  2:3865203/1‑100 (MQ=255)
                                                           cGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATg                                    >  2:81567/1‑100 (MQ=255)
                                                           cGTTGCCACTACCGCGCGCTGCGGCATTTTTGCGCTGCTCATTGTCGAAACGCTCCCGAGGTTTACCGTCCTTGCCAGCCTGTTTGCCCTGCCGCACATg                                    >  2:2276232/1‑100 (MQ=255)
                                                             ttGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATgcg                                  <  2:4730481/100‑1 (MQ=255)
                                                             ttGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATgcg                                  <  1:5013082/100‑1 (MQ=255)
                                                                gCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCg                               >  1:1904575/1‑100 (MQ=255)
                                                                gCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCg                               >  2:1278458/1‑100 (MQ=255)
                                                                 cACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGt                              >  2:1961574/1‑100 (MQ=255)
                                                                 cACTACCGCGCGCTGGGGCATTTTTGCGCTACTCATTGTCAAAACACTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGt                              >  1:4839292/1‑100 (MQ=255)
                                                                  aCTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTg                             >  2:5041967/1‑100 (MQ=255)
                                                                   cgaCCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTgg                            <  2:1682858/98‑1 (MQ=255)
                                                                   cTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGTCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTgg                            >  2:5151822/1‑100 (MQ=255)
                                                                   cTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTgg                            >  2:1183960/1‑100 (MQ=255)
                                                                    tACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTGGt                           <  2:3962405/100‑1 (MQ=255)
                                                                      ccGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTGGTTa                         <  1:2178914/100‑1 (MQ=255)
                                                                      ccGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTGGTTa                         <  1:5567002/100‑1 (MQ=255)
                                                                      ccGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTGGTTa                         <  2:3921483/100‑1 (MQ=255)
                                                                      ccGCGCGCTGGGGCATTTTTGCGCTGCCCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTGGTTa                         >  1:4906710/1‑100 (MQ=255)
                                                                       cgcgcgCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTGGTTAt                        >  1:3423625/1‑100 (MQ=255)
                                                                       cgcgcgCTGGGGCAGTTTTGCGCTGCCCATTGTCGAAACGCTGGCGTGGTTTACCGTGGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTGGTTAt                        <  2:3831863/100‑1 (MQ=255)
                                                                        gcgcgcTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTGGTTATc                       <  2:1497052/100‑1 (MQ=255)
                                                                         cgcgcTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTGGTTATCa                      <  1:2701262/100‑1 (MQ=255)
                                                                            gcTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTGGTTATCAACg                   <  1:2686252/100‑1 (MQ=255)
                                                                              tGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTGGTTATCAACGTc                 <  1:418230/100‑1 (MQ=255)
                                                                                 ggCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTGGTTATCAACGTCTgg              >  1:4297321/1‑100 (MQ=255)
                                                                                  gCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTGGTTATCAACGTCTGGc             <  2:2870143/100‑1 (MQ=255)
                                                                                   cATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACGATCGTTGCCAGCCTGTTTGCCCTGCCGGAAATGCGCCGTGGTTGTCACCTTCTGGCg            <  2:4528771/100‑1 (MQ=255)
                                                                                     tttttGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTGGTTATCAACGTCTGGCGaa          >  1:181568/1‑100 (MQ=255)
                                                                                      ttttGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTGGTTATCAACGTCTGGCGaan         <  1:2366890/100‑2 (MQ=255)
                                                                                       tttGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTGGTTATCAACGTCTGGCGAAGt        >  2:3529203/1‑100 (MQ=255)
                                                                                         tGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTGGTTATCAACGTCTGGCGAAGTgg      <  1:2078051/100‑1 (MQ=255)
                                                                                          gcgcTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTGGTTATCAACGTCTGGCGAAGTGGa     <  2:4908746/100‑1 (MQ=255)
                                                                                           cgcTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTGGTTATCAACGTCTGGCGAAGTGGAt    <  1:1158003/100‑1 (MQ=255)
                                                                                           cgcTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTGGTTATCAACGTCTGGCGAAGTGGAt    >  1:3969280/1‑100 (MQ=255)
                                                                                           cgcTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTGGTTATCAACGTCTGGCGAAGTGGAt    <  2:2071902/100‑1 (MQ=255)
                                                                                            gctgctCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTGGTTATCAACGTCTGGCGAAGTGGAtt   <  2:37082/100‑1 (MQ=255)
                                                                                            gctgctCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTGGTTATCAACGTCTGGCGAAGTGGAtt   >  1:5276639/1‑100 (MQ=255)
                                                                                             cTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTGGTTATCAACGTCTGGCGAAGTGGATTg  >  2:4602787/1‑100 (MQ=255)
                                                                                                |                                                                                                
GCTAATCCGAAAGCGATTATCTACTTTGGCTCGGTGTTCTCATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGATCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTGGTTATCAACGTCTGGCGAAGTGGATTG  >  NC_000913/4008126‑4008318

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