Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,545,666 A→G D81G (GAC→GGC) 
D81G (GAC→GGC) 
b4525 →
fdnG →
hypothetical protein
formate dehydrogenase‑N, nitrate‑inducible, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,545,6660AG82.0% 520.3 / 115.5 255D81G (GAC→GGC) 
D81G (GAC→GGC) 
b4525
fdnG
hypothetical protein
formate dehydrogenase‑N, nitrate‑inducible, alpha subunit
Reads supporting (aligned to +/- strand):  ref base A (30/16);  new base G (118/91);  total (148/107)
Fisher's exact test for biased strand distribution p-value = 3.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.46e-02

Only 100 of 255 total aligned reads displayed.

CATACTGTTCCGTAGGTTGCGGGCTATTGATGTATAGCCTGGGTGATGGCGCAAAAAACGCCAGAGAAGCGATTTATCACATTGAAGGTGACCCGGATCATCCGGTAAGCCGTGGTGCGCTGTGCCCGAAAGGGGCCGGTTTGCTGGATTACGTCAACAGTGAAAACCGTCTGCGCTACCCG  >  NC_000913/1545576‑1545757
                                                                                          |                                                                                           
caTACTGTTCCGTAGGTTGCGGGCTATTGATGTATAGCCTGGGTGATGGCGCAAAAAACGCCAGAGAAGCGATTTATCACATTGAACGTGGcc                                                                                           <  2:4724144/93‑1 (MQ=255)
  tACTGTTCCGTAGGTTGCGGGCTATTGATGTATAGCCTGGGTGATGGCGCAAAAAACGCCAGAGAAGCGATTTATCACATTGAAGGTGGCCCg                                                                                         >  1:3148338/1‑93 (MQ=255)
  tACTGTTCCGTAGGTTGCGGGCTATTGATGTATAGCCTGGGTGATGGCGCAAAAAACGCCAGAGAAGCGATTTATCACATTGAAGGTGGCCCg                                                                                         >  1:1409700/1‑93 (MQ=255)
   aCTGTTCCGTAGGTTGCGGGCTATTGATGTATAGCCTGGGTGATGGCGCAAAAAACGCCAGAGAAGCGATTTATCACATTGAAGGTGGCCCgg                                                                                        >  1:6358258/1‑93 (MQ=255)
    cTGTTCCGTAGGTTGCGGGCTATTGATGTATAGCCTGGGTGATGGCGCAAAAAACGCCAGAGAAGCGATTTATCACATTGAAGGTGGCCCGGa                                                                                       >  1:6570516/1‑93 (MQ=255)
    cTGTTCCGTAGGTTGCGGGCTATTGATGTATAGCCTGGGTGATGGCGCAAAAAACGCCAGAGAAGCGATTTATCACATTGAAGGTGGCCCGGa                                                                                       <  2:3092541/93‑1 (MQ=255)
    cTGTTCCGTAGGTTGCGGGCTATTGATGTATAGCCTGGGTGATGGCGCAAAAAACGCCAGAGAAGCGATTTATCACATTGAAGGTGGCCCGGa                                                                                       >  1:1353558/1‑93 (MQ=255)
    cTGTTCCGTAGGTTGCGGGCTATTGATGTATAGCCTGGGTGATGGCGCAAAAAACGCCAGAGAAGCGATTTATCACATTGAAGGTGGCCCGGa                                                                                       <  1:6463259/93‑1 (MQ=255)
    cTGTTCCGTAGGTTGCGGGCTATTGATGTATAGCCTGGGTGATGGCGCAAAAAACGCCAGAGAAGCGATTTATCACATTGAAGGTGACCCGGa                                                                                       >  1:2794691/1‑93 (MQ=255)
     tGTTCCGTAGGTTGCGGGCTATTAATGTATAGCCTGGGTGATGGCGCAAAAAACGCCAGAGAAGCGATTTATCACATTGAAGGTGGCCCGGat                                                                                      >  1:4599913/1‑93 (MQ=255)
       ttCCGTAGGTTGCGGGCTATTGATGTATAGCCTGGGTGATGGCGCAAAAAACGCCAGAGAAGCGATTTATCACATTGAAGGTGGCTCGGatca                                                                                    <  1:6733939/93‑1 (MQ=255)
        tCCGTAGGTTGCGGGCTATTGATGTATAGCCTGGGTGATGGCGCAAAAAACGCCAGAGAAGCGATTTATCACATTGAAGGTGGCCCGGatcat                                                                                   <  2:1913658/93‑1 (MQ=255)
         ccGTAGGTTGCGGGCTATTGATGTATAGCCTGGGTGATGGCGCAAAAAACGCCAGAGAAGCGATTTATCACATTGAAGGTGGCCCGGatcatc                                                                                  <  1:1661820/93‑1 (MQ=255)
           gTAGGTTGCGGGCTATTGATGTATAGCCTGGGTGATGGCGCAAAAAACGCCAGAGAAGCGATTTATCACATTGAAGGTGGCCCGGATCATCCg                                                                                >  2:1349243/1‑93 (MQ=255)
           gTAGGTTGCGGGCTATTGATGTATAGCCTGGGTGATGGCGCAAAAAACGCCAGAGAAGCGATTTATCACATTGAAGGTGGCCCGGATCATCCg                                                                                >  1:1220184/1‑93 (MQ=255)
              ggTTGCGGGCTATTGATGTATAGCCTGGGTGATGGCGCAAAAAACGCCAGAGAAGCGATTTATCACATTGAAGGTGGCCCGGATCATCCGGTa                                                                             <  2:6409489/93‑1 (MQ=255)
               gTTGCGGGCTATTGATGTATAGCCTGGGTGATGGCGCAAAAAACGCCAGAGAAGCGATTTATCACATTGAAGGTGGCCCGGATCATCCGGTaa                                                                            >  1:3897813/1‑93 (MQ=255)
               gTTGCGGGCTATTGATGTATAGCCTGGGTGATGGCGCAAAAAACGCCAGAGAAGCGATTTATCACATTGAAGGTGACCCGGATCATCCGGTaa                                                                            <  2:4563165/93‑1 (MQ=255)
               gTTGCGGGCTATTGATGTATAGCCTGGGTGATGGCGCAAAAAACGCCAGAGAAGCGATTTATCACATTGAAGGTGACCCGGATCATCCGGTaa                                                                            <  1:3732190/93‑1 (MQ=255)
                 tGCGGGCTATTGATGTATAGCCTGGGTGATGGCGCAAAAAACGCCAGAGAAGCGATTTATCACATTGAAGGTGACCCGGATCATCCGGTAAGc                                                                          <  1:7522395/93‑1 (MQ=255)
                     ggCTATTGATGTATAGCCTGGGTGATGGCGCAAAAAACGCCAGAGAAGCGATTTATCACATTGAAGGTGGCCCGGATCATCCGGTAAGCCgtg                                                                      <  2:6432909/93‑1 (MQ=255)
                     ggCTATTGATGTATAGCCTGGGTGATGGCGCAAAAAACGCCAGAGAAGCGATTTATCACATTGAAGGTGGCCCGGATCATCCGGTAAGCCgtg                                                                      >  1:7100124/1‑93 (MQ=255)
                     ggCTATTGATGTATAGCCTGGGTGATGGCGCAAAAAACGCCAGAGAAGCGATTTATCACATTGAAGGTGGCCCGGATCATCCGGTAAGCCgtg                                                                      >  2:1986443/1‑93 (MQ=255)
                       cTATTGATGTATAGCCTGGGTGATGGCGCAAAAAACGCCAGAGAAGCGATTTATCACATTGAAGGTGGCCCGGATCATCCGGTAAGCCgtggt                                                                    <  2:5592979/93‑1 (MQ=255)
                       cTATTGATGTATAGCCTGGGTGATGGCGCAAAAAACGCCAGAGAAGCGATTTATCACATTGAAGGTGACCCGGATCATCCGGTAAGCCgtggt                                                                    >  1:5770184/1‑93 (MQ=255)
                       cTATTGATGTATAGCCTGGGTGATGGCGCAAAAAACGCCAGAGAAGCGATTTATCACATTGAAGGTGACCCGGATCATCCGGTAAGCCgtggt                                                                    >  1:2702290/1‑93 (MQ=255)
                        tATTGATGTATAGCCTGGGTGATGGCGCAAAAAACGCCAGAGAAGCGATTTATCACATTGAAGGTGGCCCGGATCATCCGGTAAGCCgtggtg                                                                   >  2:3023169/1‑93 (MQ=255)
                        tATTGATGTATAGCCTGGGTGATGGCGCAAAAAACGCCAGAGAAGCGATTTATCACATTGAAGGTGACCCGGATCATCCGGTAAGCCgtggtg                                                                   >  2:5168151/1‑93 (MQ=255)
                         aTTGATGTATAGCCTGGGTGATGGCGCAAAAAACGCCAGAGAAGCGATTTATCACATTGAAGGTGGCCCGGATCATCCGGTAAGCCGTGGTgc                                                                  >  1:526757/1‑93 (MQ=255)
                          ttGATGTATAGCCTGGGTGATGGCGCAAAAAACGCCAGAGAAGCGATTTATCACATTGAAGGTGGCCCGGATCATCCGGTAAGCCGTGGTgcg                                                                 >  1:2242516/1‑93 (MQ=255)
                           tGATGTATAGCCTGGGTGATGGCGCAAAAAACGCCAGAGAAGCGATTTATCACATTGAAGGTGGCCCGGATCATCCGGTAAGCCGTGGTgcgc                                                                <  2:294723/93‑1 (MQ=255)
                            gATGTATAGCCTGGGTGATGGCGCAAAAAACGCCAGAGAAGCGATTTATCACATTGAAGGTGGCCCGGATCATCCGGTAAGCCGTGGTGCGCt                                                               >  1:3402551/1‑93 (MQ=255)
                            gATGTATAGCCTGGGTGATGGCGCAAAAAACGCCAGAGAAGCGATTTATCACATTGAAGGTGACCCGGATCATCCGGTAAGCCGTGGTGCGCt                                                               >  2:7195993/1‑93 (MQ=255)
                             aTGTATAGCCTGGGTGATGGCGCAAAAAACGCCAGAGAAGCGATTTATCACATTGAAGGTGGCCCGGATCATCCGGTAAGCCGTGGTGCGCtg                                                              <  1:3641257/93‑1 (MQ=255)
                              tGTATAGCCTGGGTGATGGCGCAAAAAACGCCAGAGAAGCGATTTATCACATTGAAGGTGGCCCGGATCATCCGGTAAGCCGTGGTGCGCtgt                                                             <  2:5323824/93‑1 (MQ=255)
                              tGTATAGCCTGGGTGATGGCGCAAAAAACGCCAGAGAAGCGATTTATCACATTGAAGGTGGCCCGGATCATCCGGTAAGCCGTGGTGCGCtgt                                                             <  1:984766/93‑1 (MQ=255)
                              tGTATAGCCTGGGTGATGGCGCAAAAAACGCCAGAGAAGCGATTTATCACATTGAAGGTGGCCCGGATCATCCGGTAAGCCGTGGTGCGCtgt                                                             >  1:4673945/1‑93 (MQ=255)
                               gTATAGCCTGGGTGATGGCGCAAAAAACGCCAGAGAAGCGATTTATCACATTGAAGGTGGCCCGGATCATCCGGTAAGCCGTGGTGCGCtgtg                                                            >  1:1887498/1‑93 (MQ=255)
                                tataGCCTGGGTGATGGCGCAAAAAACGCCAGAGAAGCGATTTATCACATTGAAGGTGGCCCGGATCATCCGGTAAGCCGTGGTGCGCTGTGc                                                           >  2:3258696/1‑93 (MQ=255)
                                 ataGCCTGGGTGATGGCGCAAAAAACGCCAGAGAAGCGATTTATCACATTGAAGGTGGCCCGGATCATCCGGTAAGCCGTGGTGCGCTGTGcc                                                          >  2:706186/1‑93 (MQ=255)
                                  taGCCTGGGTGATGGCGCAAAAAACGCCAGAGAAGCGATTTATCACATTGAAGGTGGCCCGGATCATCCGGTAAGCCGTGGTGCGCTGTGccc                                                         >  1:5536248/1‑93 (MQ=255)
                                  taGCCTGGGTGATGGCGCAAAAAACGCCAGAGAAGCGATTTATCACATTGAAGGTGGCCCGGATCATCCGGTAAGCCGTGGTGCGCTGTGccc                                                         >  2:1845772/1‑93 (MQ=255)
                                  taGCCTGGGTGATGGCGCAAAAAACGCCAGAGAAGCGATTTATCACATTGAAGGTGGCCCGGATCATCCGGTAAGCCGTGGTGCGCTGTGccc                                                         >  2:6865484/1‑93 (MQ=255)
                                    gCCTGGGTGATGGCGCAAAAAACGCCAGAGAAGCGATTTATCACATTGAAGGTGGCCCGGATCATCCGGTAAGCCGTGGTGCGCTGTGCCCGa                                                       >  2:4432287/1‑93 (MQ=255)
                                        gggTGATGGCGCAAAAAACGCCAGAGAAGCGATTTATCACATTGAAGGTGGCCCGGATCATCCGGTAAGCCGTGGTGCGCTGTGCCCGAAAgg                                                   <  2:5098734/93‑1 (MQ=255)
                                        gggTGATGGCGCAAAAAACGCCAGAGAAGCGATTTATCACATTGAAGGTGACCCGGATCATCCGGTAAGCCGTGGTGCGCTGTGCCCGAAAgg                                                   >  1:3298633/1‑93 (MQ=255)
                                          gTGATGGCGCAAAAAACGCCAGAGAAGCGATTTATCACATTGAAGGTGGCCCGGATCATCCGGTAAGCCGTGGTGCGCTGTGCCCGAAAgggg                                                 <  2:2800550/93‑1 (MQ=255)
                                          gTGATGGCGCAAAAAACGCCAGAGAAGCGATTTATCACATTGAAGGTGACCCGGATCATCCGGTAAGCCGTGGTGCGCTGTGCCCGAAAgggg                                                 <  2:3491762/93‑1 (MQ=255)
                                           tGATGGCGCAAAAAACGCCAGAGAAGCGATTTATCACATTGAAGGTGGCCCGGATCATCCGGTAAGCCGTGGTGCGCTGTGCCCGAAAGGGGc                                                <  2:2311207/93‑1 (MQ=255)
                                           tGATGGCGCAAAAAACGCCAGAGAAGCGATTTATCACATTGAAGGTGGCCCGGATCATCCGGTAAGCCGTGGTGCGCTGTGCCCGAAAGGGGc                                                >  2:110775/1‑93 (MQ=255)
                                           tGATGGCGCAAAAAACGCCAGAGAAGCGATTTATCACATTGAAGGTGGCCCGGATCATCCGGTAAGCCGTGGTGCGCTGTGCCCGAAAGGGGc                                                >  1:3580837/1‑93 (MQ=255)
                                           tGATGGCGCAAAAAACGCCAGAGAAGCGATTTATCACATTGAAGGTGGCCCGGATCATCCGGTAAGCCGTGGTGCGCTGTGCCCGAAAGGGGc                                                <  1:6242603/93‑1 (MQ=255)
                                            gATGGCGCAAAAAACGCCAGAGAAGCGATTTATCACATTGAAGGTGGCCCGGATCATCCGGTAAGCCGTGGTGCGCTGTGCCCGAAAGGGGcc                                               <  2:6952090/93‑1 (MQ=255)
                                             aTGGCGCAAAAAACGCCAGAGAAGCGATTTATCACATTGAAGGTGGCCCGGATCATCCGGTAAGCCGTGGTGCGCTGTGCCCGAAAGGGGCCg                                              <  1:729890/93‑1 (MQ=255)
                                             aTGGCGCAAAAAACGCCAGAGAAGCGATTTATCACATTGAAGGTGGCCCGGATCATCCGGTAAGCCGTGGTGCGCTGTGCCCGAAAGGGGCCg                                              >  2:6731300/1‑93 (MQ=255)
                                                gcgcAAAAAACGCCAGAGAAGCGATTTATCACATTGAAGGTGACCCGGATCATCCGGTAAGCCGTGGTGCGCTGTGCCCGAAAGGGGCCGGtt                                           <  2:5781301/93‑1 (MQ=255)
                                                     aaaaaCGCCAGAGAAGCGATTTATCACATTGAAGGTGGCCCGGATCATCCGGTAAGCCGTGGTGCGCTGTGCCCGAAAGGGGCCGGTTTGCTg                                      <  1:1046722/93‑1 (MQ=255)
                                                     aaaaaCGCCAGAGAAGCGATTTATCACATTGAAGGTGGCCCGGATCATCCGGTAAGCCGTGGTGCGCTGTGCCCGAAAGGGGCCGGTTTGCTg                                      <  2:3887547/93‑1 (MQ=255)
                                                     aaaaaCGCCAGAGAAGCGATTTATCACATTGAAGGTGGCCCGGATCATCCGGTAAGCCGTGGTGCGCTGTGCCCGAAAGGGGCCGGTTTGCTg                                      <  1:4345686/93‑1 (MQ=255)
                                                     aaaaaCGACAGAGAAGCGATTTATCACATTGAAGGTGGCCCGGATCATCCGGTAAGCCGTGGTGCGCTGTGCCCGAAAGGGGCCGGTTTGCTg                                      <  2:6058003/93‑1 (MQ=255)
                                                        aaCGCCAGAGAAGCGATTTATCACATTGAAGGTGACCCGGATCATCCGGTAAGCCGTGGTGCGCTGTGCCCGAAAGGGGCCGGTTTGCTGGAt                                   <  2:3785394/93‑1 (MQ=255)
                                                         aCGCCAGAGAAGCGATTTATCACATTGAAGGTGACCCGGATCATCCGGTAAGCCGTGGTGCGCTGTGCCCGAAAGGGGCCGGTTTGCTGGAtt                                  >  1:3997957/1‑93 (MQ=255)
                                                           gCCAGAGAAGCGATTTATCACATTGAAGGTGGCCCGGATCATCCGGTAAGCCGTGGTGCGCTGTGCCCGAAAGGGGCCGGTTTGCTGGATTAc                                <  2:5373585/93‑1 (MQ=255)
                                                           gCCAGAGAAGCGATTTATCACATTGAAGGTGGCCCGGATCATCCGGTAAGCCGTGGTGCGCTGTGCCCGAAAGGGGCCGGTTTGCTGGATTAc                                >  2:4285638/1‑93 (MQ=255)
                                                            ccAGAGAAGCGATTTATCACATTGAAGGTGGCCCGGATCATCCGGTAAGCCGTGGTGCGCTGTGCCCGAAAGGGGCCGGTTTGCTGGATTACg                               >  2:6640467/1‑93 (MQ=255)
                                                             cAGAGAAGCGATTTATCACATTGAAGGTGGCCCGGATCATCCGGTAAGCCGTGGTGCGCTGTGCCCGAAAGGGGCCGGTTTGCTGGATTACGt                              <  2:4097950/93‑1 (MQ=255)
                                                              agagAAGCGATTTATCACATTGAAGGTGGCCCGGATCATCCGGTAAGCCGTGGTGCGCTGTGCCCGAAAGGGGCCGGTTTGCTGGATTACGTc                             >  1:4996472/1‑93 (MQ=255)
                                                              agagAAGCGATTTATCACATTGAAGGTGGCCCGGATCATCCGGTAAGCCGTGGTGCGCTGTGCCCGAAAGGGGCCGGTTTGCTGGATTACGTc                             >  1:4027093/1‑93 (MQ=255)
                                                                 gaAGCGATTTATCACATTGAAGGTGGCCCGGATCATCCGGTAAGCCGTGGTGCGCTGTGCCCGAAAGGGGCCGGTTTGCTGGATTACGTCAAc                          >  2:6482243/1‑93 (MQ=255)
                                                                   aGCGATTTATCACATTGAAGGTGGCCCGGATCATCCGGTAAGCCGTGGTGCGCTGTGCCCGAAAGGGGCCGGTTTGCTGGATTACGTCAACAg                        >  2:1717011/1‑93 (MQ=255)
                                                                   aGCGATTTATCACATTGAAGGTGGCCCGGATCATCCGGTAAGCCGTGGTGCGCTGTGCCCGAAAGGGGCCGGTTTGCTGGATTACGTCAACAg                        <  1:5190629/93‑1 (MQ=255)
                                                                   aGCGATTTATCACATTGAAGGTGGCCCGGATCATCCGGTAAGCCGTGGTGCGCTGTGCCCGAAAGGGGCCGGTTTGCTGGATTACGTCAACAg                        >  2:3377707/1‑93 (MQ=255)
                                                                    gCGATTTATCACATTGAAGGTGGCCCGGATCATCCGGTAAGCCGTGGTGCGCTGTGCCCGAAAGGGGCCGGTTTGCTGGATTACGTCAACAGt                       >  2:872552/1‑93 (MQ=255)
                                                                    gCGATTTATCACATTGAAGGTGGCCCGGATCATCCGGTAAGCCGTGGTGCGCTGTGCCCGAAAGGGGCCGGTTTGCTGGATTACGTCAACAGt                       <  1:1182458/93‑1 (MQ=255)
                                                                      gATTTATCACATTGAAGGTGGCCCGGATCATCCGGTAAGCCGTGGTGCGCTGTGCCCGAAAGGGGCCGGTTTGCTGGATTACGTCAACAGTGa                     >  1:2202191/1‑93 (MQ=255)
                                                                      gATTTATCACATTGAAGGTGGCCCGGATCATCCGGTAAGCCGTGGTGCGCTGTGCCCGAAAGGGGCCGGTTTGCTGGATTACGTCAACAGTGa                     >  2:5913421/1‑93 (MQ=255)
                                                                      gATTTATCACATTGAAGGTGACCCGGATCATCCGGTAAGCCGTGGTGCGCTGTGCCCGAAAGGGGCCGGTTTGCTGGATTACGTCAACAGTGa                     >  2:6766178/1‑93 (MQ=255)
                                                                       aTTTATCACATTGAAGGTGGCCCGGATCATCCGGTAAGCCGTGGTGCGCTGTGCCCGAAAGGGGCCGGTTTGCTGGATTACGTCAACAGTGaa                    <  1:2908147/93‑1 (MQ=255)
                                                                       aTTTATCACATTGAAGGTGACCCGGATCATCCGGTAAGCCGTGGTGCGCTGTGCCCGAAAGGGGCCGGTTTGCTGGATTACGTCAACAGTGaa                    >  2:6289187/1‑93 (MQ=255)
                                                                        tttATCACATTGAAGGTGGCCCGGATCATCCGGTAAGCCGTGGTGCGCTGTGCCCGAAAGGGGCCGGTTTGCTGGATTACGTCAACAGTGaaa                   >  2:3567607/1‑93 (MQ=255)
                                                                        tttATCACATTGAAGGTGGCCCGGATCATCCGGTAAGCCGTGGTGCGCTGTGCCCGAAAGGGGCCGGTTTGCTGGATTACGTCAACAGTGaaa                   >  1:6063806/1‑93 (MQ=255)
                                                                         ttATCACATTGAAGGTGGCCCGGATCATCCGGTAAGCCGTGGTGCGCTGTGCCCGAAAGGGGCCGGTTTGCTGGATTACGTCAACAGTGaaaa                  >  1:111745/1‑93 (MQ=255)
                                                                         ttATCACATTGAAGGTGGCCCGGATCATCCGGTAAGCCGTGGTGCGCTGTGCCCGAAAGGGGCCGGTTTGCTGGATTACGTCAACAGTGaaaa                  <  1:509350/93‑1 (MQ=255)
                                                                         ttATCACATTGAAGGTGGCCCGGATCATCCGGTAAGCCGTGGTGCGCTGTGCCCGAAAGGGGCCGGTTTGCTGGATTACGTCAACAGTGaaaa                  <  1:7137502/93‑1 (MQ=255)
                                                                         ttATCACATTGAAGGTGGCCCGGATCATCCGGTAAGCCGTGGTGCGCTGTGCCCGAAAGGGGCCGGTTTGCTGGATTACGTCAACAGTGaaaa                  <  2:1305037/93‑1 (MQ=255)
                                                                          tATCACATTGAAGGTGGCCCGGATCATCCGGTAAGCCGTGGTGCGCTGTGCCCGAAAGGGGCCGGTTTGCTGGATTACGTCAACAGTGAAAAc                 <  2:3649682/93‑1 (MQ=255)
                                                                          tATCACATTGAAGGTGGCCCGGATCATCCGGTAAGCCGTGGTGCGCTGTGCCCGAAAGGGGCCGGTTTGCTGGATTACGTCAACAGTGAAAAc                 <  2:7317532/93‑1 (MQ=255)
                                                                          tATCACATTGAAGGTGACCCGGATCATCCGGTAAGCCGTGGTGCGCTGTGCCCGAAAGGGGCCGGTTTGCTGGATTACGTCAACAGTGAAAAc                 <  1:7379375/93‑1 (MQ=255)
                                                                           aTCACATTGAAGGTGGCCCGGATCATCCGGTAAGCCGTGGTGCGCTGTGCCCGAAAGGGGCCGGTTTGCTGGATTACGTCAACAGTGAAAAcc                <  2:4924049/93‑1 (MQ=255)
                                                                             cacaTTGAAGGTGGCCCGGATCATCCGGTAAGCCGTGGTGCGCTGTGCCCGAAAGGGGCCGGTTTGCTGGATTACGTCAACAGTGAAAACCGt              <  1:4279928/93‑1 (MQ=255)
                                                                             cacaTTGAAGGTGGCCCGGATCATCCGGTAAGCCGTGGTGCGCTGTGCCCGAAAGGGGCCGGTTTGCTGGATTACGTCAACAGTGAAAACCGt              <  2:16430/93‑1 (MQ=255)
                                                                              acaTTGAAGGTGGCCCGGATCATCCGGTAAGCCGTGGTGCGCTGTGCCCGAAAGGGGCCGGTTTGCTGGATTACGTCAACAGTGAAAACCGTc             <  1:4972490/93‑1 (MQ=255)
                                                                              acaTTGAAGGTGGCCCGGATCATCCGGTAAGCCGTGGTGCGCTGTGCCCGAAAGGGGCCGGTTTGCTGGATTACGTCAACAGTGAAAACCGTc             <  1:6210282/93‑1 (MQ=255)
                                                                              acaTTGAAGGTGACCCGGATCATCCGGTAAGCCGTGGTGCGCTGTGCCCGAAAGGGGCCGGTTTGCTGGATTACGTCAACAGTGAAAACCGTc             >  2:6601853/1‑93 (MQ=255)
                                                                               caTTGAAGGTGACCCGGATCATCCGGTAAGCCGTGGTGCGCTGTGCCCGAAAGGGGCCGGCTTGCTGGATTACGTCAACAGTGAAAACCGTCt            >  2:2009835/1‑93 (MQ=255)
                                                                                aTTGAAGGTGGCCCGGATCATCCGGTAAGCCGTGGTGCGCTGTGCCCGAAAGGGGCCGGTTTGCTGGATTACGTCAACAGTGAAAACCGTCTg           <  2:850926/93‑1 (MQ=255)
                                                                                aTTGAAGGTGACCCGGATCATCCGGTAAGCCGTGGTGCGCTGTGCCCGAAAGGGGCCGGTTTGCTGGATTACGTCAACAGTGAAAACCGTCTg           >  1:1727157/1‑93 (MQ=255)
                                                                                 ttGAAGGTGACCCGGATCATCCGGTAAGCCGTGGTGCGCTGTGCCCGAAAGGGGCCGGTTTGCTGGATTACGTCAACAGTGAAAACCGTCTgc          >  1:1534432/1‑93 (MQ=255)
                                                                                   gAAGGTGGCCCGGATCATCCGGTAAGCCGTGGTGCGCTGTGCCCGAAAGGGGCCGGTTTGCTGGATTACGTCAACAGTGAAAACCGTCTgcgc        >  1:2654772/1‑93 (MQ=255)
                                                                                      ggTGACCCGGATCATCCGGTAAGCCGTGGTGCGCTGTGCCCGAAAGGGGCCGGTTTGCTGGATTACGTCAACAGTGAAAACCGTCTGCGCTAc     >  1:2087395/1‑93 (MQ=255)
                                                                                          |                                                                                           
CATACTGTTCCGTAGGTTGCGGGCTATTGATGTATAGCCTGGGTGATGGCGCAAAAAACGCCAGAGAAGCGATTTATCACATTGAAGGTGACCCGGATCATCCGGTAAGCCGTGGTGCGCTGTGCCCGAAAGGGGCCGGTTTGCTGGATTACGTCAACAGTGAAAACCGTCTGCGCTACCCG  >  NC_000913/1545576‑1545757

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