Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,303,6410AT27.2% 170.1 / 71.6 93intergenic (‑144/‑264)napF/ecoferredoxin‑type protein, role in electron transfer to periplasmic nitrate reductase NapA/ecotin, a serine protease inhibitor
Reads supporting (aligned to +/- strand):  ref base A (33/34);  new base T (15/10);  total (49/44)
Fisher's exact test for biased strand distribution p-value = 4.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAATAAACCGCATAATTAATAAGCCATTTTTATAGCCGCTAAGATATTAAAGGATGTGTCAAAGATGCATACCCCGATCGGGGTAAAACCTGTAGCAGGATCAAAAGAGTGGCGGGAAGCGTGGCAAAAACGGGCTTTTGCTCACATTTCAAATGGTTATAAATATATTTATATAGCGATTGATTCACCAGAGA  >  NC_000913/2303542‑2303735
                                                                                                   |                                                                                              
aaaTAAACCGCATAATTAATAAGCCATTTTTATAGCCGCTAAGATATTAAAGGATGTGTCAAAGATGCATACCCCGATCGGGGTAAAACCTGTAGCAGGa                                                                                                <  2:582510/100‑1 (MQ=255)
aaaTAAACCGCATAATTAATAAGCCATTTTTATAGCCGCTAAGATATTAAAGGATGTGTCAAAGATGCATACCCCGATCGGGGTAAAACCTGTAGCAGGa                                                                                                <  1:1239345/100‑1 (MQ=255)
  aTAAACCGCATAATTAATAAGCCATTTTTATAGCCGCTAAGATATTAAAGGATGTGTCAAAGATGCATACCCCGATCGGGGTAAAACCTGTAGCAGGTTc                                                                                              <  1:523191/100‑1 (MQ=255)
  aTAAACCGCATAATTAATAAGCCATTTTTATAGCCGCTAAGATATTAAAGGATGTGTCAAAGATGCATACCCCGATCGGGGTAAAACCTGTAGCAGGATc                                                                                              <  2:1099904/100‑1 (MQ=255)
       ccGCATAATTAATAAGCCATTTTTATAGCCGCTAAGATATTAAAGGATGTGTCAAAGATGCATACCCCGATCGGGGTAAAACCTGTAGCAGGTTCAAAag                                                                                         <  1:537582/100‑1 (MQ=255)
       ccGCATAATTAATAAGCCATTTTTATAGCCGCTAAGATATTAAAGGATGTGTCAAAGATGCATACCCCGATCGGGGTAAAACCTGTAGCAGGATCAAAag                                                                                         <  1:1621042/100‑1 (MQ=255)
        cGCATAATTAATAAGCCATTTTTATAGCCGCTAAGATATTAAAGGATGTGTCAAAGATGCATACCCCGATCGGGGTAAAACCTGTAGCAGGATCAAAaga                                                                                        >  1:1440829/1‑100 (MQ=255)
         gCATAATTAATAAGCCATTTTTATAGCCGCTAAGATATTAAAGGATGTGTCAAAGATGCATACCCCGATCGGGGTAAAACCTGTAGCAGGATCAAAagag                                                                                       >  1:1037104/1‑100 (MQ=255)
             aattaatAAGCCATTTTTATAGCCGCTAAGATATTAAAGGATGTGTCAAAGATGCATACCCCGATCGGGGTAAAACCTGTAGCAGGTTCAAAAGAGTGGc                                                                                   >  1:269403/1‑100 (MQ=255)
               ttaatAAGCCATTTTTATAGCCGCTAAGATATTAAAGGATGTGTCAAAGATGCATACCCCGATCGGGGTAAAACCTGTAGCAGGATCAAAAGAGTGGCgg                                                                                 <  2:871763/100‑1 (MQ=255)
                 ccgatctcATTTTTATAGCCGCTAAGATATTAAAGGATGTGTCAAAGATGCATACCCCGATCGGGGTAAAACCTGTAGCAGGATCAAAAGAGTGGCGGGa                                                                               <  1:1233789/93‑1 (MQ=255)
                    aaGCCATTTTTATCTCCGCTAAGATATTAACGGGTGTGTCAATGATGCATACCCCGTTCTGGGTAAACCCTGTAGCAGGCTTAAAAGACTGGCGGGAAGc                                                                            >  2:483899/1‑100 (MQ=255)
                    aaGCCATTTTTATAGCCGCTAAGATATTAAAGGATGTGTCAAAGATGCATACCCCGATCGGGGTAAAACCTGTAGCAGGATCAAAAGAGTGGCGGGAAg                                                                             >  2:2318564/1‑99 (MQ=255)
                    aaGCCATTTTTATAGCCGCTAAGATATTAAAGGATGTGTCAAAGATGCATACCCCGATCGGGGTAAAACCTGTAGCAGGATCAAAAGAGTGGCGGGAAGc                                                                            <  1:501241/100‑1 (MQ=255)
                    aaGCCATTTTTATAGCCGCTAAGATATGAAAGGATGTGTCAAAGATGCATACCCCGATCGGGGTAAAACCTGTAGCAGGATCAAAAGAGTGGCGGGAAGc                                                                            >  1:1806853/1‑100 (MQ=255)
                     aGCCATTTTTATAGCCGCTAAGATATTAAAGGATGTGTCAAAGATGCATACCCCGATCGGGGTAAAACCTGTAGCAGGATCAAAAGAGTGGCGGGAAGCg                                                                           >  1:2047813/1‑100 (MQ=255)
                       ccATTTTTATAGCCGCTAAGATATTAAAGGATGTGTCAAAGATGCATACCCCGATCGGGGTAAAACCTGTAGCAGGATCAAAAGAGTGGCGGGAAGCGTg                                                                         <  1:1878211/100‑1 (MQ=255)
                       ccATTTTTATAGCCGCTAAGATATTAAAGGATGTGTCAAAGATGCATACCCCGATCGGGGTAAAACCTGTAGCAGGATCAAAAGAGCGGCGGGAAGCGTg                                                                         <  2:1440825/100‑1 (MQ=255)
                        cATTTTTATAGCCGCTAAGATATTAAAGGATGTGTCAAAGATGCATACCCCGATCGGGGTAAAACCTGTAGCAGGATCAAAAGAGTGGCGGGAAgatcgg                                                                        >  2:1233786/1‑95 (MQ=255)
                        cATTTTTATAGCCGCTAAGATATTAAAGGATGTGTCAAAGATGCATACCCCGATCGGGGTAAAACCTGTAGCAGGATCAAAAGAGTGGCGGGAAGCGTgg                                                                        >  1:1917403/1‑100 (MQ=255)
                         aTTTTTATAGCCGCTAAGATATTAAAGGATGTGTCAAAGATGCATACCCCGATCGGGGTAAAACCTGTAGCAGGATCAAAAGAGTGGCGGGAAGCGTGGc                                                                       >  2:1867974/1‑100 (MQ=255)
                          tttttATAGCCGCTAAGATATTAAAGGATGTGTCAAAGATGCATACCCCGATCGGGGTAAAACCTGTAGCAGGATCAAAAGAGTGGCGGGAAGCGTGGCa                                                                      <  1:2037245/100‑1 (MQ=255)
                             ttATAGCCGCTAAGATATTAAAGGATGTGTCAAAGATGCATACCCCGATCGGGGTAAAACCTGTAGCAGGATCAAAAGAGTGGCGGGAAGCGTGGCaaaa                                                                   <  2:738541/100‑1 (MQ=255)
                               ataGCCGCTAAGATATTAAAGGATGTGTCAAAGATGCATACCCCGATCGGGGTAAAACCTGTAGCAGGTTCAAAAGAGTGGCGGGAAGCGTGGCAAAAAc                                                                 >  1:1265743/1‑100 (MQ=255)
                                   ccGCTAAGATATTAAAGGATGTGTCAAAGATGCATACCCCGATCGGGGTAAAACCTGTAGCAGGTTCAAAAGAGTGGCGGGAAGCGTGGCAAAAACGGGc                                                             <  1:461078/100‑1 (MQ=255)
                                    cGCTAAGATATTAAAGGATGTGTCAAAGATGCATACCCCGATCGGGGTAAAACCTGTAGCAGGTTCAAAAGAGTGGCGGGAAGCGTGGCAAAAACGGGCt                                                            <  2:1084707/100‑1 (MQ=255)
                                     gCTAAGATATTAAAGGATGTGTCAAAGATGCATACCCCGATCGGGGTAAAACCTGTAGCAGGATCAAAAGAGTGGCGGGAAGCGTGGCAAAAACGGGCtt                                                           >  1:1129043/1‑100 (MQ=255)
                                      cTAAGATATTAAAGGATGTGTCAAAGATGCATACCCCGATCGGGGTAAAACCTGTAGCAGGTTCAAAAGAGTGGCGGGAAGCGTGGCAAAAACGGGCttt                                                          >  2:1347493/1‑100 (MQ=255)
                                      cTAAGATATTAAAGGATGTGTCAAAGATGCATACCCCGATCGGGGTAAAACCTGTAGCAGGATCAAAAGAGTGGCGGGAAGCGTGGCAAAAACGGGCttt                                                          >  1:1167600/1‑100 (MQ=255)
                                       tAAGATATTAAAGGATGTGTCAAAGATGCATACCCCGATCGGGGTAAAACCTGTAGCAGGATCAAAAGAGTGGCGGGAAGCGTGGCAAAAACGGGCtttt                                                         <  1:137928/100‑1 (MQ=255)
                                        aaGATATTAAAGGATGTGTCAAAGATGCATACCCCGATCGGGGTAAAACCTGTAGCAGGTTCAAAAGAGTGGCGGGAAGCGTGGCAAAAACGGGCTTTTg                                                        <  2:762599/100‑1 (MQ=255)
                                        aaGATATTAAAGGATGTGTCAAAGATGCATACCCCGATCGGGGTAAAACCTGTAGCAGGATCAAAAGAGTGGCGGGAAGCGTGGCAAAAACGGGCTTTTg                                                        <  1:1867980/100‑1 (MQ=255)
                                          gATATTAAAGGATGTGTCAAAGATGCATACCCCGATCGGGGTAAAACCTGTAGCAGGATCAAAAGAGTGGCGGGAAGCGTGGCAAAAACGGGCTTTTGCt                                                      >  2:822271/1‑100 (MQ=255)
                                           atatTAAAGGATGTGTCAAAGATGCATACCCCGATCGGGGTAAAACCTGTAGCAGGTTCAAAAGCGTGGCGGGAAACGTGGCAAAAACGGGCTTTTGCTc                                                     >  2:570697/1‑100 (MQ=255)
                                            tatTAAAGGTTGTGACAAAGATGCATACCCCGATCGGGGTAAAACCTATAGCAGGATCAAAAGAGTGGTGAAAAAAAAAAAAAAAACGGGCTTTTGCTc                                                     <  2:1531760/99‑1 (MQ=255)
                                              ttAAAGGATGTGTCAAAGATGCATACCCCGATCGGGGTAAAACCTGTAGCAGGATCAAAAGAGTGGCGGGAAGCGTGGCAAAAACGGGCTTTTGCTcaca                                                  >  1:1963053/1‑100 (MQ=255)
                                               tAAAGGATGTGTCAAAGATGCATACCCCGATCGGGGTAAAACCTGTAGCAGGTTCAAAAGAGTGGCGGGAAGCGTGGCAAAAACGGGCTTTTGCTCACAt                                                 >  1:1293112/1‑100 (MQ=255)
                                                aaaGGATGTGTCAAAGATGCATACCCCGATCGGGGTAAAACCTGTAGCAGGATCAAAAGAGTGGCGGGAAGCGTGGCAAAGAAGGGCTTTTGCTcac                                                   >  1:1964357/1‑97 (MQ=255)
                                                 aaGGATGTGTCAAAGATGCATACCCCGATCGGGGTAAAACCTGTAGCAGGTTCAAAAGAGTGGCGGGAAGCGTGGCAAAAACGGGCTTTTGCTCACAttt                                               >  2:972273/1‑100 (MQ=255)
                                                 aaGGATGTGTCAAAGATGCATACCCCGATCGGGGTAAAACCTGTAGCACGATCAAAAGAGTGGCGGGAAGCGTGGCAAAAACGGGCTTTTGCTCACAtt                                                <  1:797747/99‑1 (MQ=255)
                                                  aGGATGTGTCAAAGATGCATACCCCGATCGGGGTAAAACCTGTAGCAGGATCAAAAGAGTGGCGGGAAGCGTGGCAAAAACGGGCTTTTGCTCACATTTc                                              <  1:629551/100‑1 (MQ=255)
                                                   ggATGTGTCAAAGATGCCTACCCCGATCGGGGTAAAACCTGTAGCAGGTTCAAAAGAGTGGCGGGAAGCGTGGCAAACACGGGCGTTTGCTCACATTTCa                                             >  2:708820/1‑100 (MQ=255)
                                                   ggATGTGTCAAAGATGCATACCCCGATCGGGGTAAAACCTGTAGCAGGTTCAAAAGAGTGGCGGGAAGCGTGGCAAAAACGGGCTTTTGCTCACATTTCa                                             <  2:1102084/100‑1 (MQ=255)
                                                   ggATGTGTCAAAGATGCATACCCCGATCGGGGTAAAACCTGTAGCAGGTTCAAAAGAGTGGCGGGAAGCGTGGCAAAAACGGGCTTTTGCTCACATTTCa                                             <  1:943142/100‑1 (MQ=255)
                                                   ggATGTGTCAAAGATGCATACCCCGATCGGGGTAAAACCTGTAGCAGGATCAAAAGAGTGGCGGGAAGCGTGGCAAAAACGGGCTTTTGCTCACATTTCa                                             >  1:1512058/1‑100 (MQ=255)
                                                      tgtgTCAAAGATGCATACCCCGATCGGGGTAAAACCTGTAGCAGGATCAAAAGAGTGGCGGGAAGCGTGGCAAAAACGGGCTTTTGCTCACATTTCAAAt                                          <  2:2185395/100‑1 (MQ=255)
                                                      tgtgTCAAAGATGCATACCCCGATCGGGGTAAAACCTGTAGCAGGATCAAAAGAGTGGCGGGAAGCGTGGCAAAAACGGGCTTTTGCTCACATTTCAAAt                                          <  1:1717419/100‑1 (MQ=255)
                                                       gtgtCAAAGATGCATACCCCGATCGGGGTAAAACCTGTAGCAGGATCAAAAGAGTGGCGGGAAGCGTGGCAAAAACGGGCTTTTGCTCACATTTCAAATg                                         >  1:287573/1‑100 (MQ=255)
                                                       gtgggaaaGATGCATATGGTGATCGTGGTAATAGCTGTAGTTGGATGAAAAGAGTGGTGGTAAGTGTTGCAAAAACGGGCTTTTGCTCATATTTCAAATg                                         <  2:1964503/95‑1 (MQ=255)
                                                         gtCAAAGATGCATACCCCGATCGGGGTAAAACCTGTAGCAGGTTCAAAAGAGTGGCGGGAAGCGTGGCAAAAACGGGCTTTTGCTCACATTTCAAATGGt                                       >  1:989879/1‑100 (MQ=255)
                                                         gtCAAAGATGCATACCCCGATCGGGGTAAAACCTGTAGCAGGATCAAAAGAGTGGCGGGAAGCGTGGCAAAAACGGGCTTTTGCTCACATTTCAAATGGt                                       >  1:133644/1‑100 (MQ=255)
                                                         gtCAAAGATGCATACCCCGATCGGGGTAAAACCTGTAGCAGGATCAAAAGAGTGGCGGGAAGCGTGGCAAAAACGGGCTTTTGCTCACATTTCAAATGGt                                       >  1:2282252/1‑100 (MQ=255)
                                                         gtCAAAGATGCATACCCCGATCGGGGTAAAACCTGTAGCAGGATCAAAAGAGTGGCGGGAAGCGTGGCAAAAACGGGCTTTTGCTCACATTTCAAATGGt                                       <  2:1167597/100‑1 (MQ=255)
                                                         gtCAAAGATGCATACCCCGATCGGGGTAAAACCTGTAGCAGGATCAAAAGAGTGGCGGGAAGCGTGGCAAAAACGGGCTTTTGCTCACATTTCAAATGGt                                       >  1:162732/1‑100 (MQ=255)
                                                         gtCAAAGATGCATACCCCGATCAGGGTAAAACCTGTAGCAGGATCAAAAGAGTGGCGGGAAGCGTGGCAAAAACGGGCTTTTGCTCACATTTCAAATGGt                                       <  1:326774/100‑1 (MQ=255)
                                                         ctCAAAGATGCATACCCCGATCGGGGTAAAACCTGTAGCAGGTTCAAAAGAGTGGCGGGAAGCGTGGTAAAAACGGGCTTTTGCTCACATTTCAAATGGt                                       <  1:2214318/99‑1 (MQ=255)
                                                          tCAAAGATGCATACCCCGATCGGGGTAAAACCTGTAGCAGGTTCAAAAGAGTGGCGGGAAGCGTGGCAAAAACGGGCTTTTGCTCACATTTCAAATGGta                                      >  2:2214310/1‑99 (MQ=255)
                                                          tCAAAGATGCATACCCCGATCGGGGTAAAACCTGTAGCAGGATCAAAAGAGTGGCGGGAAGCGTGGCAAAAACGGGCTTTTGCTCACATTTCAAATGGtt                                      <  2:631926/100‑1 (MQ=255)
                                                          tCAAAGATGCATACCCCGATCGGGGTAAAACCTGTAGCAGGATCAAAAGAGTGGCGGGAAGCGTGGCAAAAACGGGCTTTTGCTCACATTTCAAATGGtt                                      <  1:638701/100‑1 (MQ=255)
                                                          tCAAAGATGCATACCCCGATCGGGGTAAAACCTGTAGCAGGATCAAAAGAGTGGCGGGAAGCGTGGCAAAAACGGGCTTTTGCTCACATTTCAAATGGtt                                      >  1:2167815/1‑100 (MQ=255)
                                                           cAAAGATGCATACCCCGATCGGGGTAAAACCTGTAGCAGGTTCAAAAGAGTGGCGGGCAGCGTGGCAAAAACGGGGTTTTGCTCACATTTCAAATGGTta                                     >  1:2282326/1‑100 (MQ=255)
                                                           cAAAGATGCATACCCCGATCGGGGTAAAACCTGTAGCAGGTTCAAAAGAGTGGCGGGAAGCGTGGCAAAAACGGGCTTTTGCTCACATTTCAAATGGtt                                      <  2:2224661/99‑1 (MQ=255)
                                                           cAAAGATGCATACCCCGATCGGGGTAAAACCTGTAGCAGGATCAAAAGAGTGGCGGGAAGCGTGGCAAAAACGGGGTTTTGCTCACATTTCAAATGGTta                                     >  2:464220/1‑100 (MQ=255)
                                                             aaGATGCATACCCCGATCGGGGTAAAACCTGTAGCAGGATCAAAAGAGTGGCGGGAAGCGTGGCAAAAACGGGCTTTTGCTCACATTTCAAATGGTtata                                   >  1:1202413/1‑100 (MQ=255)
                                                                 tGCATACCCCGATCGGGGTAAAACCTGTAGCAGGTTCAAAAGAGTGGCGGGAAGCGTGGCAAAAACGGGCTTTTGCTCACATTTCAAATGGTTATAAata                               >  1:757758/1‑100 (MQ=255)
                                                                   cATACCCCGATCGGGGTAAAACCTGTAGCAGGTTCAAAAGAGTGGCGGGAAGCGTGGCAAAAACGGGCTTTTGCTCACATTTCAAATGGTTATAAatata                             >  1:502728/1‑100 (MQ=255)
                                                                   cATACCCCGATCGGGGTAAAACCTGTAGCAGGATCAAAAGAGTGGCGGGAAGCGTGGCAAAAACGGGCTTTTGCTCACATTTCAAATGGTTATAAatata                             <  1:2318573/100‑1 (MQ=255)
                                                                   cATACCCCGATCGGGGTAAAACCTGTAGCAGGATCAAAAGAGTGGCGGGAAGCGTGGCAAAAACGGGCTTTTGCTCACATTTCAAATGGTTATAAatata                             >  1:1675252/1‑100 (MQ=255)
                                                                       ccccGATCGGGGTAAAACCTGTAGCAGGTTCAAAAGAGTGGCGGGAAGCGTGGCAAAAACGGGCTTTTGCTCACATTTCAAATGGTTATAAATATATTta                         >  1:655890/1‑100 (MQ=255)
                                                                       ccccGATCGGGGTAAAACCTGTAGCAGGATCAAAAGAGTGGCGGGAAGCGTGGCAAAAACGGGCTTTTGCTCACATTTCAAATGGTTATAAATATATTta                         >  1:607679/1‑100 (MQ=255)
                                                                       ccccGATCGGGGTAAAACCTGTAGCAGGATCAAAAGAGTGGCGGGAAGCGTGGCAAAAACGGGCTTTTGCTCACATTTCAAATGGTTATAAATATATTta                         >  1:2039951/1‑100 (MQ=255)
                                                                        cccGATCGGGGTAAAACCTGTAGCAGGATCAAAAGAGTGGCGGGAAGCGTGGCAAAAACGGGCTTTTGCTCACATTTCAAATGGTTATAAATATATTtat                        <  2:1917396/100‑1 (MQ=255)
                                                                           gATCGGGGTAAAAACTGTAGCAGGATCAAAAAAGTGGCGGGAAGCGTGGCAAAAACGGGCTTTTGCTCACATTTCAAATGGTTATAAATATATTtatata                     >  1:1573197/1‑100 (MQ=255)
                                                                            aTCGGGGTAAAACCTGTAGCAGGATCAAAAGAGTGGCGGGAAGCGTGGCAAAAACGGGCTTTTGCTCACATTTCAAATGGTTATAAATATATTTATATAg                    <  2:287571/100‑1 (MQ=255)
                                                                             tCGGGGTAAAACCTGTAGCAGGATCAAAAGAGTGGCGAAAAAAAAAAAAAAAACGGGCTTTTGCTCACATTTCAAATGGTTATAAATATATTTATATAGc                   <  2:1747782/100‑1 (MQ=255)
                                                                              cGGGGTAAAACCTGTAGCAGGATCAAAAGAGTGGCGGGAAGCGTGGCAAAAACGGGCTTTTGCTCACCTTTCAAATGGTTATAAATATATTTATATAGc                   >  1:2028463/1‑99 (MQ=255)
                                                                              cGGGGTAAAACCTGTAGCAGGATCAAAAGAGTGGCGGGAAGCGTGGCAAAAACGGGCTTTTGCTCACATTTCAAATGGTTATAAATATATTTATATAGCg                  <  1:1819269/100‑1 (MQ=255)
                                                                               ggggTAAAACCTGTAGCAGGATCAAAAGAGTGGCGGGAAGCGTGGCAAAAACGGGCTTTTGCTCACATTTCAAATGGTTATAAATATATTTATATAGCga                 >  1:1230729/1‑100 (MQ=255)
                                                                                 ggTAAAACCTGTAGCAGGATCAAAAGAGTGGCGGGAAGCGTGGCAAAAACGGGCTTTTGCTCACATTTCAAATGGTTATAAATATATTTATATAGCgatt               <  2:1350935/100‑1 (MQ=255)
                                                                                  gTAAAACCTGTAGCAGGATCAAAAGAGTGGCGGGAAGCGTGGCAAAAACGGGCTTTTGCTCACATTTCAAATGGTTATAAATATATTTATATAGCgattg              >  1:1437056/1‑100 (MQ=255)
                                                                                  gTAAAACCTGTAGCAGGATCAAAAGAGTGGCGGGAAGCGTGGCAAAAACGGGCTTTTGCTCACATTTCAAATGGTTATAAATATATTTATATAGCgattg              >  1:813349/1‑100 (MQ=255)
                                                                                    aaaaCCTGTAGCAGGTTCAAAAGAGTGGCGGGAAGCGTGGCAAAAACGGGCTTTTGCTCACATTTCAAATGGTTATAAATATATTTATATAGCgattgat            >  1:1023131/1‑100 (MQ=255)
                                                                                    aaaaCCTGTAGCAGGATCAAAAGAGTGGCGGGAAGCGTGGCAAAAACGGGCTTTTGCTCACATTTCAAATGGTTATAAATATATTTATATAGCgattgat            <  1:1999279/100‑1 (MQ=255)
                                                                                    aaaaCCTGTAGCAGGATCAAAAGAGTGGCGGGAAGCGTGGCAAAAACGGGCTTTTGCTCACATTTCAAATGGTTATAAATATATTTATATAGCgattgat            <  2:53740/100‑1 (MQ=255)
                                                                                    aaaaCCTGTAGCAGGATCAAAAGAGTGGCGGGAAGCGTGGCAAAAACGGGCTTTTGCTCACATTTCAAATGGTTATAAATATATTTATATAGCgattgat            <  1:1285696/100‑1 (MQ=255)
                                                                                     aaaCCTGTAGCAGGATCAAAAGAGTGGCGGGAAGCGTGGCAAAAACGGGCTTTTGCTCACATTTCAAATGGTTATAAATATATTTATATAGCgattgatt           >  2:2079471/1‑100 (MQ=255)
                                                                                     aaaCCTGTAGCAGGATCAAAAGAGTGGCGGGAAGCGTGGCAAAAACGGGCTTTTGCTCACATTTCAAATGGTTATAAATATATTTATATAGCgattgatt           <  1:189220/100‑1 (MQ=255)
                                                                                          tGTAGCAGGTTCAAAAGAGTGGCGGGAAGCGTGGCAAAAACGGGCTTTTGCTCACATTTCAAATGGTTATAAATATATTTATATAGCGATTGATTCACCa      >  1:143714/1‑100 (MQ=255)
                                                                                          tGTAGCAGGTTCAAAAGAGTGGCGGGAAGCGTGGCAAAAACGGGCTTTTGCTCACATTTCAAAGGGTTATAAATATATTTATATAGCGATTGATTCAcc       <  2:269401/99‑1 (MQ=255)
                                                                                          tGTAGCAGGATCAAAAGAGTGGCGGGAAGCGTGGCAAAAACGGGCTTTTGCTCACATTTCAAATGGTTATAAATATATTTATATAGCGATTGATTCACCa      >  2:1993431/1‑100 (MQ=255)
                                                                                           gTAGCAGGATCAAAAGAGTGGCCGGAAGCGTCGCAAATACGGGCTTTTGCTCACATTTCAAATGGTTATAAATATATTTATATAGCGATTGATTCACCag     <  2:1573192/100‑1 (MQ=255)
                                                                                            tAGCAGGATCAAAAGAGTGGCGGGAAGCGTGGCAAAAACGGGCTTTTGCTCACATTTCAAATGGTTATAAATATATTTATATAGCGATTGATTCACCaga    <  2:2039944/100‑1 (MQ=255)
                                                                                             aGCAGGATCAAAAGAGTGGCGGGAAGCGTGGCAAAAACGGGCTTTTGCTCACATTTCAAATGGTTATAAATATATTTATATAGCGATTGATTCACCagag   <  2:2047806/100‑1 (MQ=255)
                                                                                              gCAGGATCAAAAGAGTGGCGGGAAGCGTGGCAAAAACGGGCTTTTGCTCACATTTCAAATGGTTATAAATATATTTATATAGCGATTGATTCACCAgaga  >  1:1284032/1‑100 (MQ=255)
                                                                                                 ngATCAAAAGAGTGGCGGGAAGCGTGGCAAAAACGGGCTTTTGCTCACATTTCAAATGGGTATAAATATATTTATATAGCGATTGATTCACCAgaga  >  1:996117/2‑97 (MQ=255)
                                                                                                   |                                                                                              
AAATAAACCGCATAATTAATAAGCCATTTTTATAGCCGCTAAGATATTAAAGGATGTGTCAAAGATGCATACCCCGATCGGGGTAAAACCTGTAGCAGGATCAAAAGAGTGGCGGGAAGCGTGGCAAAAACGGGCTTTTGCTCACATTTCAAATGGTTATAAATATATTTATATAGCGATTGATTCACCAGAGA  >  NC_000913/2303542‑2303735

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: