Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,299,266 T→C N503S (AAC→AGC)  napA ← putative nitrate reductase 3

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,299,2660TC98.3% 554.3 / ‑2.2 179N503S (AAC→AGC) napAputative nitrate reductase 3
Reads supporting (aligned to +/- strand):  ref base T (2/1);  new base C (79/97);  total (81/98)
Fisher's exact test for biased strand distribution p-value = 5.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

Only 100 of 193 total aligned reads displayed.

CCCTCTTTCTCTACCCACATTGCGGTCGGCAGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGTTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACATTC  >  NC_000913/2299176‑2299355
                                                                                          |                                                                                         
cccTCTTTCTCTACCCACATTGCGGTCGGCAGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTg                                                                                         >  2:6310407/1‑93 (MQ=255)
 ccTCTTTCTCTACCCACATTGCGGTCGGCAGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTgg                                                                                        >  2:3296217/1‑92 (MQ=255)
  ctctTTCTCTACCCACATTGCGGTCGGCAGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTgcg                                                                                       >  1:8419207/1‑93 (MQ=255)
   tctTTCTCTACCCACATTGCGGTCGGCAGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTgcgc                                                                                      <  2:1378059/93‑1 (MQ=255)
      ttCTCTACCCACATTGCGGTCGGCAGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGa                                                                                   <  1:7029310/93‑1 (MQ=255)
        ctctACCCACATTGCGGTCGGCAGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGGGGAtg                                                                                 >  2:2126954/1‑92 (MQ=255)
        ctctACCCACATTGCGGTCGGCAGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATc                                                                                 <  1:96996/93‑1 (MQ=255)
         tctACCCACATTGCGGTCGGCAGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATcg                                                                                >  1:9159919/1‑93 (MQ=255)
         tctACCCACATTGCGGTCGGCAGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATcg                                                                                <  1:8710828/93‑1 (MQ=255)
           tACCCACATTGCGGTCGGCAGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATcgcg                                                                              >  2:6800938/1‑93 (MQ=255)
           tACCCACATTGCGGTCGGCAGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGACACGGATCGGAGGCGAAGATGAAGCTGCGCGGATcgcc                                                                              >  1:384338/1‑92 (MQ=255)
            aCCCACATTGCGGTCGGCAGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCgcgc                                                                             <  2:8034583/93‑1 (MQ=255)
             cccACATTGCGGTCGGCAGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGcc                                                                            <  1:6100990/93‑1 (MQ=255)
             cccACATTGCGGTCGGCAGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGcc                                                                            <  1:5828938/93‑1 (MQ=255)
              ccACATTGCGGTCGGCAGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCa                                                                           >  1:8141726/1‑93 (MQ=255)
               cacaTTGCGGTCGGCAGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAg                                                                          <  1:6419417/93‑1 (MQ=255)
               cacaTTGCGGTCGGCAGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAg                                                                          <  2:2691345/93‑1 (MQ=255)
               aacaTTGCGGTCGGCAGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAg                                                                          >  2:100034/2‑93 (MQ=255)
                  aTTGCGGTCGGCAGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGccc                                                                       <  1:3936753/93‑1 (MQ=255)
                  aTTGCGGTCGGCAGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGccc                                                                       <  1:3625884/93‑1 (MQ=255)
                  aTTGCGGTCGGCAGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGccc                                                                       >  1:4543227/1‑93 (MQ=255)
                   ttGCGGTCGGCAGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCg                                                                      <  1:8156397/93‑1 (MQ=255)
                      cGGTCGGCAGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCa                                                                   <  2:9012007/93‑1 (MQ=255)
                      cGGTCGGCAGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCa                                                                   <  2:5294059/93‑1 (MQ=255)
                      cGGTCGGCAGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCa                                                                   >  1:535949/1‑93 (MQ=255)
                       ggTCGGCAGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCAt                                                                  <  2:1786740/93‑1 (MQ=255)
                       ggTCGGCAGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCAt                                                                  <  1:7792736/93‑1 (MQ=255)
                        gTCGGCAGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATa                                                                 <  2:4665565/93‑1 (MQ=255)
                        gTCGGCAGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATa                                                                 >  2:8215627/1‑93 (MQ=255)
                            gCAGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGAACGCGCCAGCCCGGGATACGct                                                             >  2:2452660/1‑93 (MQ=255)
                             cAGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCCGCCCGGCATACGctc                                                            >  2:5228122/1‑93 (MQ=255)
                             cAGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGctc                                                            <  2:9806433/93‑1 (MQ=255)
                             cAGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGctc                                                            >  2:5778072/1‑93 (MQ=255)
                             cAGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGccc                                                            >  2:5841114/1‑91 (MQ=255)
                             cAGGATCAAGTCGGCGGCCAGCGCACGGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGctc                                                            <  1:7388075/93‑1 (MQ=255)
                              aGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGctct                                                           <  2:8397029/93‑1 (MQ=255)
                              aGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGctct                                                           >  2:2554546/1‑93 (MQ=255)
                               ggATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCtt                                                          >  1:4880191/1‑93 (MQ=255)
                                gATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTc                                                         >  1:649402/1‑93 (MQ=255)
                                gATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTc                                                         >  2:7501260/1‑93 (MQ=255)
                                   cAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCAtt                                                      >  1:9096500/1‑93 (MQ=255)
                                   cAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACACTCTTCAtt                                                      <  1:4776131/93‑1 (MQ=255)
                                    aaGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATATGCTCTTCATTg                                                     >  2:8942506/1‑93 (MQ=255)
                                    aaGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTg                                                     <  1:2630080/93‑1 (MQ=255)
                                       tCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATATGCTCTTCATTGATg                                                  >  1:8526209/1‑93 (MQ=255)
                                       tCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATg                                                  >  1:3670856/1‑93 (MQ=255)
                                        cggcggCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGt                                                 <  2:5486780/93‑1 (MQ=255)
                                        cggcggCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGt                                                 <  1:2820686/93‑1 (MQ=255)
                                        cggcggCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGt                                                 <  1:9522270/93‑1 (MQ=255)
                                         ggcggcCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGtt                                                <  1:2903880/93‑1 (MQ=255)
                                         ggcggcCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGtt                                                <  2:4138052/93‑1 (MQ=255)
                                          gcggcCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGTTGCGCGGATCGCGCCAGCCCGGCCTACGCTCTTCATTGATTTTc                                               >  1:180227/1‑93 (MQ=255)
                                          gcggcCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTc                                               <  2:4940754/93‑1 (MQ=255)
                                          gcggcCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTc                                               <  1:5579967/93‑1 (MQ=255)
                                          gcggcCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTc                                               <  2:985502/93‑1 (MQ=255)
                                           cggtCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTcg                                              <  2:6776320/93‑1 (MQ=255)
                                            ggcCAGCGCGCTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTcgg                                             <  1:5866148/93‑1 (MQ=255)
                                              ccAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCCGCATACGCTCTTCATTGATGTTcggcc                                           >  1:4974538/1‑93 (MQ=255)
                                                aGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTcggcccg                                         >  1:2069674/1‑93 (MQ=255)
                                                   gcACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTcggcccggcc                                      <  2:9146743/93‑1 (MQ=255)
                                                    cACTGACTGTCGGATACGGATGGGAGTCGATGGTGAAGTTGCGCGGCTGGCGCCACGCCGGCAGCAGCGTTTCCGTGATGTTcggaccggccc                                     >  1:6431043/1‑92 (MQ=255)
                                                    cACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCt                                     <  2:4562291/93‑1 (MQ=255)
                                                     actgactgTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTg                                    >  1:4791548/1‑93 (MQ=255)
                                                       tgactgTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCa                                  <  2:652585/93‑1 (MQ=255)
                                                        gactgTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCAt                                 >  2:1726347/1‑93 (MQ=255)
                                                           tgTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCAtgtt                              >  2:850378/1‑93 (MQ=255)
                                                            gTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCAtgttg                             >  1:9448065/1‑93 (MQ=255)
                                                            gTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCAtgttg                             >  1:5505561/1‑93 (MQ=255)
                                                            gTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCAtgttg                             <  1:888993/93‑1 (MQ=255)
                                                            gTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCAtgttg                             <  1:9542553/93‑1 (MQ=255)
                                                             tCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCAtgttgt                            <  2:35207/93‑1 (MQ=255)
                                                             tCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCAtgttgt                            <  2:2897623/93‑1 (MQ=255)
                                                                gATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTgg                         >  2:7355913/1‑93 (MQ=255)
                                                                gATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTgg                         >  1:6327084/1‑93 (MQ=255)
                                                                gATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTgg                         <  1:6919092/93‑1 (MQ=255)
                                                                 aTACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGt                        >  1:1941331/1‑93 (MQ=255)
                                                                  tACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTa                       <  2:3700273/93‑1 (MQ=255)
                                                                     ggATCGGAGCCGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTacac                    <  1:3356127/93‑1 (MQ=255)
                                                                     ggATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTacac                    <  2:1127052/93‑1 (MQ=255)
                                                                      gATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTAcaca                   >  1:4448836/1‑93 (MQ=255)
                                                                           gagaCGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTc              >  1:7517017/1‑93 (MQ=255)
                                                                             gaCGATGATGAGGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCa            >  1:244878/1‑93 (MQ=255)
                                                                               cGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGt          <  1:7304804/93‑1 (MQ=255)
                                                                                gatgatGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTa         <  2:400376/93‑1 (MQ=255)
                                                                                 atgatgAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTaa        >  1:1395678/1‑93 (MQ=255)
                                                                                 atgatgAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTaa        >  2:2318403/1‑93 (MQ=255)
                                                                                 atgatgAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTaa        >  1:1585052/1‑93 (MQ=255)
                                                                                 atgatgAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTaa        >  1:2249940/1‑93 (MQ=255)
                                                                                  tgatgaAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTaaa       <  1:1158248/93‑1 (MQ=255)
                                                                                   gatgaAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAAc      <  1:1099532/93‑1 (MQ=255)
                                                                                    atgaAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACa     <  1:4245173/93‑1 (MQ=255)
                                                                                    atgaAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACa     <  1:4177432/93‑1 (MQ=255)
                                                                                    atgaAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACa     >  2:4461406/1‑93 (MQ=255)
                                                                                     tgaAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACAt    <  2:4396319/93‑1 (MQ=255)
                                                                                      taaGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACAtt   <  2:6954252/92‑1 (MQ=255)
                                                                                      gaAGTTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACAtt   <  2:6708341/93‑1 (MQ=255)
                                                                                      gaAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACAtt   <  2:2477934/93‑1 (MQ=255)
                                                                                       aaGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACATTc  <  1:1479410/93‑1 (MQ=255)
                                                                                       aaGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACATTc  >  2:6288841/1‑93 (MQ=255)
                                                                                       aaGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACATTc  <  2:507897/93‑1 (MQ=255)
                                                                                          |                                                                                         
CCCTCTTTCTCTACCCACATTGCGGTCGGCAGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGTTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACATTC  >  NC_000913/2299176‑2299355

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: