Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,449,5550CT20.6% 144.3 / 29.2 68W165* (TGG→TGA) yfcODUF2544 family putative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base C (30/24);  new base T (7/7);  total (37/31)
Fisher's exact test for biased strand distribution p-value = 7.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

TCCCATCTTTCATCAAATACGTTGAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAG  >  NC_000913/2449458‑2449651
                                                                                                 |                                                                                                 
tCCCATCTTTCATCAAATACGTTGAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCa                                                                                                 >  2:175904/1‑100 (MQ=255)
  ccATCTTTCATCAAATACGTTGAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTTCATg                                                                                               >  2:367099/1‑100 (MQ=255)
  ccATCTTTCATCAAATACGTTGAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTTCATg                                                                                               <  2:1530953/100‑1 (MQ=255)
  ccATCTTTCATCAAATACGTTGAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATg                                                                                               <  2:206969/100‑1 (MQ=255)
   cATCTTTCATCAAATACGTTGAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGa                                                                                              >  2:1343159/1‑100 (MQ=255)
      cTTTCATCAAATACGTTGAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAAc                                                                                           >  2:1282823/1‑100 (MQ=255)
      cTTTCATCAAATACGTTGAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAAc                                                                                           <  2:1548728/100‑1 (MQ=255)
         tcatcaAATACGTTGAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTc                                                                                        <  2:894474/100‑1 (MQ=255)
          catcaAATACGTTGAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCa                                                                                       >  2:734482/1‑100 (MQ=255)
             caAATACGTTGAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTTCATGAAACGTCATaa                                                                                    >  2:421861/1‑100 (MQ=255)
             caAATACGTTGAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATaa                                                                                    >  1:277880/1‑100 (MQ=255)
              aaaTACGTTGAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAg                                                                                   >  1:37932/1‑100 (MQ=255)
               aaTACGTTGAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTTCATGAAACGTCATAAgg                                                                                  <  1:1238579/100‑1 (MQ=255)
                aTACGTTGAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGCAACGTCATAAGGt                                                                                 <  1:293124/100‑1 (MQ=255)
                 tACGTTGAAGGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCAGAAGGtt                                                                                >  2:1158814/1‑100 (MQ=255)
                 tACGTTGAAAGCGTTTTCTCAACCGTACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGtt                                                                                <  2:330062/100‑1 (MQ=255)
                 tACGTTGAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTTAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACTTCATAAGGt                                                                                 >  2:522415/1‑99 (MQ=255)
                 tACGTTGAAAGCGTTTTCTCAACCACACCGCCCATCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGtt                                                                                >  2:1023664/1‑100 (MQ=255)
                  aCGTTGAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTg                                                                               <  2:1279484/100‑1 (MQ=255)
                   cGTTGAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTTCATGAAACGTCATAAGGTTGa                                                                              <  1:865119/100‑1 (MQ=255)
                     ttGAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTTCATGAAACGTCATAAGGTTGAAc                                                                            <  1:1275517/100‑1 (MQ=255)
                     ttGAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAAc                                                                            >  1:190960/1‑100 (MQ=255)
                            cGTTTTCTCAACCACACCGCCTGTCGACTGCAGGGTGAGCAGTTTTGACTTATGGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAAc                                                                      <  2:874289/99‑1 (MQ=255)
                              ttttCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGa                                                                   >  1:1363306/1‑100 (MQ=255)
                                ttCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTTCATGAAACGTCATAAGGTTGAACAGAAACGGATg                                                                 >  1:59264/1‑100 (MQ=255)
                                 tctcAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTTCATGAAACGTCATAAGGTTGAACAGAAACGGATgg                                                                >  2:1336475/1‑100 (MQ=255)
                                  ctcAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGa                                                               >  1:729701/1‑100 (MQ=255)
                                  ctcAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGa                                                               >  2:654408/1‑100 (MQ=255)
                                    cAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACa                                                             >  1:848016/1‑100 (MQ=255)
                                    cAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACa                                                             <  1:1282832/100‑1 (MQ=255)
                                    cAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACa                                                             >  2:398166/1‑100 (MQ=255)
                                         acacCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAAtcg                                                        <  2:275699/100‑1 (MQ=255)
                                         acacCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAAtcg                                                        >  2:871713/1‑100 (MQ=255)
                                          cacCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAAtcgt                                                       >  2:244242/1‑100 (MQ=255)
                                             cGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTTCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAAtcgtgtc                                                    <  1:421862/100‑1 (MQ=255)
                                              gCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACTGATGGACAAAtcgtgt                                                     <  2:660714/98‑1 (MQ=255)
                                               cccGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAAtcgtgtcgt                                                  >  2:342762/1‑100 (MQ=255)
                                               cccGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAAtcgtgtcgt                                                  <  1:1248591/100‑1 (MQ=255)
                                                  gTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAc                                               <  2:1055572/100‑1 (MQ=255)
                                                  gTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAc                                               >  1:109581/1‑100 (MQ=255)
                                                     gACTGCAGGGTGCGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TAt                                           <  1:231029/100‑1 (MQ=255)
                                                     gACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TAt                                           >  1:1435763/1‑100 (MQ=255)
                                                      aCTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATc                                          >  1:281726/1‑100 (MQ=255)
                                                       cTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCa                                         >  2:202051/1‑100 (MQ=255)
                                                        tGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAg                                        >  1:400159/1‑100 (MQ=255)
                                                         gCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAgg                                       <  2:1046290/100‑1 (MQ=255)
                                                         gCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAgg                                       <  2:958/100‑1 (MQ=255)
                                                           aGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTTCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGcc                                     >  1:417709/1‑100 (MQ=255)
                                                              gTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCt                                  <  1:342763/100‑1 (MQ=255)
                                                                     gTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGt                           >  2:111886/1‑100 (MQ=255)
                                                                      ttttGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTc                          <  1:283734/100‑1 (MQ=255)
                                                                        ttGACTTACTGACATAATTCTCATTTCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACGTATCAGGGGCCTGATA‑GTCAc                        <  2:996840/100‑1 (MQ=37)
                                                                         tgccTTACTGACATAATTCTCATTTCATGAAACGTCATAAGGTTGAACAGAAACGGATGGATAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCAc                        <  1:811007/96‑1 (MQ=255)
                                                                         tGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAAGCGTGTCGTGAC‑TATCAGGccccagctgag                            >  2:1482341/1‑88 (MQ=255)
                                                                           aCTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGcc                     <  1:1363830/100‑1 (MQ=255)
                                                                              tACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGGACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGccccc                  <  1:76387/100‑1 (MQ=255)
                                                                                  gACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGga              >  2:1505808/1‑100 (MQ=255)
                                                                                   aCATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGgag             >  1:868260/1‑100 (MQ=255)
                                                                                     aTAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGagag           <  1:1531016/100‑1 (MQ=255)
                                                                                       aaTTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGc         <  2:797628/100‑1 (MQ=255)
                                                                                        aTTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCt        <  1:797182/100‑1 (MQ=255)
                                                                                          tctcATTTCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTc      >  2:105267/1‑100 (MQ=255)
                                                                                          tctcATTTCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTc      >  2:105190/1‑100 (MQ=255)
                                                                                           ctcATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCa     <  2:109581/100‑1 (MQ=255)
                                                                                             cATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCaga   >  1:305167/1‑100 (MQ=255)
                                                                                             cATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCaga   >  1:1072341/1‑100 (MQ=255)
                                                                                              aTTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCATGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCagag  >  1:1340000/1‑100 (MQ=255)
                                                                                              aTTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCagag  <  1:1394823/100‑1 (MQ=255)
                                                                                                 |                                                                                                 
TCCCATCTTTCATCAAATACGTTGAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAG  >  NC_000913/2449458‑2449651

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: ‑