Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,564,601 A→G 40.1% F378S (TTT→TCT)  dosP ← oxygen sensor, c‑di‑GMP phosphodiesterase, heme‑regulated; cold‑ and stationary phase‑induced bioflim regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,564,6010AG40.1% 45.8 / 116.1 106F378S (TTT→TCT) dosPoxygen sensor, c‑di‑GMP phosphodiesterase, heme‑regulated; cold‑ and stationary phase‑induced bioflim regulator
Reads supporting (aligned to +/- strand):  ref base A (35/27);  new base G (23/19);  total (60/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.29e-01

Only 100 of 109 total aligned reads displayed.

CCGATGAGATACACCACGGGAGAGACGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGATAAAGGCGC  >  NC_000913/1564505‑1564700
                                                                                                |                                                                                                   
ccGATGAGATACACCACGGGAGAGACGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGAtt                                                                                                  <  1:1284564/100‑1 (MQ=255)
ccGATGAGATACACCACGGGAGAGACGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAAt                                                                                                   >  2:2731626/1‑99 (MQ=255)
ccGATGAGATACACCACGGGAGAGACGGCTCTGTCGAGAAGGTGATCGAGGTCGTTGTGCTAAATATTACAATTTAGCAGACCGGTCATCGGATCAAAtt                                                                                                  >  1:4050685/1‑100 (MQ=255)
 cgttgatgtACACCAAGGGATAGACGGCTTTGTCGATCATGTCATCGAGGTAATTGGGCAGGTTATTGCGATTTGGCAGACCAGTCATCGGCTCAAATTg                                                                                                 <  1:662449/92‑1 (MQ=255)
  gATGAGATACACCACGGGAGAGACGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTgg                                                                                                <  2:2820305/100‑1 (MQ=255)
  gATGAGATACACCACGGGAGAGACGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTgg                                                                                                >  1:2981722/1‑100 (MQ=255)
  gATGAGATACACCACGGGAGAGACGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTg                                                                                                 <  2:3448182/99‑1 (MQ=255)
    tGAGATACACCACGGGAGAGACGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGAt                                                                                              <  1:3231846/100‑1 (MQ=255)
     gagaTACACCACGGGAGAGACGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATg                                                                                             <  1:778639/100‑1 (MQ=255)
       taTACACCACGTGAGAGTCGGCTTTGTCGACCCGGTCATCGCGGTAATTGTGCAGGTTCTTGCGATTTGGCACACCGGTCATCGGATCAAATCgg                                                                                                >  1:1906131/2‑95 (MQ=255)
        aTACACCACGGGAGAGACGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAgt                                                                                          <  2:823504/100‑1 (MQ=255)
         tACACCACGGGAGAGACGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAgtt                                                                                         >  2:3544192/1‑100 (MQ=255)
          acacCACGGGAGAGACGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAgttg                                                                                        <  1:1352160/100‑1 (MQ=255)
           caccacGGGAGAGACGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAgttgt                                                                                       >  1:3571327/1‑100 (MQ=255)
           caccacGGGAGAGACGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGGTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGAAGAggtgt                                                                                       >  1:256811/1‑100 (MQ=255)
             ccacGGGAGAGACGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTc                                                                                     <  2:753177/100‑1 (MQ=255)
             ccacGGGAGAGACGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTc                                                                                     <  2:3979428/100‑1 (MQ=255)
              cacGGGAGAGACGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCa                                                                                    <  2:3729586/100‑1 (MQ=255)
              cacGGGAGAGACGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCa                                                                                    <  1:883581/100‑1 (MQ=255)
               acGGGAGAGACGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCaa                                                                                   <  1:1677909/100‑1 (MQ=255)
                  ggAGAGACGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAatat                                                                                >  2:1134836/1‑100 (MQ=255)
                    agagaCGGCTTTGTCGACCAGGACATCGAGGTAATTGTGCAGGTAATTGCGATGTGGCAGACCGGCAATCGGAACAAATCGGACGCGTTGCTc                                                                                     >  1:3453553/1‑93 (MQ=255)
                      agaCGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGAATGAGTTGTTCAATATGCTg                                                                            <  1:1083558/100‑1 (MQ=255)
                         cGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACg                                                                         >  2:1033685/1‑100 (MQ=255)
                           gCTTTGTCGACGAGGTCATCGAGGTAATTGTACAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATTACTTTTTCAGTATGCTGACGGc                                                                       <  1:1470743/100‑1 (MQ=255)
                           gCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGc                                                                       >  1:1562603/1‑100 (MQ=255)
                           gCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGc                                                                       >  1:921634/1‑100 (MQ=255)
                           gCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGc                                                                       >  2:3794398/1‑100 (MQ=255)
                             tttGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCtt                                                                     >  1:2126443/1‑100 (MQ=255)
                                gTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGCttttt                                                                  >  1:1970876/1‑100 (MQ=255)
                                   gACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTcc                                                               >  2:2068424/1‑100 (MQ=255)
                                   gACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTcc                                                               >  1:235848/1‑100 (MQ=255)
                                   gACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTcc                                                               <  2:2151559/100‑1 (MQ=255)
                                     ccAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTg                                                             <  2:2971529/100‑1 (MQ=255)
                                      cAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGt                                                            <  1:4239491/100‑1 (MQ=255)
                                         gTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTcc                                                         <  1:3323540/100‑1 (MQ=255)
                                           cATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTTTTGtt                                                           >  1:1438625/1‑96 (MQ=255)
                                           cATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTccag                                                       >  1:2382026/1‑100 (MQ=255)
                                           cATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTccag                                                       >  2:1737869/1‑100 (MQ=255)
                                           cATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTAGATGAGTTGTTAACTATGCTGACGGCTTTTTTCCTg                                                             >  1:3521465/1‑94 (MQ=255)
                                             tCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTccagcg                                                     >  1:2468058/1‑100 (MQ=255)
                                             tCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTccagcg                                                     >  2:823175/1‑100 (MQ=255)
                                              cGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCACATTCGATGAGTTGTTCAATATGCTCACGCACTTTTTACTGATacagaga                                                    >  2:1846551/1‑97 (MQ=255)
                                              cGAGGTAATTGTGCAGGTTATAGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTcga                                                        >  1:3544220/1‑94 (MQ=255)
                                               gAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTccagcgcc                                                   >  2:2918087/1‑100 (MQ=255)
                                                  gTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTccagcgccagc                                                >  2:2725675/1‑100 (MQ=255)
                                                       tgtgCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGc                                           <  1:639205/100‑1 (MQ=255)
                                                       tgtgCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGc                                           >  2:2121675/1‑100 (MQ=255)
                                                        gtgCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGcc                                          >  1:301404/1‑100 (MQ=255)
                                                         tgCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCa                                         <  1:2068420/100‑1 (MQ=255)
                                                          gCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCat                                        >  2:1904732/1‑100 (MQ=255)
                                                           cAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCata                                       <  1:234349/100‑1 (MQ=255)
                                                             ggTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATg                                     <  2:1952945/100‑1 (MQ=255)
                                                             ggTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATg                                     >  2:1100308/1‑100 (MQ=255)
                                                              gTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATgc                                    >  2:2275041/1‑100 (MQ=255)
                                                                tATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATgctg                                  <  1:2061969/100‑1 (MQ=255)
                                                                 aTTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATgctgg                                 >  2:486109/1‑100 (MQ=255)
                                                                  ttGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATgctggc                                <  1:2178600/100‑1 (MQ=255)
                                                                  ttGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATgctggc                                >  2:4137016/1‑100 (MQ=255)
                                                                  ttGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGAAGGCTTTTTTCCTGTTCTAGCGCCAGCGCGGCCATATgctggc                                >  2:2107308/1‑100 (MQ=255)
                                                                   tGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATgctggct                               >  2:2274914/1‑100 (MQ=255)
                                                                   tGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATgctggct                               <  1:32899/100‑1 (MQ=255)
                                                                    gCGATTTGGCAGACCGGTCATCGGGTCAAATGGGATGAGGTGTTCAATATGCTGACTGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCat                                        >  1:1183698/1‑90 (MQ=255)
                                                                    gCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATgctggctg                              <  1:2465451/100‑1 (MQ=255)
                                                                     cGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGa                             <  2:3548448/100‑1 (MQ=255)
                                                                     cGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGa                             <  2:294868/100‑1 (MQ=255)
                                                                     cGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGa                             <  2:46186/100‑1 (MQ=255)
                                                                      gATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGa                             >  1:738351/1‑99 (MQ=255)
                                                                       aTTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCCGCGCGGCCATATGCTGGCTGata                           >  1:529416/1‑100 (MQ=255)
                                                                       aTTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGata                           >  2:570761/1‑100 (MQ=255)
                                                                         ttGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTTCTTTTTCAGCGCCAGCGCGGCCATATGCTGGCTGata                           >  1:1326642/1‑98 (MQ=255)
                                                                           ggCAGACCTGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTg                       >  1:3221422/1‑100 (MQ=255)
                                                                            gCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGc                      >  2:2596765/1‑100 (MQ=255)
                                                                             cAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGcc                     <  1:2918075/100‑1 (MQ=255)
                                                                             cAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGcc                     <  2:2346210/100‑1 (MQ=255)
                                                                             cAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGcc                     >  1:2615042/1‑100 (MQ=255)
                                                                                aCCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCCcg                  >  1:3758914/1‑100 (MQ=255)
                                                                                 ccGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCAcgc                 <  1:2691154/100‑1 (MQ=255)
                                                                                 ccGGCCGTCGGATTACACCGGATGAGTTGGTCAATATGCTCACGGCTTTTATCCTGTTCCATCGCCAGGCCGGCCATTTGCTGGCTGATATCTCTCAcgc                 >  1:3422323/1‑100 (MQ=255)
                                                                                   ggTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGt               >  2:3875812/1‑100 (MQ=255)
                                                                                    gTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGtt              <  2:2815248/100‑1 (MQ=255)
                                                                                    gTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGtt              >  1:2429085/1‑100 (MQ=255)
                                                                                     tCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTc             <  2:57083/100‑1 (MQ=255)
                                                                                     tCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTc             <  2:3546652/100‑1 (MQ=255)
                                                                                      cATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCg            <  2:1700608/100‑1 (MQ=255)
                                                                                       aTCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGa           >  2:3182498/1‑100 (MQ=255)
                                                                                       aTCGGATCAAATTGGATAAGTTGTTCGATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGTCACGCGTTCGa           <  1:145239/100‑1 (MQ=255)
                                                                                        tCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGAt          >  1:998045/1‑100 (MQ=255)
                                                                                         cGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGATa         <  2:1604460/100‑1 (MQ=255)
                                                                                          ggATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGATaa        >  2:471123/1‑100 (MQ=255)
                                                                                           gATCAGATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGATaaa       >  1:3195223/1‑100 (MQ=255)
                                                                                           gATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGATaaa       >  2:3863886/1‑100 (MQ=255)
                                                                                           gATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGATaaa       <  1:4153402/100‑1 (MQ=255)
                                                                                            aTCAGATTGGATGAGTTGTTCAATATGCTGACGGCTTGTTTCCCGTTCCAGCGCCACCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGATAAAg      <  2:416218/100‑1 (MQ=255)
                                                                                             tCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGATAAAgg     <  2:744725/100‑1 (MQ=255)
                                                                                             tCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGATAAAgg     <  2:1859584/100‑1 (MQ=255)
                                                                                              cAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGATAAAGgc    >  2:877486/1‑100 (MQ=255)
                                                                                               aaaTTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGATAAAGgcg   >  1:1439780/1‑100 (MQ=255)
                                                                                                aaTTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTTTTTTTTTTTTTCCAGCGCGGCCATATGCTGGCTTATATCTGCCACGCGTTCGATAAAGgcg   >  1:2918159/1‑99 (MQ=255)
                                                                                                aaTTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGATAAAGgcgc  >  1:302400/1‑100 (MQ=255)
                                                                                                |                                                                                                   
CCGATGAGATACACCACGGGAGAGACGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGATAAAGGCGC  >  NC_000913/1564505‑1564700

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