Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,604,3440CT35.1% 73.8 / 30.9 57S62F (TCT→TTT) bglJDNA‑binding transcriptional regulator BglJ
Reads supporting (aligned to +/- strand):  ref base C (6/31);  new base T (20/0);  total (26/31)
Fisher's exact test for biased strand distribution p-value = 1.90e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCAT  >  NC_000913/4604208‑4604481
                                                                                                                                        |                                                                                                                                         
tGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttcttt                                                                                                                                        <  1:1547270/140‑1 (MQ=255)
    tCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTttttctttttct                                                                                                                                    <  2:786937/140‑1 (MQ=255)
    tCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTttttctttttct                                                                                                                                    <  2:682690/140‑1 (MQ=255)
    tCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTttttctttttct                                                                                                                                    <  1:1548422/140‑1 (MQ=255)
    tCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTttttctttttct                                                                                                                                    <  2:159935/140‑1 (MQ=255)
      aTGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGc                                                                                                                                  <  2:1420001/140‑1 (MQ=255)
              taTCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGaa                                                                                                                          >  1:451586/1‑140 (MQ=255)
                tCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGt                                                                                                                        <  1:1360905/140‑1 (MQ=255)
                tCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGt                                                                                                                        <  1:1403528/140‑1 (MQ=255)
                   gTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGTggg                                                                                                                     >  2:1270822/1‑138 (MQ=255)
                   gTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAggggg                                                                                                                     >  1:903590/1‑136 (MQ=255)
                   gTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTTTGCCATGAGAAgaggg                                                                                                                     >  2:1232767/1‑136 (MQ=255)
                    tATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGTGAgg                                                                                                                    >  2:829960/1‑139 (MQ=255)
                     aTTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCACGGAAAGTGAggg                                                                                                                   >  2:1164908/1‑138 (MQ=255)
                     aTTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAgcg                                                                                                                   <  1:1218441/140‑1 (MQ=255)
                     aTTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAGTTTTTTTTTTTTTGCCCTGAGAAAAGAgcg                                                                                                                   >  1:1490564/1‑140 (MQ=255)
                       tGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTTTACCAGTCAGGAGTGATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCa                                                                                                                 <  2:553515/140‑1 (MQ=255)
                        gagaGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCag                                                                                                                <  2:1250515/140‑1 (MQ=255)
                        gagaGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCag                                                                                                                <  1:907070/140‑1 (MQ=255)
                         agagTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCaga                                                                                                               <  2:1121013/140‑1 (MQ=255)
                          gagTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTCCCCTGAGAAGAGAGCGCagag                                                                                                              >  1:1537130/1‑140 (MQ=255)
                              ttATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAATGAGCGCAGAGAggg                                                                                                          >  2:1218441/1‑140 (MQ=255)
                              ttATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAggg                                                                                                          <  2:781999/140‑1 (MQ=255)
                              ttATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAggg                                                                                                          <  2:518863/140‑1 (MQ=255)
                                aTTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGGGAGCGCagggaggggt                                                                                                        >  1:594092/1‑138 (MQ=255)
                                       aaGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTg                                                                                                 <  2:52704/140‑1 (MQ=255)
                                                         tataAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGAt                                                                               <  1:1279489/140‑1 (MQ=255)
                                                              gCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGt                                                                          <  1:1033834/140‑1 (MQ=255)
                                                               cTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTCTTTTTTTTCCCTGCGGAGAGAGCGCAAAGAGGGGATATCTTGTCTGTCTTAACCAGCGATTaaact                                                                         >  2:756329/1‑137 (MQ=255)
                                                               cTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGtt                                                                         <  1:1270980/140‑1 (MQ=255)
                                                                tCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGttt                                                                        >  2:1280590/1‑140 (MQ=255)
                                                                tCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCCTGAGAAGTGAGCGCAGAGAGGGGTTATCTTTCCTGACTGAACTGGCGATTAAGttt                                                                        >  2:1279489/1‑140 (MQ=255)
                                                                tCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCCTTTTTTCTTTTTCTTCCACGAGAAGAGAGCGCAGAGAGGGATTATCTTTCCTGACTGAACTTGCGATTaaaatt                                                                        >  2:132949/1‑136 (MQ=255)
                                                                         ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGTGAGCGCAGAGGGGGGTTATCTTTCCTGACTGAACTGGCGATTAAGTTTCCGCGTAcc                                                               >  2:1556235/1‑140 (MQ=255)
                                                                         ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTAcc                                                               <  2:1212320/140‑1 (MQ=255)
                                                                          taCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTAccc                                                              <  1:740285/140‑1 (MQ=255)
                                                                          taCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTAccc                                                              <  2:451586/140‑1 (MQ=255)
                                                                            ccAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCgg                                                            <  1:1328100/140‑1 (MQ=255)
                                                                                tCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTTAGAAATGGGCGCAGGGAGGGGGTTTTTTTCCTGACTGAACTGGCGGTTTAATTTTCGCGTACCCGGCggt                                                        >  2:765128/1‑138 (MQ=255)
                                                                                tCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGAAGGGAGCGCAGAGGGGGGTTATTTTGCCTGACTGAACTGGCGATTAAGTTTTCGCGGACCCGGCggt                                                        >  2:190154/1‑138 (MQ=255)
                                                                                tCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTCCCATGAGAAGTGAGCGCAGAGAGGGGTTATCTTGCCTGACTGAACTGGGGGTTAAGTTTTCCCGTACCCcgcggt                                                        >  2:1328100/1‑138 (MQ=255)
                                                                                 cAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGAAGTGAGCGCAGAGAGGGGTTATCTTGCCTGGCTTAACTGGCGGTTAAGTTTTCGCGTACCCCGCGgtt                                                       >  2:718775/1‑140 (MQ=255)
                                                                                         aTTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGGGAGCGCAGAGGGGGGTTATCTTGCCTGACTGAACTGGCGATTTAAGTTCCGCGTACCCGGCGTTTAGTTAttt                                               >  2:36795/1‑139 (MQ=255)
                                                                                         aTTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTTTTCCCTTAGAAGGGAGGGGAGGGAGGGGTTATCTTTCCTGACTTAACTGGGGGTTAAGGTTTCGCGGACCCCGCGGTTTGTttttt                                               >  2:556414/1‑139 (MQ=255)
                                                                                           ttCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCg                                             <  1:757105/140‑1 (MQ=255)
                                                                                             cAGGATGCCATGTCGCGGATCTCGGTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGga                                           <  1:1357463/140‑1 (MQ=255)
                                                                                              aGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGgat                                          <  1:190154/140‑1 (MQ=255)
                                                                                              aGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGgat                                          <  1:756329/140‑1 (MQ=255)
                                                                                                gATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGgatga                                        <  1:986630/140‑1 (MQ=255)
                                                                                                 aTGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGgatgat                                       <  1:718775/140‑1 (MQ=255)
                                                                                                 aTGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGgatgat                                       <  1:559069/140‑1 (MQ=255)
                                                                                                     cATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGAAATGAGCGCCGAGAGGGGTTATTTTTTCTGACTGAACTGGCGATTAAGTTTTCGCGTACCCGGCGGTTTGTTTTTTCGGGTGATGata                                   >  2:1146301/1‑140 (MQ=255)
                                                                                                     cATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTTTTTCCTGTGAAGTGAGCGGGGCGGGGGGNTTTCTNTCCTGGCTGAACTGGGGGTTAAAGTTTCGCGGACCCGGCGGTTTGTTATTGCGggtggtgatt                                   >  2:824069/1‑139 (MQ=255)
                                                                                                       tGTCGCGGATCTCGTTTGCGGCGGTCTTTTTTTCTTTTTTTGTCCTGAGAAGGTAGTGCGGAGGGGGGTTATATTGTCTGTCTGAAGTGGCGGTTAAGGTTCCGCGTAACCGGGGGTTTGTTATTGCGggtgatgtgatt                                 >  2:803099/1‑135 (MQ=255)
                                                                                                       tGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGAAGGGAGCGCAGGGGGGGGTTATCTTGCCTGACTTAACTGGCGGTTAAGTTTTCGCGTAACCGGCGTTTAGTTATTGCGGATGGTGGTAtt                                 >  2:802790/1‑140 (MQ=255)
                                                                                                          cgcgGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTTCCATGAGAAGTGAGCGCAGAGAGGGGTTATTTTGCCTGACTGAACTGGCGATTTAGTTTCCGCGTACCCCGCGTTTAGTTATTTCGGATGATGGTATTTaa                              >  2:559069/1‑140 (MQ=255)
                                                                                                           gcgGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAg                             <  1:803099/140‑1 (MQ=255)
                                                                                                           gcgGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAg                             <  1:1270822/140‑1 (MQ=255)
                                                                                                           gcgGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAg                             <  1:230267/140‑1 (MQ=255)
                                                                                                              gATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGAAGGGAGCGCAGAGGGGGGTTTTCTTGCCTGACTGAACTGGCGGTTAAGTTTTCGCGTACCCCGCGGTTTGTTTTTTCGGGTGATGATATTTAAGCTc                          >  2:692793/1‑140 (MQ=255)
                                                                                                               aTCTCGTTTGCGGCGGTCATTTTTTCTTTTTTTTCCCTGAGAAGTGAGCGCGGAGAGGGGTTATCTTGCCTGACTGAACTGGGGGTTAAGTTTTCGCGTACCCGGCGTTTTGTTATTTTGGATGATGATATTTAAGCTCg                         >  2:1121834/1‑140 (MQ=255)
                                                                                                                    gTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGa                    <  1:824069/140‑1 (MQ=255)
                                                                                                                        gcggcgGTCATTTTTTTTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGGTTATTTTGCCTGACTGAAATGGCGGTTAAGTTTCCGCGTACCCGGCGGTTTGTTTTTTCGGGTGGTGATATTGAAGGTCGGGTGGTTgg                >  2:1127188/1‑140 (MQ=255)
                                                                                                                          ggcggTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGTGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGtt              <  2:1490564/140‑1 (MQ=255)
                                                                                                                               tCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTg         <  1:1164908/140‑1 (MQ=255)
                                                                                                                                      ttctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCAt  >  1:995055/1‑140 (MQ=255)
                                                                                                                                      ttctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCAt  <  1:802790/140‑1 (MQ=255)
                                                                                                                                      ttctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCAt  <  1:556414/140‑1 (MQ=255)
                                                                                                                                      ttctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCAt  <  1:36795/140‑1 (MQ=255)
                                                                                                                                      ttctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCAt  <  2:995055/140‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                         
TGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCAT  >  NC_000913/4604208‑4604481

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑