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,0450GA83.8% 677.2 / 87.6 253P857L (CCG→CTG) narZcryptic nitrate reductase 2 alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (14/27);  new base A (108/104);  total (122/131)
Fisher's exact test for biased strand distribution p-value = 6.02e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.65e-01

Only 100 of 259 total aligned reads displayed.

GCCACAATATTCGGCGCGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCGGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAATGCTTTGGCGATACCTTT  >  NC_000913/1537955‑1538136
                                                                                          |                                                                                           
gCCACAATATTCGGCGCGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCGga                                                                                           >  2:6962476/1‑93 (MQ=255)
gCCACAATATTCGGCGCGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAga                                                                                           >  2:2792482/1‑93 (MQ=255)
gCCACAATATTCGGCGCGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAga                                                                                           >  2:5140516/1‑93 (MQ=255)
gCCACAATATTCGGCGCGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAga                                                                                           >  1:3107836/1‑93 (MQ=255)
  cacaATATTCGGCGCGGTTTTGCCTGGAATTAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAgaga                                                                                         >  2:6648293/1‑93 (MQ=255)
  cacaATATTCGGCGCGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAgaga                                                                                         >  1:3295055/1‑93 (MQ=255)
  cacaATATTCGGCGCGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAgaga                                                                                         <  1:1965149/93‑1 (MQ=255)
  cacaATATTCGGCGCGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAgaga                                                                                         <  2:5203519/93‑1 (MQ=255)
   acaATATTCGGCGCGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCagagag                                                                                        >  1:5718698/1‑93 (MQ=255)
     aaTATTCGGCGCGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTc                                                                                      >  1:6723564/1‑93 (MQ=255)
     aaTATTCGGCGCGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTc                                                                                      >  2:3810764/1‑93 (MQ=255)
      atatTCGGCGCGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCa                                                                                     <  2:1384323/93‑1 (MQ=255)
       tatTCGGCGCGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCAt                                                                                    >  2:6024797/1‑93 (MQ=255)
       tatTCGGCGCGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCAt                                                                                    >  2:3580680/1‑93 (MQ=255)
           cGGCGCGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCGGAGAGTCATGCAg                                                                                <  1:5910125/93‑1 (MQ=255)
           cGGCGCGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCGGAGAGTCATGCAg                                                                                <  2:5359179/93‑1 (MQ=255)
           cGGCGCGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAg                                                                                >  1:13508/1‑93 (MQ=255)
             gcgcGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCGGAGAGTCATGCAGTa                                                                              <  1:1327856/93‑1 (MQ=255)
               gcgGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAgg                                                                            >  2:1442600/1‑93 (MQ=255)
               gcgGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAgg                                                                            >  1:5357776/1‑93 (MQ=255)
               gcgGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAgg                                                                            <  1:3571209/93‑1 (MQ=255)
                 ggTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAgggg                                                                          <  1:2225683/93‑1 (MQ=255)
                  gTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCGGAGAGTCATGCAGTAGGGGt                                                                         <  2:2628600/93‑1 (MQ=255)
                   ttttGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGtt                                                                        <  2:343726/93‑1 (MQ=255)
                    tttGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTg                                                                       <  2:415863/93‑1 (MQ=255)
                       gCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTaa                                                                    >  2:4664276/1‑93 (MQ=255)
                         cTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATa                                                                  >  2:4616912/1‑93 (MQ=255)
                         cTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATa                                                                  <  2:6162601/93‑1 (MQ=255)
                          tGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATAc                                                                 >  2:3885396/1‑93 (MQ=255)
                          tGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATAc                                                                 >  1:5543052/1‑93 (MQ=255)
                          tGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATAc                                                                 <  2:5784067/93‑1 (MQ=255)
                           ggAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATAcc                                                                <  1:3936791/93‑1 (MQ=255)
                           ggAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATAcc                                                                <  1:574710/93‑1 (MQ=255)
                            gAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCa                                                               <  2:2607843/93‑1 (MQ=255)
                             aaTCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCAc                                                              <  1:1578797/93‑1 (MQ=255)
                              aTCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCGGAGAGTCATGCAGTAGGGGTTGTAATACCACg                                                             >  2:7010727/1‑93 (MQ=255)
                               tCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGt                                                            <  2:3912697/93‑1 (MQ=255)
                               tCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGt                                                            >  2:6544187/1‑93 (MQ=255)
                                cAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTc                                                           <  1:6307481/93‑1 (MQ=255)
                                 aGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCg                                                          <  2:5342613/93‑1 (MQ=255)
                                 aGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCg                                                          >  1:891304/1‑93 (MQ=255)
                                   aTCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATCCAGTAGGGGTTGTAATACCACGTCGGt                                                        <  2:6422186/93‑1 (MQ=255)
                                    tCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGtt                                                       <  1:2851473/93‑1 (MQ=255)
                                    tCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGtt                                                       <  1:4482582/93‑1 (MQ=255)
                                      gCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTc                                                     >  1:6533358/1‑93 (MQ=255)
                                      gCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTc                                                     <  2:2744286/93‑1 (MQ=255)
                                       cATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCGGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCt                                                    >  2:7463549/1‑93 (MQ=255)
                                       cATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCt                                                    >  2:1989965/1‑93 (MQ=255)
                                        aTTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCtt                                                   <  1:4185712/93‑1 (MQ=255)
                                        aTTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCtt                                                   >  1:5493431/1‑93 (MQ=255)
                                        aTTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCtt                                                   >  2:6697238/1‑93 (MQ=255)
                                        aTTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCtt                                                   >  1:474271/1‑93 (MQ=255)
                                         ttCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCttt                                                  >  1:3782805/1‑93 (MQ=255)
                                         ttCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCttt                                                  <  1:1676076/93‑1 (MQ=255)
                                         ttCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCttt                                                  >  1:4270704/1‑93 (MQ=255)
                                         ttCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTAATGCAGTAGGGGTTGTAATACCACGTCGGTTTCttt                                                  <  1:2584297/93‑1 (MQ=255)
                                            gCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGcc                                               >  1:2807856/1‑93 (MQ=255)
                                             ccTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCGGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCa                                              >  2:5847557/1‑93 (MQ=255)
                                             ccTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCa                                              >  1:2417216/1‑93 (MQ=255)
                                             ccTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCa                                              >  1:7086508/1‑93 (MQ=255)
                                              cTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCaa                                             >  2:1706265/1‑93 (MQ=255)
                                               tttGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAg                                            >  2:6488717/1‑93 (MQ=255)
                                                 tGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGAt                                          <  1:815642/93‑1 (MQ=255)
                                                   cgcCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCGGAGAGTCATGCAGTGGGGGTTTTAATCCCACGTCGGTTTCTTTGCCAAGATGa                                        <  1:6557362/93‑1 (MQ=255)
                                                    gcCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGAc                                       <  1:6001285/93‑1 (MQ=255)
                                                    gcCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGAc                                       >  1:6169779/1‑93 (MQ=255)
                                                     ccAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGAcc                                      >  1:3012893/1‑93 (MQ=255)
                                                       aGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCa                                    >  2:401359/1‑93 (MQ=255)
                                                       aGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCa                                    <  1:1535169/93‑1 (MQ=255)
                                                        gTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCAc                                   >  2:562000/1‑93 (MQ=255)
                                                        gTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCAc                                   >  1:3548797/1‑93 (MQ=255)
                                                        gTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCAc                                   >  1:4605885/1‑93 (MQ=255)
                                                        gTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCAc                                   >  1:5255949/1‑93 (MQ=255)
                                                         tCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACg                                  >  2:2997685/1‑93 (MQ=255)
                                                         tCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACg                                  >  1:626479/1‑93 (MQ=255)
                                                             gCACTTCACACGGCTGTGAGAGCTCTGCCGGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGcaca                              >  2:1202565/1‑93 (MQ=255)
                                                              cACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTTCCAAGATGACCCACGCacac                             >  1:1044494/1‑93 (MQ=255)
                                                                  tCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGc                         <  1:1017385/93‑1 (MQ=255)
                                                                  tCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGc                         <  2:4909890/93‑1 (MQ=255)
                                                                    acacGGCTGTGAGAGCTCTGCCGGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGa                       >  2:6846026/1‑93 (MQ=255)
                                                                      acGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGaaa                     <  2:5696027/93‑1 (MQ=255)
                                                                        ggCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAAt                   <  2:7366758/93‑1 (MQ=255)
                                                                           tgtgAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAATGCt                <  2:2261613/93‑1 (MQ=255)
                                                                           tgtgAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAATGCt                >  1:7252032/1‑93 (MQ=255)
                                                                             tgAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAATGCttt              >  1:4206097/1‑93 (MQ=255)
                                                                             tgAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAATGCttt              <  1:6753572/93‑1 (MQ=255)
                                                                              gagaGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAATGCTTTg             <  2:2956007/93‑1 (MQ=255)
                                                                               agagCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAATGCTTTgg            <  2:835751/93‑1 (MQ=255)
                                                                                gagCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAATGCTTTGGc           >  2:6395689/1‑93 (MQ=255)
                                                                                gagCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAATGCTTTGGc           >  1:4081039/1‑93 (MQ=255)
                                                                                gagCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAATGCTTTGGc           <  1:4891378/93‑1 (MQ=255)
                                                                                gagCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAATGCTTTGGc           >  2:434197/1‑93 (MQ=255)
                                                                                gagCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAATGCTTTGGc           <  2:4444558/93‑1 (MQ=255)
                                                                                  gCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAATGCTTTGGCGa         >  1:540914/1‑93 (MQ=255)
                                                                                  gCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAATGCTTTGGCGa         <  2:2348066/93‑1 (MQ=255)
                                                                                     ctGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAATGCTTTGGCGATAc      <  2:1654512/93‑1 (MQ=255)
                                                                                      tGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAATGCTTTGGCGATAcc     >  1:7381659/1‑93 (MQ=255)
                                                                                      tGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAATGCTTTGGCGATAcc     >  2:3074225/1‑93 (MQ=255)
                                                                                       gCCGGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAATGCTTTGGCGATACCt    <  2:400530/93‑1 (MQ=255)
                                                                                         cGGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAATGCTTTGGCGATACCttt  <  1:5333417/93‑1 (MQ=255)
                                                                                          |                                                                                           
GCCACAATATTCGGCGCGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCGGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAATGCTTTGGCGATACCTTT  >  NC_000913/1537955‑1538136

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