Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,049,885 T→C intergenic (+134/‑13) polA → / → spf 5' to 3' DNA polymerase and 3' to 5'/5' to 3' exonuclease/Spot 42 sRNA antisense regulator of galK translation, Hfq‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,049,8850TC97.4% 229.5 / ‑4.3 78intergenic (+134/‑13)polA/spf5' to 3' DNA polymerase and 3' to 5'/5' to 3' exonuclease/Spot 42 sRNA antisense regulator of galK translation, Hfq‑dependent
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base C (33/41);  total (34/44)
Fisher's exact test for biased strand distribution p-value = 5.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

CTAAGATTCGCCTGAACATGCCTTTTTTCGTAAGTAAGCAACATAAGCTGTCACGTTTTGTGATGGCTATTAGAAATTCCTATGCAACAACTG‑AAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAA  >  NC_000913/4049748‑4050021
                                                                                                                                          |                                                                                                                                        
cTAAGATTCGCCTGAACATTCCTTTTTTCGTAAGTAAGCAACATAAGCTGTCACGTTTTGTGATGGCTATTAGAAATTCCTATGCAACAACTG‑AAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGCAAAGt                                                                                                                                     <  1:6277/143‑1 (MQ=255)
 tAAGATTCGCCTGAACATGCCTTTTTTCGTAAGTAAGCAACATAAGCTGTCACGTTTTGTGATGGCTATTAGAAATTCCTATGCAACAACTG‑AAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGCAAAGtt                                                                                                                                    >  2:550127/1‑143 (MQ=255)
   aGATTCGCCTGAACATGCCTTTTTTCGTAAGTAAGCAACATAAGCTGTCACGTTTTGTGATGGCTATTAGAAATTCCTATGCAACAACTG‑AAAAAAAATTACACAAAGTGCTTTCTGAACTGAACAAAAAAGAGAAAAGTCAg                                                                                                                                  >  1:358632/1‑143 (MQ=255)
         gCCTGAACATGCCTTTTTTCGTAAGTAAGCAACATAAGCTGTCACGTTTTGTGATGGCTATTAGAAATTCCTATGCAACAACTGAAAAAAAAATTACACAAAGTGCTTTCTTAACTGAACAAAACAGAGCAAACTTGGTcgcg                                                                                                                             >  1:18323/1‑143 (MQ=255)
         gCCTGAACATGCCTTTTTTCGTAAGTAAGCAACATAAGCTGTCACGTTTTGTGATGGCTATTAGAAATTCCTATGCAACAACTG‑AAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGCAAAGTTAGTCGCGt                                                                                                                            >  2:217780/1‑143 (MQ=255)
                   gCCTTTTTTCGTAAGTAAGCAACATAAGCTGTCACGTTTTGTGATGGCTATTAGAAATTCCTATGCAACAACTG‑AAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACaga                                                                                                                  >  1:1309724/1‑143 (MQ=255)
                   gCCTTTTTTCGTAAGTAAGCAACATAAGCTGTCACGTTTTGTGATGGCTATTAGAAATTCCTATGCAACAACTG‑AAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACaga                                                                                                                  >  1:1309944/1‑143 (MQ=255)
                   gCCTTTTTTCGTAAGTAAGCAACATAAGCTGTCACGTTTTGTGATGGCTATTAGAAATTCCTATGCAACAACTG‑AAAAAAAATTACAAAAAGTGCTGTCTGAACTGAACAAAAAAGAGCAAAGTTAGTAGCGCAGGGTACaga                                                                                                                  >  2:821609/1‑143 (MQ=255)
                   gCCTTCTCTCGTAAGTAAGCAACATAACCTGTCACGTTTTGTGATGGCTATTAGAAATTCCTATGCAACAACTG‑AAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACaga                                                                                                                  <  2:594466/143‑1 (MQ=255)
                    ccTTTTTTCGTAAGTAAGCAACAAAAGCAGTCACATACTATGATCGCTATTAGAAATTCCTATGCAACAACTG‑AAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACagag                                                                                                                 <  2:196015/143‑1 (MQ=255)
                      ttttttCGTAAGTAAGCAACATAAGCTGTCACGTTTTGTGATGGCTATTAGAAATTCCTATGCAACAACTG‑AAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACAGAGGt                                                                                                               >  2:1287247/1‑143 (MQ=255)
                         tttCGTAAGTAAGCAACTTAAGCTGTCACGGTTTGTGATGGCTATTAGAAATTCCTATGCAACAACTG‑AAAAAAAATTACAAAAAGTGCTTTCTTAACTGAACAAAAAAGAGCAAAGTTAGTCTCGTAGGGTACAGAGGTAAg                                                                                                            <  1:353856/143‑1 (MQ=255)
                           tCGTAAGTAAGCAACATAAGCTGTCACGTATTGTGATGGCTATTAGAAATTCCAATGCAACAACTG‑AAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGCAAAGGTAGACGCGTGGGGTGCAGAGGTAAGAt                                                                                                          >  1:462770/1‑143 (MQ=255)
                            cgtatttaataaaCATAAGCTGTCAAGTTTTGTAATGGCGATTAAAAATTCCTATTCAACAACTG‑AAAAAAAAGTACAAAAAATACCTGCTGAACAGAACAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATg                                                                                                         <  2:1344202/132‑1 (MQ=255)
                             gtaagtaaGCAACATAAGCTGTCACGTTTTGTGATGGCTATTAGAAATTCCTATGCAACAACTG‑AAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGCAAAGTTCGTCGCGTAGGGTACAGAGGTAAGATGt                                                                                                        >  2:949665/1‑143 (MQ=255)
                             gtaagtaaGCAACATAAGCTGGCACGTTTTGTGATGGCTATTAGAAATTCCTATGCAACAACTG‑AAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGt                                                                                                        <  2:970645/143‑1 (MQ=255)
                              gaagtaagGAACATAAGCTGTGACGGTTTGTGATGGCTATTAGAAATTCCTATGCAACAACTG‑AAAAAAAAGTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGtt                                                                                                       <  1:916495/142‑1 (MQ=255)
                               aagtaagCAACATAAGCTGTCACGTTTTGTGATGGCTATTAGAAATTCCTATGCAACAACTG‑AAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTc                                                                                                      <  2:713840/143‑1 (MQ=255)
                               aagtaagCAACATAAGCTGTCACGTTTTGTGATGGCTATTAGAAATTCCTATGCAACAACTG‑AAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGCAAACGTAGTCGCGTCAGCCACCGAGGAAAGATGTTc                                                                                                      >  1:157323/1‑143 (MQ=255)
                                agtaagCAACATAAGCTGTCACGTTTTGTGATGGCTATTAGAAATTCCTATGCAACAACTG‑AAAAAAAATTACAAAATGTGCTTTCTGAACTGAACAAAAAAGGGCAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCt                                                                                                     <  1:345760/143‑1 (MQ=255)
                                     gCAACATAAGCTGTCACGTTTTGTGATGGCTATTAGAAATTCCTATGCAACAACTG‑AAAAGAAGTTACAACAAGTGCTTTCTGAACTGAACACCAAAGCGCACAGTTGGTCGCGTAGGGTACAGCGGTGTGATGTTCtcgctt                                                                                                >  2:189462/1‑138 (MQ=255)
                                     gCAACATAAGCTGTCACGTTTTGTGATGGCTATTAGAAATTCCTATGCAACAACTG‑AAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTCCtctcct                                                                                                >  2:184726/1‑141 (MQ=255)
                                           tAAGCTGCACCGTTTTGTGATGGCTATTAGAAATTCCTATGCAACAACTG‑AAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGAc                                                                                          <  1:776232/143‑1 (MQ=255)
                                            aaGCTGTCACGTTTTGTGATGGCTATTAGAAATTCCTATGCAACAACTG‑AAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGAcc                                                                                         >  2:829433/1‑143 (MQ=255)
                                               cTGTCACGTTTTGTGATGGCTATTAGAAATTCCTATTCAATAACTG‑AAAAAAAAGTACAAATAGTGCTTTCTGAACTGAACAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCttt                                                                                      <  1:1085229/143‑1 (MQ=255)
                                               cTGTCACGTTTTGTGATGGCTATTAGAAATTCCTATGCAACAACTG‑AAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACCGAGGTACGATGTTCTATCTTTCAGACCttt                                                                                      >  2:1132815/1‑143 (MQ=255)
                                               cTGTCACGTTTTGTGATGGCTATTAGAAATTCCTATGCAACAACTG‑AAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCttt                                                                                      <  1:993568/143‑1 (MQ=255)
                                               cTGTCACGTTTTGTGATGGCTATTAGAAATTCCTATGCAACAACTG‑AAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCttt                                                                                      <  2:401259/143‑1 (MQ=255)
                                                 gTCACGTTTTGTGATGGCTATTAGAAATTCCTATGCAACAACTG‑AAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTa                                                                                    <  1:350177/143‑1 (MQ=255)
                                                 gTCACGTTATGAAGTGGCTATTAGAAATTCCTATGCAAAAACTG‑AAAAAAAATGACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTa                                                                                    <  1:974250/143‑1 (MQ=255)
                                                  tCACGTTTTGTGATGGCTATTAGAAATTCCTATGCAACAACTG‑AAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTGAc                                                                                   >  1:1313516/1‑143 (MQ=255)
                                                   cACGTATTGTGATGGCTATGAGAAATTCCGAGGCAACAACTG‑AAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGCCAAGGTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACt                                                                                  <  2:1354330/143‑1 (MQ=255)
                                                     cggtgtgggATGGCGATTAGAAGTTCGGATGCAACAGGTG‑AAAAAAAATTACAAAAAGTGCTGGCTGAACTGAACAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTc                                                                                <  1:924391/135‑1 (MQ=255)
                                                     cGTTTTGTGATGGCTATTAGAAATTCCTATGCAACAACTG‑AAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTc                                                                                <  1:311825/143‑1 (MQ=255)
                                                       ttttGTGATGGCTATTAGAAATTCCTATGCAACAACTG‑AAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCAc                                                                              <  1:217780/143‑1 (MQ=255)
                                                        tttGTGATGGCTATTAGAAATTCCTATGCAACAACTG‑AAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACg                                                                             <  1:193396/143‑1 (MQ=255)
                                                        tttGTGATGGCTATTAGAAATTCCTATGCAACAACTG‑AAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACg                                                                             >  1:586115/1‑143 (MQ=255)
                                                         ttGTGATGGCTATTAGACATTCCTATGAAACAACAG‑AAAAAAAATTACAAAAAGTGCCTTCTGAACTGAACAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGt                                                                            <  1:1296583/143‑1 (MQ=255)
                                                           gtgATGGCTATTAGAAATTCCTATGCAACAACTG‑AAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGCAAAGTTAGTCGCGTAGGGAAAAGCGGTAAGACGTTCTATCTTTCAGACCTTTTACTTCACGTaa                                                                          >  1:235316/1‑143 (MQ=255)
                                                            tgATGGCTATTAGAAATTCCTATGCAACAACTG‑AAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAAt                                                                         <  1:829433/143‑1 (MQ=255)
                                                                 gCTATTAGAAATTCCTATGCAAAAAATG‑AAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGAt                                                                    <  1:118411/143‑1 (MQ=255)
                                                                   tATTAGAAATTCCTATGCAACAACTG‑AAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGAttt                                                                  <  1:231224/143‑1 (MQ=255)
                                                                     ttAGAAATTCCTATGCAACAACTG‑AAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTgg                                                                <  2:1329584/143‑1 (MQ=255)
                                                                              ccTATGCAACAACTG‑AAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATAtt                                                       <  2:1309944/143‑1 (MQ=255)
                                                                              ccTATGCAACAACTG‑AAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATAtt                                                       <  2:1309724/143‑1 (MQ=255)
                                                                               cgaTGCAACAACTG‑AAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATAttt                                                      <  1:1297857/141‑1 (MQ=255)
                                                                                   gCAACAACTG‑AAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGCAAAGTCAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCCGACCTTTTACTTCACGCAATCGGATTTCGCTCAATATTTTAgc                                                  >  1:846355/1‑143 (MQ=255)
                                                                                    caacgaCTG‑AAAAAAAATTACAAAAAGTGCTTGCTGAACTGAACAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAgcc                                                 <  1:322017/143‑1 (MQ=255)
                                                                                    caacaaCTG‑AAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCCCGTAATCGGATTTGGCTGCATATTTTAgcc                                                 >  2:1307942/1‑143 (MQ=255)
                                                                                     aacaacTG‑AAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAgccg                                                >  2:932385/1‑143 (MQ=255)
                                                                                        aacTG‑AAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGccc                                             >  1:593557/1‑143 (MQ=255)
                                                                                        aacTG‑AAAAAAAAGTACAAGAAGTGCTTTCTGAACTGAACAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGccc                                             <  1:1202584/143‑1 (MQ=255)
                                                                                           tG‑AAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATATTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATACTTTAGCCGCCCcag                                          >  1:750472/1‑143 (MQ=255)
                                                                                               aaaaaaaTTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCcagtca                                       <  2:715366/143‑1 (MQ=255)
                                                                                               aaaaaaaTTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCcagtca                                       >  2:825045/1‑143 (MQ=255)
                                                                                                aaaaaaTTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATTTTAAACACCCCCCACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCcagtcag                                      >  2:329596/1‑143 (MQ=255)
                                                                                                aaaaaaTTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCcagtcag                                      >  2:340133/1‑143 (MQ=255)
                                                                                                  aaaaTTACGAAAAGTTCTGTCCGAACTGAACAAAAAAGAGCAAAGTTAGACGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTCATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTa                                    <  1:84034/143‑1 (MQ=255)
                                                                                                   aaaTTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTaa                                   <  2:1260918/143‑1 (MQ=255)
                                                                                                      ttACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGa                                <  2:352675/143‑1 (MQ=255)
                                                                                                         ccataattGCGTTCTGAACTGAACAAAACACAGCAGAGTTAGTCGCGCAGGGTACAGAGGTATGATGTTCTATCTTTCAGACCTTTTTCTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTg                             <  1:315858/136‑1 (MQ=255)
                                                                                                         cAAAAACTGCTTTCTGAAGTGAACAAAAAAGCGCAAACTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTg                             <  1:1249399/143‑1 (MQ=255)
                                                                                                          aaaaaGTGCTTTCTGAACTGAACAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTgg                            <  1:1185903/143‑1 (MQ=255)
                                                                                                          aaaaaGTGCTTTCTGAACTGAACAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTgg                            <  2:1313045/143‑1 (MQ=255)
                                                                                                          aaaaaGTGCTTTCTGAACTGAACAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTgg                            >  2:1004179/1‑143 (MQ=255)
                                                                                                           aagaGTGCTTTCTGAACGGAACAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTggg                           <  1:994752/140‑1 (MQ=255)
                                                                                                            aaaGTGCTTTCTGAACTGAACAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTgggg                          >  2:1178103/1‑143 (MQ=255)
                                                                                                             aaGTGCTTTCTGAACTGAACAAAAAATAGCAAAGTTAGGCGCGTAGGGTACAGAGGTAAGATGTTCTATCTGTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCATTACTGACTGGGGc                         <  1:941392/143‑1 (MQ=255)
                                                                                                                    ttCTGAACTGAACAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGtttttt                  >  1:28832/1‑143 (MQ=255)
                                                                                                                         aactgaacAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTgg             <  1:149349/143‑1 (MQ=255)
                                                                                                                         aactgaacAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTCTTTATTgg             >  1:1044535/1‑143 (MQ=255)
                                                                                                                         aactgaacAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGATTTTTAGTgg             >  2:39741/1‑143 (MQ=255)
                                                                                                                         aaatgaaaAAAAAAGAGTAAAGGTAGTCGCGTAGGGTACAGAGGTAAAATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTgg             <  1:410957/140‑1 (MQ=255)
                                                                                                                           cggaacAAAAAAGAGATAAGTTAGTCGCGTAGGGGACAGAGGTAAGATGTTCTATCTTGCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGc           <  2:1201980/141‑1 (MQ=255)
                                                                                                                              aacAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAACCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGGAATGACCGGGGCGTTTTTTATTTGGCgaa        >  1:917571/1‑143 (MQ=255)
                                                                                                                               acAAAAAAGAGCAAAGTTCGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTGCACGTAATCGGATGTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCgaaa       <  1:1304266/143‑1 (MQ=255)
                                                                                                                                cAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTCTATTGGGCgaaac      >  1:736368/1‑142 (MQ=255)
                                                                                                                                    aaaGAGCAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGACCGAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGaaagaaaa  >  2:265831/1‑143 (MQ=255)
                                                                                                                                          |                                                                                                                                        
CTAAGATTCGCCTGAACATGCCTTTTTTCGTAAGTAAGCAACATAAGCTGTCACGTTTTGTGATGGCTATTAGAAATTCCTATGCAACAACTG‑AAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAA  >  NC_000913/4049748‑4050021

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