Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL606273,5810TG31.5% 146.2 / 19.7 73G136G (GGT→GGGyagUconserved inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (14/36);  new base G (23/0);  total (37/36)
Fisher's exact test for biased strand distribution p-value = 1.48e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AGTGTTTGCTGTCGGTTATTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCATGGGACTGACACCACCTCTGTTTGATCTTCCG  >  REL606/273506‑273656
                                                                           |                                                                           
aGTGTTTGCTGTCGGTTATTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGt                                                                             <  2:878870/76‑1 (MQ=255)
 gtgtTTGCTGTCGGTTATTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTg                                                                            <  2:1468337/76‑1 (MQ=255)
  tgtTTGCTGTCGGTTATTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGc                                                                           >  2:1092214/1‑76 (MQ=255)
     ttGCTGTCGGTTATTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCttt                                                                        <  2:1183938/76‑1 (MQ=255)
     ttGCTGTCGGTTATTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGCgtt                                                                        >  1:908703/1‑76 (MQ=255)
      tGCTGTCGGTTATTCTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTa                                                                       <  2:622278/76‑1 (MQ=255)
       gCTGTCGGTTATTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAg                                                                      <  1:923380/76‑1 (MQ=255)
       gCTGTCGGTTATTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAg                                                                      <  2:1511647/76‑1 (MQ=255)
       gCTGTCGGTTATTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAg                                                                      <  2:1265638/76‑1 (MQ=255)
       gCTGTCGGTTATTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGCgtttg                                                                      >  2:370833/1‑74 (MQ=255)
        cTGTCGGTTATTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGGtttttc                                                                     >  2:790744/1‑73 (MQ=255)
        cTGTCGGTTATTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGc                                                                     <  1:1148935/76‑1 (MQ=255)
         tGTCGGTTATTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGGTTTTGcc                                                                    >  2:1271900/1‑76 (MQ=255)
           tCGGTTATTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCa                                                                  <  1:883054/76‑1 (MQ=255)
           tCGGTTATTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGCTTTAGCCCa                                                                  >  1:808844/1‑76 (MQ=255)
            cGGTTATTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCaa                                                                 <  2:893285/76‑1 (MQ=255)
            cGGTTATTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGCTTTAGCCCaa                                                                 >  2:637797/1‑76 (MQ=255)
             ggTTATTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGCTTTTGCCCCAc                                                                >  2:287992/1‑76 (MQ=255)
             ggTTATTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGCTTTTGCCCCAc                                                                >  2:947207/1‑76 (MQ=255)
               ttATTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGCTTTTGCCCAACtt                                                              >  1:995024/1‑76 (MQ=255)
                 aTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGCTTTAGCCCAACcttt                                                            >  1:287023/1‑76 (MQ=255)
                  ttGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGCTTTAGCCCAACttttt                                                           >  1:614092/1‑76 (MQ=255)
                   tgtgtgGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGCTTTTGCCCAACtttttt                                                          >  2:1551012/1‑75 (MQ=255)
                     tgtgGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGCTTTTGCCCAACTTTTTGtt                                                        >  1:226606/1‑76 (MQ=255)
                     tgtgGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGCTTTTGCCCAACTTTTTGtt                                                        >  2:378888/1‑76 (MQ=255)
                      gtgGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTTGCCCACCCTttttttt                                                       >  2:1149387/1‑75 (MQ=255)
                      gtgGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGCTTTTGCCCAACttttttttt                                                       >  2:586483/1‑75 (MQ=255)
                       tgGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGCTTTAGCCCAACTTTTTTTTCa                                                      >  2:441176/1‑76 (MQ=255)
                        ggTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGCTTTTGCCCCACttttttttttt                                                     >  1:996446/1‑73 (MQ=255)
                         gTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTTGCCCCACTTTTTTTTtata                                                    >  1:337007/1‑76 (MQ=255)
                           cGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTTGCAGGTGCTTTAGCCCAACTTTTTGTTCATATg                                                  <  2:1460640/76‑1 (MQ=255)
                           cGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATg                                                  <  1:342177/76‑1 (MQ=255)
                             ctcaaGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGCTTTTGCCCAACTTTTTTTTTATATGAt                                                >  1:278133/4‑76 (MQ=255)
                               gAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGCTTTAGCCCAACTTTTTTTTCATATGAttt                                              >  2:316026/1‑76 (MQ=255)
                                aaGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGAtttc                                             >  1:202632/1‑76 (MQ=255)
                                aaGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGAtttc                                             <  2:663422/76‑1 (MQ=255)
                                   tATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTGATATGAtttcttt                                          <  1:704367/76‑1 (MQ=255)
                                   tATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGAtttcttt                                          <  2:514330/76‑1 (MQ=255)
                                   tATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGAtttcttt                                          <  2:1223563/76‑1 (MQ=255)
                                   tATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGAtttcttt                                          <  2:889149/76‑1 (MQ=255)
                                    aTTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTTGCCCAACTTTTTTTTCATATGAtttctttc                                         >  2:1631967/1‑76 (MQ=255)
                                     tttCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTcc                                        <  2:1142670/76‑1 (MQ=255)
                                        ccAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGCTTTAGCCCAACTTTTTTTTTATATGATTTCTTTTCctc                                     >  1:19406/1‑76 (MQ=255)
                                         cAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCctct                                    <  1:673724/76‑1 (MQ=255)
                                          aaaaaTTAAAATCTTGCAGGGCTTACTGGCAGGGGGGTTTGCCCCACTCTTTTTTTTTTTTATTTTTTTTCCtctc                                   >  2:1581301/1‑76 (MQ=255)
                                           aaaaTTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCa                                  <  1:896078/76‑1 (MQ=255)
                                            aaaTTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCAt                                 <  2:1387989/76‑1 (MQ=255)
                                             aaTTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCATg                                <  1:1635763/76‑1 (MQ=255)
                                             aaTTAAACTCTGGCAGGGCTTACTGGCAGGGGCTTTTGCCCCACTTTTTTTTTATATTATTTTTTTTCCTCTCATg                                >  2:1125493/1‑76 (MQ=255)
                                             aaTTAAACTCTGGCAGGGCTTACTGGCAGGGGCTTTAGCCCAACTTTTTTTTTATATGATTTTTTTTCCTCTCATg                                >  2:1505449/1‑76 (MQ=255)
                                             aaTTAAACTCTGGCAGGGCTTACTGGCAGGGGCTTTAGCCCAACTTTTTTTTCATATGATTTCTTTCCCTCTCCTg                                >  1:227753/1‑76 (MQ=255)
                                              aTTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCATgg                               <  2:1633191/76‑1 (MQ=255)
                                              aTTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCATgg                               <  2:1595735/76‑1 (MQ=255)
                                              aTTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCATgg                               <  1:1508857/76‑1 (MQ=255)
                                               ttAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTTTTCATATGATTTCTTTCCCTCTCATggg                              >  1:184816/1‑76 (MQ=255)
                                                 aaaCTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCCACTTTTTTTTTATATGATTTTTTTTCCTCTCCTggggc                            >  1:328610/1‑74 (MQ=255)
                                                 aaaCTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCCACTTTTTTTTTATATGATTTTTTTTCCTCTCATGGGAc                            >  2:1279502/1‑76 (MQ=255)
                                                 aaaCTCTGGCAGGGCTTACTGGCAGGGGCTTTTGCCCAACTTTTTTTTTATATGATTTTTTTTCCTCTCATggggc                            >  1:436204/1‑74 (MQ=255)
                                                   aCTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCATGGGACTg                          <  1:116912/76‑1 (MQ=255)
                                                   aCTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCATGGGACTg                          <  1:1166874/76‑1 (MQ=255)
                                                    ctctGGCAGGGCTTACTGGCAGGTGCTTTAGGCTCACCCTTTGTTTATATTGTTTCTTTTCCTCTCCTTGGAccgc                         >  1:660488/1‑73 (MQ=255)
                                                    ctctGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCATGGGACTGa                         <  2:704541/76‑1 (MQ=255)
                                                    ctctGGCAGGGCTTACTGGCAGGGGCTTTTGCCCCACTTTTTTTTTATATTATTTTTTTTCCTCTCCTTGGGCTGa                         >  1:366311/1‑76 (MQ=255)
                                                    ctctGGCAGGGCTTACTGGCAGGGGCTTTTGCCCCACTTTTTTTTTATATGATTTTTTTTCCCCTCCTGGGACTGa                         >  1:1332176/1‑76 (MQ=255)
                                                    ctctGGCAGGGCTTACTGGCAGGGGCTTTTGCCCAACTTTTTTTTTATATGATTTTTTTTCCTCTCATTGGGCTGa                         >  1:660472/1‑76 (MQ=255)
                                                      ctGGCAGGGCTTACTGGCAGGGGCTTTAGCCCAACTTTTTGTTCATATGATTTTTTTCCCTCTCTTGGGGCTGaca                       >  2:1099253/1‑76 (MQ=255)
                                                        ggCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCATGGGACTGAcacc                     <  1:106843/76‑1 (MQ=255)
                                                         gCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCATGGGACTGAcacca                    <  1:441176/76‑1 (MQ=255)
                                                           aGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCATGGGACTGACaccccc                  >  2:1554830/1‑76 (MQ=255)
                                                           aGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCATGGGACTGACaccacc                  <  2:1206428/76‑1 (MQ=255)
                                                            gggCTTACTGGCAGGGGCTTTTGCCCAACTTTTTTTTTATATGATTTTTTTTCCTCTCATGGGGCTGACaccaccc                 >  2:710908/1‑75 (MQ=255)
                                                            gggCTTACTGGCAGGAGCTTTTGACCCACTTTTTGTTTCTTTTTTTTTTTTCCCTCTCTTGGGAGCGACacccccc                 >  2:1146686/1‑75 (MQ=255)
                                                              gCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCGCTCATGGGACTGACACCACctct               <  1:41823/76‑1 (MQ=255)
                                                              gCTTACTGGCAGGGGCTTTAGCCCCACTTTTTTTTTATATGATTTCTTTCCCTCTCATGGGACTGACACCACctct               >  2:129756/1‑76 (MQ=255)
                                                               cTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCATGGGACTGACACCTCCTCTg              <  2:1282332/76‑1 (MQ=255)
                                                                ttACTGGCAGGGGCTTTTGCCCAACTTTTTTTTCATATGATTTCTTTCCCTCTCATGGGACTGACACCACCCCTGt             >  2:1280733/1‑76 (MQ=255)
                                                                  aCTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCATGGGACTGACACCACCTCTGttt           <  2:272197/76‑1 (MQ=255)
                                                                  aCTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCATGGGACTGACACCACCTCTGttt           <  1:1643152/76‑1 (MQ=255)
                                                                   cTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCATGGGACTGACACCACCTCTGTTTg          >  2:105071/1‑76 (MQ=255)
                                                                    tGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCATGGGACTGACACCACCTCTGTTTGa         >  2:1666219/1‑76 (MQ=255)
                                                                     ggCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCATGGGACTGACACCACCTCTGTTTGAt        >  1:1343767/1‑76 (MQ=255)
                                                                     ggCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCATGGGACTGACACCACCTCTGTTTGAt        <  2:1484679/76‑1 (MQ=255)
                                                                     ggCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCATGGGACTGACACCACCTCTGTTTGAt        >  1:701063/1‑76 (MQ=255)
                                                                      gCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCATGGGACTGACACCACCTCTGTTTGATc       >  2:815949/1‑76 (MQ=255)
                                                                       cAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCATGGGACTGACACCACCTCTGTTTGATCt      >  2:1423739/1‑76 (MQ=255)
                                                                        aGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCATGGGACTGACACCACCTCTGTTTGATCtt     >  1:1671051/1‑76 (MQ=255)
                                                                         ggTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCATGGGACTGACACCACCTCTGTTTGATCTTc    <  1:751540/76‑1 (MQ=255)
                                                                           tGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCATGGGACTGACACCACCTCTGTTTGATCTTCCg  <  2:695162/76‑1 (MQ=255)
                                                                           |                                                                           
AGTGTTTGCTGTCGGTTATTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCATGGGACTGACACCACCTCTGTTTGATCTTCCG  >  REL606/273506‑273656

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: