Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,008,990 C→T R109C (CGT→TGT)  yigM → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,008,9900CT99.5% 1387.3 / ‑5.0 382R109C (CGT→TGT) yigMhypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base T (173/207);  total (175/207)
Fisher's exact test for biased strand distribution p-value = 2.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.42e-01

Only 100 of 387 total aligned reads displayed.

TATCTCTACCTGACGGTTTCCGAGCTGCTGCTGTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCCGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATCAGGTCACCGACCATTTCTGG  >  NC_000913/4008900‑4009079
                                                                                          |                                                                                         
 aTCTCTACCTGACGGTTTCCGAGCTGCTGCTGTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTc                                                                                        >  2:7674116/1‑93 (MQ=255)
 aTCTCTACCTGACGGTTTCCGAGCTGCTGCTGTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTc                                                                                        >  2:2041039/1‑93 (MQ=255)
 aTCTCTACCTGACGGTTTCCGAGCTGCTGCTGTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTc                                                                                        <  1:2940561/93‑1 (MQ=255)
  tctcTACCTGACGGTTTCCGAGCTGCTGCTGTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGCGTAAACGCTGTCt                                                                                       >  2:2237551/1‑93 (MQ=255)
  tctcTACCTGACGGTTTCCGAGCTGCTGCTGTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCt                                                                                       >  1:7773175/1‑93 (MQ=255)
   ctctACCTGACGGTTTCCGAGCTGCTGCTGTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTg                                                                                      <  1:5549011/93‑1 (MQ=255)
    tctACCTGACGGTTGCCGAGCTGCTGCTGTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTgc                                                                                     <  1:1620462/93‑1 (MQ=255)
      tACCTGACGGTTTCCGAGCTGCTGCTGTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTgcgc                                                                                   <  2:4061069/93‑1 (MQ=255)
      tACCTGACGGTTTCCGAGCTGCTGCTGTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTgcgc                                                                                   >  2:3964420/1‑93 (MQ=255)
        ccTGACGGTTTCCGAGCTGCTGCTGTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTg                                                                                 >  1:5682097/1‑93 (MQ=255)
        ccTGACGGTTTCCGAGCTGCTGCTGTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTg                                                                                 <  1:2745240/93‑1 (MQ=255)
         cTGACGGTTTCCGAGCTGCTGCTGTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTgg                                                                                >  1:6188101/1‑93 (MQ=255)
         cTGACGGTTTCCGAGCTGCTGCTGTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTgg                                                                                >  1:3038819/1‑93 (MQ=255)
          tGACGGTTTCCGAGCTGCTGCTGTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTggg                                                                               <  1:7291277/93‑1 (MQ=255)
           gACGGTTTCCGAGCTGCTGCTGTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTgggg                                                                              <  2:5184767/93‑1 (MQ=255)
            aCGGTTTCCGAGCTGCTGCTGTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGc                                                                             <  1:2303821/93‑1 (MQ=255)
              ggTTTCCGAGCTGCTGCTGTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTa                                                                           <  1:7919084/93‑1 (MQ=255)
              ggTTTCCGAGCTGCTGCTGTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTa                                                                           <  2:6397464/93‑1 (MQ=255)
               gTTTCCGAGCTGCTGCTGTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTAt                                                                          <  1:6863454/93‑1 (MQ=255)
               gTTTCCGAGCTGCTGCTGTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTAt                                                                          <  2:4399341/93‑1 (MQ=255)
                  tCCGAGCTGCTGCTGTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGcc                                                                       <  1:4418354/93‑1 (MQ=255)
                    cGAGCTGCTGCTGTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCtt                                                                     >  2:6864262/1‑93 (MQ=255)
                     gAGCTGCTGCTGTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCttt                                                                    >  2:3796320/1‑93 (MQ=255)
                     gAGCTGCTGCTGTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCttt                                                                    >  1:4972385/1‑93 (MQ=255)
                      aGCTGCTGCTGTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTa                                                                   <  2:1101776/93‑1 (MQ=255)
                       gctgctgctGTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAg                                                                  <  2:1804564/93‑1 (MQ=255)
                       gctgctgctGTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAg                                                                  >  2:461004/1‑93 (MQ=255)
                       gctgctgctGTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAg                                                                  <  1:2259910/93‑1 (MQ=255)
                        ctgctgctgTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAgc                                                                 <  2:13260/93‑1 (MQ=255)
                         tgctgctgTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAgcg                                                                >  1:522351/1‑93 (MQ=255)
                          gctgctgTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAgcgc                                                               >  1:3186785/1‑93 (MQ=255)
                          gctgctgTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAgcgc                                                               >  2:7840938/1‑93 (MQ=255)
                           ctgctgTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGcc                                                              <  2:6442044/93‑1 (MQ=255)
                            tgctgTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCt                                                             <  2:5613213/93‑1 (MQ=255)
                             gctgTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCtt                                                            <  1:7053430/93‑1 (MQ=255)
                             gctgTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCtt                                                            <  1:7640158/93‑1 (MQ=255)
                             gctgTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCACTGGGGCTATGCCTTTAGCGCCtt                                                            <  2:1410097/93‑1 (MQ=255)
                              ctgTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTg                                                           >  2:6014313/1‑93 (MQ=255)
                                gTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCt                                                         >  1:109620/1‑93 (MQ=255)
                                 ttCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTg                                                        >  1:4649431/1‑93 (MQ=255)
                                  tCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTgg                                                       <  1:1945716/93‑1 (MQ=255)
                                  tCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTgg                                                       <  1:2476650/93‑1 (MQ=255)
                                  tCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTgg                                                       <  1:3834109/93‑1 (MQ=255)
                                  tCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTgg                                                       <  2:7318465/93‑1 (MQ=255)
                                      cgcgCTGACGTCGCTCTACACCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGt                                                   <  1:7109802/90‑1 (MQ=255)
                                        tGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGa                                                 <  2:3495816/93‑1 (MQ=255)
                                        tGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGa                                                 <  1:1189948/93‑1 (MQ=255)
                                        tGCTGACGCCGCTCTACATCACGCTGATTTATAACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGa                                                 >  1:6929041/1‑93 (MQ=255)
                                         gCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGAt                                                <  2:3317176/93‑1 (MQ=255)
                                         gCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGAt                                                <  1:2636849/93‑1 (MQ=255)
                                         gCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGAt                                                <  2:7379664/93‑1 (MQ=255)
                                          gtGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGAtt                                               <  1:1505128/92‑1 (MQ=255)
                                           tGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTg                                              >  2:2516395/1‑93 (MQ=255)
                                           tGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTg                                              <  1:1001661/93‑1 (MQ=255)
                                           tGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTg                                              <  2:2746018/93‑1 (MQ=255)
                                           tGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTg                                              <  2:4998033/93‑1 (MQ=255)
                                           tGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTg                                              <  1:966831/93‑1 (MQ=255)
                                            gACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTgg                                             <  2:7731973/93‑1 (MQ=255)
                                            gACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTgg                                             <  2:5676238/93‑1 (MQ=255)
                                                  gcTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGCTTGGTGCCgg                                       <  2:5094272/93‑1 (MQ=255)
                                                    tctACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGa                                     >  1:3663552/1‑93 (MQ=255)
                                                    tctACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGa                                     >  2:4708418/1‑93 (MQ=255)
                                                     ctACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGat                                    <  2:3830604/93‑1 (MQ=255)
                                                     ctACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGat                                    >  2:3283457/1‑93 (MQ=255)
                                                       aCATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGatta                                  <  2:3690191/93‑1 (MQ=255)
                                                        cATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGTGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGattat                                 <  2:846946/93‑1 (MQ=255)
                                                        cATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGattat                                 <  2:7484062/93‑1 (MQ=255)
                                                         aTCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGattatt                                >  1:3291103/1‑93 (MQ=255)
                                                          tCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTc                               <  2:3130500/93‑1 (MQ=255)
                                                            aCGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGc                             <  2:2929209/93‑1 (MQ=255)
                                                              gCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTa                           >  2:7133153/1‑93 (MQ=255)
                                                              gCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTa                           <  2:5391503/93‑1 (MQ=255)
                                                              gCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTa                           <  1:6599213/93‑1 (MQ=255)
                                                               cTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTAt                          >  1:6407663/1‑93 (MQ=255)
                                                                tGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATg                         >  1:6477407/1‑93 (MQ=255)
                                                                   tttATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATc                      <  1:5233242/93‑1 (MQ=255)
                                                                   tttATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATc                      >  1:6042268/1‑93 (MQ=255)
                                                                   tttATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATc                      >  2:2068808/1‑93 (MQ=255)
                                                                   tttATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATc                      <  2:4298063/93‑1 (MQ=255)
                                                                    ttATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATCa                     <  1:1838829/93‑1 (MQ=255)
                                                                    ttATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATCa                     <  1:5409862/93‑1 (MQ=255)
                                                                    ttATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATCa                     <  2:2348187/93‑1 (MQ=255)
                                                                    ttATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATCa                     <  2:6672581/93‑1 (MQ=255)
                                                                     tATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATCAg                    <  1:4311134/93‑1 (MQ=255)
                                                                     tATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATCAg                    >  1:2808942/1‑93 (MQ=255)
                                                                     tATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATCAg                    <  2:7576712/93‑1 (MQ=255)
                                                                      aTGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATCAgg                   <  1:3472473/93‑1 (MQ=255)
                                                                      aTGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATCAgg                   >  2:2799898/1‑93 (MQ=255)
                                                                        gACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATCAGGTc                 <  2:12696/93‑1 (MQ=255)
                                                                           atcatGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATCAGGTCAcc              >  1:3986649/1‑93 (MQ=255)
                                                                           atcatGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATCAGGTCAcc              <  2:255595/93‑1 (MQ=255)
                                                                              atGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATCAGGTCACCGAc           >  1:4851454/1‑93 (MQ=255)
                                                                               tGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATCAGGTCACCGAcc          <  1:3396945/93‑1 (MQ=255)
                                                                               tGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATCAGGTCACCGAcc          >  1:5916929/1‑93 (MQ=255)
                                                                                  gTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATCAGGTCACCGACCAtt       <  2:586400/93‑1 (MQ=255)
                                                                                  gTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTAGCTATGATCAGGTCACCGACCAtt       <  1:3132319/93‑1 (MQ=255)
                                                                                   tAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATCAGGTCACCGACCAttt      >  1:4480807/1‑93 (MQ=255)
                                                                                   tAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATCAGGTCACCGACCAttt      >  1:7553572/1‑93 (MQ=255)
                                                                                      gcgcTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATCAGGTCACCGACCATTTCTg   >  1:4365280/1‑93 (MQ=255)
                                                                                       cgcTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATCAGGTCACCGACCATTTCTgg  <  1:828016/93‑1 (MQ=255)
                                                                                          |                                                                                         
TATCTCTACCTGACGGTTTCCGAGCTGCTGCTGTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCCGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATCAGGTCACCGACCATTTCTGG  >  NC_000913/4008900‑4009079

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