Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,279,780 T→C W232R (TGG→CGG)  narG → nitrate reductase 1 alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,279,7800TC96.9% 293.6 / ‑0.1 98W232R (TGG→CGG) narGnitrate reductase 1 alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (1/2);  new base C (44/51);  total (45/53)
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.49e-01

Only 100 of 102 total aligned reads displayed.

TCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGACGCG  >  NC_000913/1279690‑1279868
                                                                                          |                                                                                        
tCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCCgg                                                                                        <  2:799013/93‑1 (MQ=255)
 cGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTggg                                                                                       >  2:4006983/1‑93 (MQ=255)
 cGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCCggg                                                                                       <  2:4396661/93‑1 (MQ=255)
   ggTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCCGGGGc                                                                                     <  2:1338023/93‑1 (MQ=255)
   ggTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCCGGGGc                                                                                     >  2:1526656/1‑93 (MQ=255)
     tGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCCGGGGCGa                                                                                   <  2:2879900/93‑1 (MQ=255)
     tGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCCGGGGCGa                                                                                   <  2:4733589/93‑1 (MQ=255)
        aCGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCCGGGGCGAGCa                                                                                >  2:2762494/1‑93 (MQ=255)
          gCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCCGGGGCGAGCaaa                                                                              <  1:2763543/93‑1 (MQ=255)
          gCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCCGGGGCGAGCaaa                                                                              <  1:2763544/93‑1 (MQ=255)
           cTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCCGGGGCGAGCaaat                                                                             >  2:1059540/1‑92 (MQ=255)
           cTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCCGGGGCGAGCAAAc                                                                             <  2:4645624/93‑1 (MQ=255)
           cTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCCGGGGCGAGCAAAc                                                                             <  2:4052151/93‑1 (MQ=255)
            tATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCCGGGGCGAGCAAACt                                                                            >  2:1578897/1‑93 (MQ=255)
            tATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCCGGGGCGAGCAAACt                                                                            >  2:1559451/1‑93 (MQ=255)
            tATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCCGGGGCGAGCAAACt                                                                            >  1:3748475/1‑93 (MQ=255)
            tATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCCGGGGCGAGCAAACt                                                                            <  1:280077/93‑1 (MQ=255)
             aTCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCCGGGGCGAGCAAACTg                                                                           >  1:2464260/1‑93 (MQ=255)
             aTCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCCGGGGCGAGCAAACTg                                                                           >  1:3500167/1‑93 (MQ=255)
              tctctcGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCCGGGGCGAGCAAACTGa                                                                          <  1:3149096/93‑1 (MQ=255)
              tctctcGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCCGGGGCGAGCAAACTGa                                                                          >  2:3706072/1‑93 (MQ=255)
               ctctcGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCCGGGGCGAGCAAACTGAc                                                                         >  1:3989935/1‑93 (MQ=255)
               ctctcGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCCGGGGCGAGCAAACTGAc                                                                         >  1:3401238/1‑93 (MQ=255)
                tctcGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCCGGGGCGAGCAAACTGACg                                                                        >  1:1481736/1‑93 (MQ=255)
                   cGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCCGGGGCGAGCAAACTGACGTAc                                                                     >  2:3350776/1‑93 (MQ=255)
                       gATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCCGGGGCGAGCAAACTGACGTACCGGa                                                                 >  2:493072/1‑93 (MQ=255)
                        aTTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCCGGGGCGAGCAAACTGACGTACCGGaa                                                                <  1:2723332/93‑1 (MQ=255)
                        aTTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCCGGGGCGAGCAAACTGACGTACCGGaa                                                                <  1:2723217/93‑1 (MQ=255)
                         ttGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCCGGGGCGAGCAAACTGACGTACCGGAAt                                                               <  2:1952353/93‑1 (MQ=255)
                         ttGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCCGGGACGAGCAAACTGACGTACCGGAAt                                                               >  2:2587005/1‑93 (MQ=255)
                          tGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCCGGGGCGAGCAAACTGACGTACCGGAATc                                                              <  2:1156613/93‑1 (MQ=255)
                          tGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCCGGGGCGAGCAAACTGACGTACCGGAATc                                                              <  2:43523/93‑1 (MQ=255)
                             cGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCCGGGGCGAGCAAACTGACGTACCGGAATctgc                                                           <  1:1844814/93‑1 (MQ=255)
                                tACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCCGGGGCGAGCAAACTGACGTACCGGAATCTGCTGa                                                        >  2:1241127/1‑93 (MQ=255)
                                 aCTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCCGGGGCGAGCAAACTGACGTACCGGAATCTGCTGAc                                                       <  1:305724/93‑1 (MQ=255)
                                 aCTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCCGGGGCGAGCAAACTGACGTACCGGAATCTGCTGAc                                                       <  1:3415014/93‑1 (MQ=255)
                                  cTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCCGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACt                                                      <  2:4659501/93‑1 (MQ=255)
                                     gCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCTTCTCCGCAAACCCGGGGCGAGCAAACTGACGTACCGGAACCTCCTGACTgat                                                   >  1:2693237/1‑91 (MQ=255)
                                     gCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCCGGGGCGACCAAACTGACGTACCGGAATCTGTGGACTGGt                                                   >  1:2693238/1‑93 (MQ=25)
                                      cTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCCGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTa                                                  >  2:2840338/1‑93 (MQ=255)
                                       ttAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCCGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTAc                                                 >  1:3670382/1‑93 (MQ=255)
                                         taGCTTCTACGAGTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACaa                                               <  2:3678029/92‑1 (MQ=255)
                                         aaGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCCGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACaa                                               >  1:2516141/1‑93 (MQ=255)
                                            cTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCCGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTc                                            >  1:1440737/1‑93 (MQ=255)
                                              tCTACTACTGTTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCCGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCa                                          <  1:531221/93‑1 (MQ=255)
                                              tCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCCGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCa                                          <  2:929844/93‑1 (MQ=255)
                                                tACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCCGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGc                                        >  1:2761419/1‑93 (MQ=255)
                                                   gACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCCGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTaa                                     >  1:384714/1‑92 (MQ=255)
                                                   gACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCCGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTAc                                     >  1:384713/1‑93 (MQ=255)
                                                   caCTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCCGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTAc                                     <  2:1294601/92‑1 (MQ=255)
                                                    gcTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTAca                                    <  1:4034931/92‑1 (MQ=255)
                                                    aCTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCCGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTAca                                    >  2:309504/1‑93 (MQ=255)
                                                     cTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCCGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTAcat                                   >  1:792190/1‑93 (MQ=255)
                                                     cTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCCGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTAcat                                   >  2:1595003/1‑93 (MQ=255)
                                                       ggTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCCGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTAcatca                                 >  2:2288857/1‑93 (MQ=255)
                                                        gTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCCGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTAcatcat                                <  1:1282981/93‑1 (MQ=255)
                                                        gTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCCGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTAcatcat                                <  2:1656977/93‑1 (MQ=255)
                                                         tACTGCGACTTGCCTCCTGCGTCTCCGCAAACCCGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACatcatc                               <  1:331256/93‑1 (MQ=255)
                                                         tACTGCGACTTGCCTCCTGCGTCTCCGCAAACCCGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACatcatc                               <  1:374966/93‑1 (MQ=255)
                                                           cTGCGACTTGCCTCCTGCGTCTCCGCAAACCCGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGc                             <  2:3728753/93‑1 (MQ=255)
                                                              cGACTTGCCTCCTGCGTCTCCGCAAACCCGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTg                          >  2:1122675/1‑93 (MQ=255)
                                                              cGACTTGCCTCCTGCGTCTCCGCAAACCCGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTg                          >  1:3418097/1‑93 (MQ=255)
                                                              cGACTTGCCTCCTGCGTCTCCGCAAACCCGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTg                          <  1:322251/93‑1 (MQ=255)
                                                               gACTTGCCTCCTGCGTCTCCGCAAACCCGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTgg                         <  1:4362772/93‑1 (MQ=255)
                                                               gACTTGCCTCCTGCGTCTCCGCAAACCCGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTgg                         <  2:4316475/93‑1 (MQ=255)
                                                                 cTTGCCTCCTGCGTCTCCGCAAACCCGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTgggg                       <  1:269745/93‑1 (MQ=255)
                                                                 cTTGCCTCCTGCGTCTCCGCAAACCCGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTgggg                       <  2:2896188/93‑1 (MQ=255)
                                                                  ttGCCTCCTGCGTCTCCGCAAACCCGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTggggg                      <  2:3887977/93‑1 (MQ=255)
                                                                  ttGCCTCCTGCGTCTCCGCAAACCCGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTggggg                      >  2:2027346/1‑93 (MQ=255)
                                                                  ctGCCTCCTGCGTCTCCGCAAACCCGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTggggg                      <  1:3584268/92‑1 (MQ=255)
                                                                   tGCCTCCTGCGTCTCCGCAAACCCGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGt                     <  2:1564743/93‑1 (MQ=255)
                                                                    gCCTCCTGCGTCTCCGCAAACCCGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTc                    <  2:634964/93‑1 (MQ=255)
                                                                     cctcctGCGTCTCCGCAAACCCGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCa                   <  1:2598694/93‑1 (MQ=255)
                                                                      ctcctGCGTCTCCGCAAACCCGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCaa                  >  2:1218815/1‑93 (MQ=255)
                                                                      ctcctGCGTCTCCGCAAACCCGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCaa                  >  2:2715878/1‑93 (MQ=255)
                                                                       tcctGCGTCTCCGCAAACCCGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCaaa                 >  2:489936/1‑93 (MQ=255)
                                                                        cctGCGTCTCCGCAAACCCGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAAc                <  1:1417720/93‑1 (MQ=255)
                                                                        cctGCGTCTCCGCAAACCCGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAAc                <  1:1931903/93‑1 (MQ=255)
                                                                         ctGCGTCTCCGCAGACCCGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTGCAACTCCAGCTACATCATCGCCTGGGGGTCAAACg               <  2:2514634/93‑1 (MQ=255)
                                                                          tGCGTCTCCGCAAACCCGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACgg              >  1:4124464/1‑92 (MQ=255)
                                                                          tGCGTCTCCGCAAACCCGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGt              >  2:2982873/1‑93 (MQ=255)
                                                                          tGCGTCTCCGCAAACCCGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGt              <  2:197828/93‑1 (MQ=255)
                                                                             gTCTCCGCAACCCCGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGcc           >  2:3509888/1‑93 (MQ=255)
                                                                             gTCTCCGCAAACCCGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGcc           <  2:123796/93‑1 (MQ=255)
                                                                             gTCTCCGCAAACCCGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGcc           >  2:3509887/1‑93 (MQ=255)
                                                                               ctcCGCAAACCCGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCga         >  2:3882417/1‑92 (MQ=255)
                                                                               ctcCGCAAACCCGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGc         <  2:2288773/93‑1 (MQ=255)
                                                                                tcCGCAAACCCGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCa        <  2:2865373/93‑1 (MQ=255)
                                                                                  cGCAAACCCGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGa      <  2:4496932/93‑1 (MQ=255)
                                                                                  cGCAAACCCGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGa      <  1:2506173/93‑1 (MQ=255)
                                                                                   gCAAACCCGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGAc     <  2:4486334/93‑1 (MQ=255)
                                                                                   gCAAACCCGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGAc     <  1:2457733/93‑1 (MQ=255)
                                                                                   ccAAACCCGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGAc     <  2:429153/92‑1 (MQ=255)
                                                                                    cAAACCCGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGAcg    <  1:1587796/93‑1 (MQ=255)
                                                                                     aaaCCCGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGAcgc   >  2:2576147/1‑93 (MQ=255)
                                                                                     aaaCCCGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATAATCGCCTGGGGGTCAAACGCGCCGCAGAcgc   >  2:1524293/1‑93 (MQ=255)
                                                                                     aaaCCCGGGGCGAGCAAACTGACGTACCAGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGAcgc   >  2:1228277/1‑93 (MQ=255)
                                                                                      aaCCCGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGAcgcg  >  1:4379769/1‑93 (MQ=255)
                                                                                      aaCCCGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGAcgcg  <  2:1421346/93‑1 (MQ=255)
                                                                                      aaCCCGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGAcgcg  >  2:854582/1‑93 (MQ=255)
                                                                                          |                                                                                        
TCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGACGCG  >  NC_000913/1279690‑1279868

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