Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,538,045 G→A P857L (CCG→CTG)  narZ ← cryptic nitrate reductase 2 alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,538,0450GA99.4% 1101.0 / ‑3.7 310P857L (CCG→CTG) narZcryptic nitrate reductase 2 alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (154/154);  total (155/155)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.07e-01

Only 100 of 324 total aligned reads displayed.

GCCACAATATTCGGCGCGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCGGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAATGCTTTGGCGATACCT  >  NC_000913/1537955‑1538134
                                                                                          |                                                                                         
gCCACAATATTCGGCGCGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAga                                                                                         <  1:4579892/93‑1 (MQ=255)
 ccACAATATTCGGCGCGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCagag                                                                                        >  2:615001/1‑93 (MQ=255)
  cacaATATTCGGCGCGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAgaga                                                                                       >  1:8158402/1‑93 (MQ=255)
  cacaATATTCGGCGCGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAgaga                                                                                       >  2:1246914/1‑93 (MQ=255)
   acaATATTCGGCGCGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCagagag                                                                                      >  2:1128109/1‑93 (MQ=255)
      atatTCGGCGCGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCa                                                                                   >  2:2101612/1‑93 (MQ=255)
       tatTCGGCGCGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGTGCACTTCACACGGCTGTGAGGGATTTGCCAGAGAGTAAt                                                                                  >  2:4484178/1‑93 (MQ=255)
       tatTCGGCGCGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCAt                                                                                  >  1:3894122/1‑93 (MQ=255)
         ttCGGCGCGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGc                                                                                <  1:22789/93‑1 (MQ=255)
            tgcgcGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGCGTCATGCAGt                                                                             >  2:4291163/2‑93 (MQ=255)
              cgcgGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAg                                                                           >  1:8678812/1‑93 (MQ=255)
              cgcgGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAg                                                                           <  2:5054014/93‑1 (MQ=255)
              cgcgGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAg                                                                           <  1:5860390/93‑1 (MQ=255)
               gcgGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAgg                                                                          <  2:4084977/93‑1 (MQ=255)
               gcgGGTTTGCCTGGAATCAGATCCCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAgg                                                                          <  1:5262459/93‑1 (MQ=255)
                cgGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAggg                                                                         >  2:3215224/1‑93 (MQ=255)
                cgGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAggg                                                                         >  2:8564213/1‑93 (MQ=255)
                  gTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGt                                                                       >  1:6039622/1‑93 (MQ=255)
                  gTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGt                                                                       <  1:8234778/93‑1 (MQ=255)
                       gCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTaa                                                                  <  1:4368788/93‑1 (MQ=255)
                         cTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATa                                                                >  1:1003206/1‑93 (MQ=255)
                          tGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATAc                                                               >  1:4051211/1‑93 (MQ=255)
                          tGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATAc                                                               >  1:7379449/1‑93 (MQ=255)
                          tGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGACATGCAGTAGGGGTTGTAATAc                                                               >  1:3564564/1‑93 (MQ=255)
                            gAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCa                                                             >  1:2029527/1‑93 (MQ=255)
                            gAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCa                                                             >  2:2359554/1‑93 (MQ=255)
                            gAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGGTGTAATACCa                                                             >  1:5937344/1‑93 (MQ=255)
                             aaTCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCAc                                                            <  2:1903096/93‑1 (MQ=255)
                             aaTCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCAc                                                            <  1:1489720/93‑1 (MQ=255)
                             aaTCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCAc                                                            <  2:2999757/93‑1 (MQ=255)
                             aaTCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCAc                                                            <  2:4029597/93‑1 (MQ=255)
                              aTCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACg                                                           >  1:5112767/1‑93 (MQ=255)
                              aTCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACg                                                           <  2:7425058/93‑1 (MQ=255)
                               tCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGt                                                          <  2:1802343/93‑1 (MQ=255)
                                cAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGTTTGTAATACCACgcg                                                         >  1:1224898/1‑91 (MQ=255)
                                cAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGAGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAGTACCACGTc                                                         <  2:2723956/93‑1 (MQ=255)
                                 aGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCg                                                        >  2:8435210/1‑93 (MQ=255)
                                  gATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCgg                                                       >  1:2555555/1‑93 (MQ=255)
                                  gATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCgg                                                       <  2:7114404/93‑1 (MQ=255)
                                  gATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCgg                                                       >  2:1548458/1‑93 (MQ=255)
                                   aTCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGt                                                      >  1:311554/1‑93 (MQ=255)
                                   aTCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGt                                                      >  2:6914237/1‑93 (MQ=255)
                                    tCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGtt                                                     <  2:5333493/93‑1 (MQ=255)
                                    tCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGtt                                                     <  2:1488296/93‑1 (MQ=255)
                                      gCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTc                                                   >  2:6265456/1‑93 (MQ=255)
                                         ttCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCttt                                                <  1:6893810/93‑1 (MQ=255)
                                         ttCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCttt                                                <  2:6853383/93‑1 (MQ=255)
                                         ttCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCttt                                                <  1:7705484/93‑1 (MQ=255)
                                         ttCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTCTCttt                                                >  1:4784400/1‑93 (MQ=255)
                                            gCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGcc                                             >  1:5341951/1‑93 (MQ=255)
                                            gCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGcc                                             <  2:5421568/93‑1 (MQ=255)
                                            gCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGcc                                             <  1:1652869/93‑1 (MQ=255)
                                            cccTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGcc                                             <  2:4227199/92‑1 (MQ=255)
                                             ccTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCa                                            >  2:5569941/1‑93 (MQ=255)
                                             ccTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCa                                            >  1:7481100/1‑93 (MQ=255)
                                             ccTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCa                                            >  2:773472/1‑93 (MQ=255)
                                              cTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCaa                                           >  1:5695316/1‑93 (MQ=255)
                                               tttGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTCCCAAg                                          >  1:1848810/1‑93 (MQ=255)
                                                ttGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGa                                         >  1:1113865/1‑93 (MQ=255)
                                                ttGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGa                                         <  2:2492997/93‑1 (MQ=255)
                                                   cgcCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACTACGTCGGTTTCTTTGCCAAGATgc                                      >  1:4948946/1‑92 (MQ=255)
                                                   cgcCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGa                                      <  1:2163143/93‑1 (MQ=255)
                                                     ccAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGAcc                                    <  2:7924561/93‑1 (MQ=255)
                                                     ccAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGAAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGAcc                                    >  1:3072171/1‑93 (MQ=255)
                                                         tCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACg                                <  2:8405132/93‑1 (MQ=255)
                                                         tCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACg                                <  2:5094053/93‑1 (MQ=255)
                                                         tCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATAACACGTCGGTTTCTTTGCCAAGATGACCCACg                                >  2:7958509/1‑93 (MQ=255)
                                                            aGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGcac                             >  2:16986/1‑93 (MQ=255)
                                                                cTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGCCCCACGCACACtt                         >  2:751772/1‑93 (MQ=255)
                                                                cTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACtt                         >  1:3318258/1‑93 (MQ=255)
                                                                cTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACtt                         >  1:7521922/1‑93 (MQ=255)
                                                                 ttCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTg                        <  2:4808861/93‑1 (MQ=255)
                                                                   cacaCGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCg                      <  2:6681612/93‑1 (MQ=255)
                                                                     cacGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGaa                    >  2:8645832/1‑93 (MQ=255)
                                                                     cacGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGaa                    <  2:100160/93‑1 (MQ=255)
                                                                       cGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGaaaa                  >  2:8200416/1‑93 (MQ=255)
                                                                       cGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGaaaa                  >  1:1295171/1‑93 (MQ=255)
                                                                        ggCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAAt                 <  2:344741/93‑1 (MQ=255)
                                                                        ggCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAAt                 <  2:8114096/93‑1 (MQ=255)
                                                                           tgtgAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAATGCt              <  1:709223/93‑1 (MQ=255)
                                                                           tgtgAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAATGCt              >  2:3697047/1‑93 (MQ=255)
                                                                            gtgAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAATGcct             >  2:4682970/1‑91 (MQ=255)
                                                                            gtgAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAATGCtt             <  1:7201163/93‑1 (MQ=255)
                                                                             tgAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAATGCttt            <  1:6223400/93‑1 (MQ=255)
                                                                             tgAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAATGCttt            <  2:5801173/93‑1 (MQ=255)
                                                                             tgAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAATGCttt            <  1:6721597/93‑1 (MQ=255)
                                                                             tgAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAATGCttt            <  2:245952/93‑1 (MQ=255)
                                                                             tgAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAATGCttt            <  2:7688533/93‑1 (MQ=255)
                                                                              gagaGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAATGCTTTg           <  1:6532423/93‑1 (MQ=255)
                                                                              gagaGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAATGCTTTg           <  1:2685439/93‑1 (MQ=255)
                                                                               agagCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAATGCTTTgg          <  2:2908694/93‑1 (MQ=255)
                                                                                gagCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAATGCTTTGGc         >  2:3814935/1‑93 (MQ=255)
                                                                                gagCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAATGCTTTGGc         >  1:4449858/1‑93 (MQ=255)
                                                                                 agCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAATGCTTTGGCg        >  2:389551/1‑93 (MQ=255)
                                                                                 agCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAATGCTTTGGCg        >  2:2590740/1‑93 (MQ=255)
                                                                                  gCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAATGCTTTGGCGa       <  1:1355347/93‑1 (MQ=255)
                                                                                   ctctGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAATGCTTTGGCGAt      >  1:3414272/1‑93 (MQ=255)
                                                                                    tctGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAATGCTTTGGCGATa     >  1:3820051/1‑93 (MQ=255)
                                                                                     ctGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAATGCTTTGGCGATAc    <  1:5489951/93‑1 (MQ=255)
                                                                                       gCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAATGCTTTGGCGATACCt  <  2:3653526/93‑1 (MQ=255)
                                                                                          |                                                                                         
GCCACAATATTCGGCGCGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCGGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAATGCTTTGGCGATACCT  >  NC_000913/1537955‑1538134

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