Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,279,054 A→G intergenic (+483/‑33) narK → / → narG nitrite extrusion protein/nitrate reductase 1 alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,279,0540AG99.4% 549.5 / NA 159intergenic (+483/‑33)narK/narGnitrite extrusion protein/nitrate reductase 1 alpha subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (83/75);  minor base C (1/0);  total (84/75)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 167 total aligned reads displayed.

TTTCATCCCCATCACTCTTGATCGTTATCAATTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCC  >  NC_000913/1278965‑1279143
                                                                                         |                                                                                         
tttCATCCCCATCACTCTTGATCGTTATCAATTCCCACGCTGTTTCAGAGTGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCggag                                                                                        <  2:3446501/93‑1 (MQ=255)
tttCATCCCCATCACTCTTGATCGTTATCAATTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCggag                                                                                        >  1:3690728/1‑93 (MQ=255)
tttCATCCCCATCACTCTTGATCGTTATCAATTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCggag                                                                                        >  1:2975897/1‑93 (MQ=255)
tttCATCCCCATCACTCTTGATCGTTATCAATTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCggag                                                                                        <  2:4501985/93‑1 (MQ=255)
tttCATCCCCATCACTCTTGATCGTTATCAATTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCggag                                                                                        >  2:4540850/1‑93 (MQ=255)
 ttCATCCCCATCACTCTTGATCGTTATCAATTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCGGAgg                                                                                       >  2:4351086/1‑93 (MQ=255)
   cATCCCCATCACTCTTGATCGTTATCAATTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCGGAGGGc                                                                                     <  2:1138567/93‑1 (MQ=255)
   cATCCCCATCACTCTTGATCGTTATCAATTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCGGAGGGc                                                                                     >  2:4885108/1‑93 (MQ=255)
   cATCCCCATCACTCTTGATCGTTATCAATTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCCGAGGGc                                                                                     >  2:4920179/1‑93 (MQ=255)
   cATCCCCATCACTCTTGATCGTTATCAATTCCCACGCTGTTTCAGAGCGTTACCTTGCACTTAAACATTAGCAATGTCGATTTATCGGAggtc                                                                                     <  1:6150367/93‑3 (MQ=255)
    aTCCCCATCACTCTTGATCGTTATCAATTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCGGAGGGcc                                                                                    >  2:5976802/1‑93 (MQ=255)
     tCCCCATCACTCTTGATCGTTATCAATTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCGGAGGGCCg                                                                                   >  2:5932398/1‑93 (MQ=255)
     tCCCCATCACTCTTGATCGTTATCAATTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCGGAGGGCCg                                                                                   >  2:5830185/1‑93 (MQ=255)
     tCCCCATCACTCTTGATCGTTATCAATTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCGGAGGGCCg                                                                                   <  2:4061429/93‑1 (MQ=255)
     tCCCCATCACTCTTGATCGTTATCAATTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCGGAGGGCCg                                                                                   >  2:6366556/1‑93 (MQ=255)
     tCCCCATCACTCTTGATCGTTATCAATTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCGGAGGGCCg                                                                                   >  2:1145124/1‑93 (MQ=255)
       cccATCACTCTTGATCGTTATCAATTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCGGAGGGCCGAc                                                                                 <  1:4407951/93‑1 (MQ=255)
       cccATCACTCTTGATCGTTATCAATTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCGGAGGGCCGAc                                                                                 >  1:2155143/1‑93 (MQ=255)
       cccATCACTCTTGATCGTTATCAATTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCGGAGGGCCGAc                                                                                 >  2:780901/1‑93 (MQ=255)
        ccATCACTCTTGATCGTTATCAATTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCGGAGGGCCGACa                                                                                >  2:1542076/1‑93 (MQ=255)
            cACTCTTGATCGTTATCAATTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCGGAGGGCCGACAGGCt                                                                            >  2:2879920/1‑93 (MQ=255)
            cACTCTTGATCGTTATCAATTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCGGAGGGCCGACAGGCt                                                                            <  1:924483/93‑1 (MQ=255)
            cACTCTTGATCGTTATCAATTCCCACGCTGTTTCAGAGCGTTACCTTGACCTTAAACATTAGCAATGTCGATTTATCGGAGGGCCGACAGGCt                                                                            >  2:6233473/1‑93 (MQ=255)
              ctctTGATCGTTATCAATTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCGGAGGGCCGACAGGCTcc                                                                          >  1:455250/1‑93 (MQ=255)
               tctTGATCGTTATCAATTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCGGAGGGCCGACAGGCTccc                                                                         >  2:1360010/1‑93 (MQ=255)
                  tGATCGTTATCAATTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCGGAGGGCCGACAGGCTCCcaca                                                                      <  1:2698512/93‑1 (MQ=255)
                  tGATCGTTATCAATTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCGGAGGGCCGACAGGCTCCcaca                                                                      <  1:5006967/93‑1 (MQ=255)
                  tGATCGTTATCAATTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCGGAGGGCCGACAGGCTCCcaca                                                                      <  1:6441987/93‑1 (MQ=255)
                  tGATCGTTATCAATTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCGGAGGGCCGACAGGCTCCcaca                                                                      >  1:3690504/1‑93 (MQ=255)
                   gATCGTTATCAATTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCGGAGGGCCGACAGGCTCCCACAg                                                                     >  1:4206834/1‑93 (MQ=255)
                    aTCGTTATCAATTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCGGAGGGCCGACAGGCTCCCACAgg                                                                    >  2:675389/1‑93 (MQ=255)
                     tCGTTATCAATTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCGGAGGGCCGACAGGCTCCCACAGga                                                                   <  2:2550791/93‑1 (MQ=255)
                     tCGTTATCAATTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCGGAGGGCCGACAGGCTCCCACAGga                                                                   >  2:5067328/1‑93 (MQ=255)
                       gTTATCAATTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCGGAGGGCCGACAGGCTCCCACAGgaga                                                                 >  1:257650/1‑93 (MQ=255)
                        ttATCAATTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCGGAGGGCCGACAGGCTCCCACAGGAGaa                                                                <  1:4431156/93‑1 (MQ=255)
                        ttATCAATTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCGGAGGGCCGACAGGCTCCCACAGGAGaa                                                                >  1:2866714/1‑93 (MQ=255)
                        ttATCAATTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCGGAGGGCCGACAGGCTCCCACAGGAGaa                                                                >  1:5408854/1‑93 (MQ=255)
                         tATCAATTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCGGAGGGCCGACAGGCTCCCACAGGAGaaa                                                               >  1:1349551/1‑93 (MQ=255)
                          aTCAATTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCGGAGGGCCGACAGGCTCCCACAGGAGaaaa                                                              <  2:196719/93‑1 (MQ=255)
                          aTCAATTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCGGAGGGCCGACAGGCTCCCACAGGAGaaaa                                                              >  1:6046883/1‑93 (MQ=255)
                               ttCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCGGAGGGCCGACAGGCTCCCACAGGAGAAAACCGAt                                                         <  2:311932/93‑1 (MQ=255)
                                  ccACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCGGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAg                                                      <  2:5840934/93‑1 (MQ=255)
                                  ccACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCGGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAg                                                      >  2:3048286/1‑93 (MQ=255)
                                    aCGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCGGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTa                                                    >  2:4650895/1‑93 (MQ=255)
                                    aCGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCGGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTa                                                    <  1:5364477/93‑1 (MQ=255)
                                        tGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCGGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAAtt                                                >  1:6171701/1‑93 (MQ=255)
                                         gTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCGGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTc                                               >  2:2669273/1‑93 (MQ=255)
                                             cAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCGGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTgg                                           >  2:5091187/1‑93 (MQ=255)
                                              agagCGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCGGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGa                                          >  1:442340/1‑93 (MQ=255)
                                              agagCGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCGGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGa                                          <  2:4204909/93‑1 (MQ=255)
                                                 gCGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCGGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCg                                       >  2:5028768/1‑93 (MQ=255)
                                                 gCGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCGGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCg                                       >  1:4816113/1‑93 (MQ=255)
                                                   gTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCGGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGt                                     <  2:4453721/93‑1 (MQ=255)
                                                    ttACCTTGCCCTTAAACATTAGCAATGTCGATTTATCGGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGtt                                    >  2:5336219/1‑93 (MQ=255)
                                                     tACCTTGCCCTTAAACATTAGCAATGTCGATTTATCGGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGttt                                   >  1:356195/1‑93 (MQ=255)
                                                     tACCTTGCCCTTAAACATTAGCAATGTCGATTTATCGGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGgtt                                   >  1:6470236/1‑93 (MQ=255)
                                                      aCCTTGCCCTTAAACATTAGCAATGTCGATTTATCGGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTc                                  <  2:2996071/93‑1 (MQ=255)
                                                      aCCTTGCCCTTAAACATTAGCAATGTCGATTTATCGGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTc                                  <  1:3515679/93‑1 (MQ=255)
                                                       ccTTGCCCTTAAACATTAGCAATGTCGATTTATCGGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCg                                 >  2:1602526/1‑93 (MQ=255)
                                                        cTTGCCCTTAAACATTAGCAATGTCGATTTATCGGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGc                                >  2:1703960/1‑93 (MQ=255)
                                                        cTTGCCCTTAAACATTAGCAATGTCGATTTATCGGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGc                                <  1:1586777/93‑1 (MQ=255)
                                                        cTTGCCCTTAAACATTAGCAATGTCGATTTATCGGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGc                                <  1:1684619/93‑1 (MQ=255)
                                                         ttGCCCTTAAACATTAGCAATGTCGATTTATCGGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCt                               <  1:5713529/93‑1 (MQ=255)
                                                         ttGCCCTTAAACATTAGCAATGTCGATTTATCGGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCt                               >  2:391548/1‑93 (MQ=255)
                                                         ttGCCCTTAAACATTAGCAATGTCGATTTATCGGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCt                               <  1:2455353/93‑1 (MQ=255)
                                                             ccTTAAACATTAGCAATGTCGATTTATCGGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACtt                           <  2:5745002/93‑1 (MQ=255)
                                                             ccTTAAACATTAGCAATGTCGATTTATCGGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACtt                           <  1:5427843/93‑1 (MQ=255)
                                                              cTTAAACATTAGCAATGTCGATTTATCGGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTc                          >  2:5710524/1‑93 (MQ=255)
                                                              cTTAAACATTAGCAATGTCGATTTATCGGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTc                          <  2:3223671/93‑1 (MQ=255)
                                                               ttAAACATTAGCAATGTCGATTTATCGGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCa                         >  1:2021603/1‑93 (MQ=255)
                                                                 aaaCATTAGCAATGTCGATTTATCGGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCaaa                       <  1:5769473/93‑1 (MQ=255)
                                                                  aaCATTAGCAATGTCGATTTATCGGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAAc                      >  2:2342428/1‑93 (MQ=255)
                                                                   aCATTAGGAATGTCGATTTATCGGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACa                     >  1:6360133/1‑93 (MQ=255)
                                                                   aCATTAGCAATGTCGATTTATCGGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACa                     >  2:385995/1‑93 (MQ=255)
                                                                     aTTAGCAATGTCGATTTATCGGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGa                   >  2:2423060/1‑93 (MQ=255)
                                                                     aTTAGCAATGTCGATTTATCGGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGa                   >  1:3850472/1‑93 (MQ=255)
                                                                     aTTAGCAATGTCGATTTACCGGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGa                   <  2:4764559/93‑1 (MQ=255)
                                                                      ttAGCAATGTCGATTTATCGGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGaa                  >  1:2570026/1‑93 (MQ=255)
                                                                       tAGCAATGTCGATTTATCGGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAg                 <  1:5582166/93‑1 (MQ=255)
                                                                        aGCAATGTCGATTTATCGGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAgg                <  1:4905472/93‑1 (MQ=255)
                                                                        aGCAATGTCGATTTATCGGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAgg                <  1:2043648/93‑1 (MQ=255)
                                                                         gCAATGTCGATTTATCGGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAggg               <  2:1808583/93‑1 (MQ=255)
                                                                         gCAATGTCGATTTATCGGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAggg               >  2:3640685/1‑93 (MQ=255)
                                                                          cAATGTCGATTTATCGGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGt              <  2:3021616/93‑1 (MQ=255)
                                                                          cAATGTCGATTTATCGGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGt              <  2:3345646/93‑1 (MQ=255)
                                                                           aaTGTCGATTTATCGGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTg             >  2:2679766/1‑93 (MQ=255)
                                                                             tGTCGATTTATCGGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGaa           <  2:6471279/93‑1 (MQ=255)
                                                                              gTCGATTTATCGGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGaaa          <  1:3627653/93‑1 (MQ=255)
                                                                              gTCGATTTATCGGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGaaa          >  2:1857037/1‑93 (MQ=255)
                                                                                cGATTTATCGGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAAcc        <  1:4041951/93‑1 (MQ=255)
                                                                                cGATTTATCGGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAAcc        <  2:5604595/93‑1 (MQ=255)
                                                                                cGATTTATCGGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAAcc        >  1:24014/1‑93 (MQ=255)
                                                                                 gATTTATCGGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCt       >  1:2355484/1‑93 (MQ=255)
                                                                                 gATTTATCGGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCt       <  1:5239554/93‑1 (MQ=255)
                                                                                  aTTTATCGGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCtt      <  1:1588069/93‑1 (MQ=255)
                                                                                    ttATCGGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTg    <  1:3330899/93‑1 (MQ=255)
                                                                                    ttATCGGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTg    <  1:1196870/93‑1 (MQ=255)
                                                                                     tATCGGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGc   >  2:3100764/1‑93 (MQ=255)
                                                                                     tATCGGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGc   >  1:1047264/1‑93 (MQ=255)
                                                                                      aTCGGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGcc  <  1:5309256/93‑1 (MQ=255)
                                                                                         |                                                                                         
TTTCATCCCCATCACTCTTGATCGTTATCAATTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCC  >  NC_000913/1278965‑1279143

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: