Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,552,388 T→C 15.2% intergenic (+397/+10) fdnI → / ← yddM formate dehydrogenase‑N, cytochrome B556 (gamma) subunit, nitrate‑inducible/putative DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,552,3880TC15.2% 347.0 / 47.7 144intergenic (+397/+10)fdnI/yddMformate dehydrogenase‑N, cytochrome B556 (gamma) subunit, nitrate‑inducible/putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (61/60);  new base C (9/13);  total (70/74)
Fisher's exact test for biased strand distribution p-value = 4.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.56e-01

Only 100 of 148 total aligned reads displayed.

TTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACAT  >  NC_000913/1552289‑1552487
                                                                                                   |                                                                                                   
ttCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGt                                                                                                     >  1:2284738/1‑100 (MQ=255)
  cAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTCAATGGCGCTTTATAAAGTGc                                                                                                   <  2:2328620/100‑1 (MQ=255)
   aaaTCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGcc                                                                                                  >  2:3721873/1‑100 (MQ=255)
    aaTCCAATAGATTGCAGTGAGCGTGTAGGCCTGATAACCGTAGCGCATCAGGCAATGTTGCGTTGGTCATCAGTTCTAAATGGCGCTTTATCAAGTGGCa                                                                                                 <  2:2573298/100‑1 (MQ=14)
    aaTCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCa                                                                                                 <  2:2835486/100‑1 (MQ=255)
      tCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCACTTTATAAAGCGCCAtt                                                                                               <  2:3595448/100‑1 (MQ=255)
       ccAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCACTTTATAAAGCGCCAttt                                                                                              <  2:140520/100‑1 (MQ=255)
        cAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCAtttt                                                                                             >  1:3646021/1‑100 (MQ=255)
        cAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCACTTTATAAAGCGCCAttta                                                                                             <  1:1007811/100‑2 (MQ=255)
        cAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCACTTTATAAAGCGCCAttta                                                                                             <  1:762652/100‑2 (MQ=255)
        cAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCACTTTATAAAGCGCCAttta                                                                                             <  2:2056524/100‑2 (MQ=255)
        cAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTCGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCAtttt                                                                                             >  2:3309085/1‑100 (MQ=255)
         aaTAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCAttttt                                                                                            >  1:178839/1‑100 (MQ=255)
         aaTAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCACTTTATAAAGCGCCAtttag                                                                                            <  1:2061529/100‑3 (MQ=255)
          tctgATTGCAGTGAACGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCAtttttt                                                                                           <  2:2533475/97‑1 (MQ=255)
          aTAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCAtttttt                                                                                           <  2:2821459/100‑1 (MQ=255)
          aTAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCAtttttt                                                                                           <  2:1743053/100‑1 (MQ=255)
          aTAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATCTCGCTTTATAAAGTGCCAtttttt                                                                                           <  1:767421/100‑1 (MQ=255)
           tAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCAttttttt                                                                                          <  2:339148/100‑1 (MQ=255)
            aGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGACGCCTTTTACAGCGCCAtttgaaac                                                                                         <  1:3650548/100‑6 (MQ=11)
            aGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCACTTTATAAAGCGCCAtttagaac                                                                                         <  1:189539/100‑6 (MQ=17)
             gATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCTTCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTAt                                                                                        <  2:720735/100‑1 (MQ=255)
             gATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTAt                                                                                        <  2:1869838/100‑1 (MQ=255)
             gATTGCAGTGAACGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCAttttttaga                                                                                        >  1:2533491/1‑97 (MQ=255)
              aTTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTAtt                                                                                       <  2:2216676/100‑1 (MQ=255)
               ttGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTg                                                                                      <  1:2678416/100‑1 (MQ=255)
                 gCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCg                                                                                    >  1:3743093/1‑100 (MQ=255)
                 gCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCg                                                                                    <  2:1518831/100‑1 (MQ=255)
                 gCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCg                                                                                    >  2:3654494/1‑100 (MQ=255)
                  cAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGt                                                                                   >  2:916474/1‑100 (MQ=255)
                       aaCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAAc                                                                                <  2:2794220/98‑1 (MQ=255)
                       aaCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCa                                                                              >  2:1977640/1‑100 (MQ=255)
                           tgtAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACg                                                                          >  1:3313988/1‑100 (MQ=255)
                           tgtAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTGTATAACGTGCCATTTTTTTATTGCGTAAcccg                                                                             >  2:2832076/1‑95 (MQ=255)
                            gtAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACgg                                                                         >  2:2736739/1‑100 (MQ=255)
                              aGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCg                                                                       <  1:657431/100‑1 (MQ=255)
                               ggCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGt                                                                      >  1:3221327/1‑100 (MQ=255)
                               ggCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGt                                                                      >  2:1631129/1‑100 (MQ=255)
                                gCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTa                                                                     <  2:1892646/100‑1 (MQ=255)
                                  cTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAAt                                                                   >  2:3549187/1‑100 (MQ=255)
                                  cTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCAATTTTTTTATTGCGTAACCAGACGGCGTAAt                                                                   >  2:1809202/1‑100 (MQ=255)
                                    gATAAGCGTAGCACATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATcg                                                                 >  2:2906750/1‑100 (MQ=255)
                                      tAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATcgcg                                                               >  1:420736/1‑100 (MQ=255)
                                      tAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATcg                                                                 >  2:741881/1‑98 (MQ=255)
                                        aGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGac                                                             <  2:975359/100‑1 (MQ=255)
                                        aGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGac                                                             <  1:2764129/100‑1 (MQ=255)
                                        aGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGac                                                             <  1:181260/100‑1 (MQ=255)
                                        aGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGac                                                             >  2:311775/1‑100 (MQ=255)
                                         gCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTATAAGGCGTCATTTTTATATTGCGTACACAGTCGGCTTAGTCGCTac                                                             >  2:2838489/1‑99 (MQ=14)
                                           gTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGAcaca                                                          >  1:3277587/1‑100 (MQ=255)
                                                gcATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCAc                                                     >  1:110372/1‑100 (MQ=255)
                                                gcATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCAc                                                     >  1:1443206/1‑100 (MQ=255)
                                                gcATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCAc                                                     >  1:3720934/1‑100 (MQ=255)
                                                 cATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACt                                                    >  2:1963978/1‑100 (MQ=255)
                                                  aTCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTg                                                   <  1:1534390/100‑1 (MQ=255)
                                                     aGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACAACCACTGtt                                                 >  2:1757815/1‑99 (MQ=255)
                                                           tGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTtctg                                          <  1:3347925/100‑1 (MQ=255)
                                                           tGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTtct                                           >  2:1059787/1‑99 (MQ=255)
                                                            ggtGCGTTGGGCAACAGTTCTATATGGCGCTTCATAAAGGGCGATTTTTTTATTGCGTAGCAGGACGGAGTAATCGCGACTTATCCGCAGTTTTTtctg                                          <  2:1454692/97‑1 (MQ=255)
                                                             ttGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTtctgcc                                        <  1:1231986/100‑1 (MQ=255)
                                                              tGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTtctgcct                                       >  1:2693369/1‑100 (MQ=255)
                                                              tGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTtctgcc                                        <  2:1038573/99‑1 (MQ=255)
                                                               gCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGCCGGCGTAATCGCGACACATTCACTGTTTTTtctgcct                                       >  2:371356/1‑99 (MQ=255)
                                                                 gTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTtctgcctctg                                    >  1:82066/1‑100 (MQ=255)
                                                                 gTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTtctgcctctg                                    >  2:470206/1‑100 (MQ=255)
                                                                       catcaGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGAc                              >  2:1471440/1‑100 (MQ=255)
                                                                        gcaataaaaaaaTGGCACTTTATAAAGCGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACt                             >  1:613827/10‑100 (MQ=37)
                                                                        atcaGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCTACACATCCACTGTTTTTTCTGCCTCTGCCAGACt                             >  2:1173456/1‑100 (MQ=255)
                                                                          caGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTcc                           >  2:545078/1‑100 (MQ=255)
                                                                           aGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCa                          <  1:1977651/100‑1 (MQ=255)
                                                                            gTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAg                         <  2:2693350/100‑1 (MQ=255)
                                                                            gTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAg                         >  1:3242741/1‑100 (MQ=255)
                                                                             aaaaaaaTGGCACTTTATAAAGCGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAgg                        >  2:2863766/5‑100 (MQ=38)
                                                                             aaaaaaaTGGCACTTTATAAAGCGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAgg                        >  1:2425884/5‑100 (MQ=38)
                                                                              aaaaaaTGGCACTTTATAAATCGCCATTTTTTTATTTCGTAACCCGGCGGCGTAAGCGCGGCACATTTCCTGTTTTTTCTGCCTCTTCCAGtctc                            >  2:3433859/4‑95 (MQ=16)
                                                                              aaaaaaTGGCACTTCATAAAGCGCCATTTTTTTAATGCGTACCCCGACTGCCTAATCGCGAAACTTCCGGTGTATTCTCTTCCTCTGCCAGACTcc                           >  2:3591755/4‑96 (MQ=16)
                                                                                 aaaTGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGtt                    <  1:87879/100‑1 (MQ=255)
                                                                                  aaTGGCACTTTATAAAGCGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGttt                   >  1:3017645/1‑100 (MQ=255)
                                                                                  aaTGGCACTTTATAAAGCGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGttt                   >  1:2242398/1‑100 (MQ=255)
                                                                                    tGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTg                 >  1:697667/1‑100 (MQ=255)
                                                                                    tGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTg                 <  1:2414891/100‑1 (MQ=255)
                                                                                     ggCGCTTTATAAAGTGCCATTTTTTTATTGCGTAAGCTGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCACGCGTTTTGc                >  2:2064054/1‑100 (MQ=255)
                                                                                     ggCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGc                >  2:2818934/1‑100 (MQ=255)
                                                                                      gcgcTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCa               <  2:2860299/100‑1 (MQ=255)
                                                                                      gcgcTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCa               <  2:1333951/100‑1 (MQ=255)
                                                                                      gcgcTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCa               <  2:420734/100‑1 (MQ=255)
                                                                                      gcgcTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCa               >  1:461869/1‑100 (MQ=255)
                                                                                      gcgcTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCa               >  1:2503695/1‑100 (MQ=255)
                                                                                      gcgcTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCa               <  1:2838509/100‑1 (MQ=255)
                                                                                      gcgcTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCa               >  2:3081460/1‑100 (MQ=255)
                                                                                         cTTTATAAAGTGCCATTTTTTTATTGCCTAACCAGACGGCGTAATCGCCACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGAt            <  2:2797483/100‑1 (MQ=255)
                                                                                            tataAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCa         >  1:2300912/1‑100 (MQ=255)
                                                                                             ataAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAg        <  2:2173850/100‑1 (MQ=255)
                                                                                              taAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTTAGc       >  2:569352/1‑100 (MQ=255)
                                                                                              taAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATAGCGACACATCCACTGTTTTTTATGCCTCTGCCAGACTCAAGGCGATTTGCAGAATCAg        >  1:317089/1‑99 (MQ=255)
                                                                                               gatctGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGcc      <  1:591789/96‑1 (MQ=255)
                                                                                               aaaGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTTTTCTTTTTCCAGACTCCAGGCGTTTTGCAGATTCAGc       >  2:2587992/1‑99 (MQ=255)
                                                                                               aaaGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGcc      <  1:3591060/100‑1 (MQ=255)
                                                                                               aaaGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGcc      <  1:3591001/100‑1 (MQ=255)
                                                                                                   tGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACAt  >  1:931987/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
TTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACAT  >  NC_000913/1552289‑1552487

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