Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,049,886 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,8860TC97.3% 225.5 / ‑4.3 76intergenic (+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/39);  total (34/42)
Fisher's exact test for biased strand distribution p-value = 4.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

CTAAGATTCGCCTGAACATGCCTTTTTTCGTAAGTAAGCAACATAAGCTGTCACGTTTTGTGATGGCTATTAGAAATTCCTATGCAACAACTGAAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAA  >  NC_000913/4049749‑4050022
                                                                                                                                          |                                                                                                                                        
cTAAGATTCGCCTGAACATTCCTTTTTTCGTAAGTAAGCAACATAAGCTGTCACGTTTTGTGATGGCTATTAGAAATTCCTATGCAACAACTGAAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGCAAAGt                                                                                                                                     <  2:6277/143‑1 (MQ=255)
 tAAGATTCGCCTGAACATGCCTTTTTTCGTAAGTAAGCAACATAAGCTGTCACGTTTTGTGATGGCTATTAGAAATTCCTATGCAACAACTGAAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGCAAAGtt                                                                                                                                    >  1:550127/1‑143 (MQ=255)
   aGATTCGCCTGAACATGCCTTTTTTCGTAAGTAAGCAACATAAGCTGTCACGTTTTGTGATGGCTATTAGAAATTCCTATGCAACAACTGAAAAAAAATTACACAAAGTGCTTTCTGAACTGAACAAAAAAGAGAAAAGTCAg                                                                                                                                  >  2:358632/1‑143 (MQ=255)
         gCCTGAACATGCCTTTTTTCGTAAGTAAGCAACATAAGCTGTCACGTTTTGTGATGGCTATTAGAAATTCCTATGCAACAACTGAAAAAAAAATTACACAAAGTGCTTTCTTAACTGAACAAAACAGAGCAAACTTGGTcgcg                                                                                                                             >  2:18323/1‑143 (MQ=255)
         gCCTGAACATGCCTTTTTTCGTAAGTAAGCAACATAAGCTGTCACGTTTTGTGATGGCTATTAGAAATTCCTATGCAACAACTGAAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGCAAAGTTAGTCGCGt                                                                                                                            >  1:217780/1‑143 (MQ=255)
                   gCCTTTTTTCGTAAGTAAGCAACATAAGCTGTCACGTTTTGTGATGGCTATTAGAAATTCCTATGCAACAACTGAAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACaga                                                                                                                  >  2:1309724/1‑143 (MQ=255)
                   gCCTTTTTTCGTAAGTAAGCAACATAAGCTGTCACGTTTTGTGATGGCTATTAGAAATTCCTATGCAACAACTGAAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACaga                                                                                                                  >  2:1309944/1‑143 (MQ=255)
                   gCCTTTTTTCGTAAGTAAGCAACATAAGCTGTCACGTTTTGTGATGGCTATTAGAAATTCCTATGCAACAACTGAAAAAAAATTACAAAAAGTGCTGTCTGAACTGAACAAAAAAGAGCAAAGTTAGTAGCGCAGGGTACaga                                                                                                                  >  1:821609/1‑143 (MQ=255)
                   gCCTTCTCTCGTAAGTAAGCAACATAACCTGTCACGTTTTGTGATGGCTATTAGAAATTCCTATGCAACAACTGAAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACaga                                                                                                                  <  1:594466/143‑1 (MQ=255)
                    ccTTTTTTCGTAAGTAAGCAACAAAAGCAGTCACATACTATGATCGCTATTAGAAATTCCTATGCAACAACTGAAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACagag                                                                                                                 <  1:196015/143‑1 (MQ=255)
                      ttttttCGTAAGTAAGCAACATAAGCTGTCACGTTTTGTGATGGCTATTAGAAATTCCTATGCAACAACTGAAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACAGAGGt                                                                                                               >  1:1287247/1‑143 (MQ=255)
                         tttCGTAAGTAAGCAACTTAAGCTGTCACGGTTTGTGATGGCTATTAGAAATTCCTATGCAACAACTGAAAAAAAATTACAAAAAGTGCTTTCTTAACTGAACAAAAAAGAGCAAAGTTAGTCTCGTAGGGTACAGAGGTAAg                                                                                                            <  2:353856/143‑1 (MQ=255)
                           tCGTAAGTAAGCAACATAAGCTGTCACGTATTGTGATGGCTATTAGAAATTCCAATGCAACAACTGAAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGCAAAGGTAGACGCGTGGGGTGCAGAGGTAAGAt                                                                                                          >  2:462770/1‑143 (MQ=255)
                             gtaagtaaGCAACATAAGCTGTCACGTTTTGTGATGGCTATTAGAAATTCCTATGCAACAACTGAAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGCAAAGTTCGTCGCGTAGGGTACAGAGGTAAGATGt                                                                                                        >  1:949665/1‑143 (MQ=255)
                             gtaagtaaGCAACATAAGCTGGCACGTTTTGTGATGGCTATTAGAAATTCCTATGCAACAACTGAAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGt                                                                                                        <  1:970645/143‑1 (MQ=255)
                              gaagtaagGAACATAAGCTGTGACGGTTTGTGATGGCTATTAGAAATTCCTATGCAACAACTGAAAAAAAAGTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGtt                                                                                                       <  2:916495/142‑1 (MQ=255)
                               aagtaagCAACATAAGCTGTCACGTTTTGTGATGGCTATTAGAAATTCCTATGCAACAACTGAAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTc                                                                                                      <  1:713840/143‑1 (MQ=255)
                               aagtaagCAACATAAGCTGTCACGTTTTGTGATGGCTATTAGAAATTCCTATGCAACAACTGAAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGCAAACGTAGTCGCGTCAGCCACCGAGGAAAGATGTTc                                                                                                      >  2:157323/1‑143 (MQ=255)
                                agtaagCAACATAAGCTGTCACGTTTTGTGATGGCTATTAGAAATTCCTATGCAACAACTGAAAAAAAATTACAAAATGTGCTTTCTGAACTGAACAAAAAAGGGCAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCt                                                                                                     <  2:345760/143‑1 (MQ=255)
                                     gCAACATAAGCTGTCACGTTTTGTGATGGCTATTAGAAATTCCTATGCAACAACTGAAAAGAAGTTACAACAAGTGCTTTCTGAACTGAACACCAAAGCGCACAGTTGGTCGCGTAGGGTACAGCGGTGTGATGTTCtcgctt                                                                                                >  1:189462/1‑138 (MQ=255)
                                     gCAACATAAGCTGTCACGTTTTGTGATGGCTATTAGAAATTCCTATGCAACAACTGAAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTCCtctcct                                                                                                >  1:184726/1‑141 (MQ=255)
                                           tAAGCTGCACCGTTTTGTGATGGCTATTAGAAATTCCTATGCAACAACTGAAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGAc                                                                                          <  2:776232/143‑1 (MQ=255)
                                            aaGCTGTCACGTTTTGTGATGGCTATTAGAAATTCCTATGCAACAACTGAAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGAcc                                                                                         >  1:829433/1‑143 (MQ=255)
                                               cTGTCACGTTTTGTGATGGCTATTAGAAATTCCTATTCAATAACTGAAAAAAAAGTACAAATAGTGCTTTCTGAACTGAACAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCttt                                                                                      <  2:1085229/143‑1 (MQ=255)
                                               cTGTCACGTTTTGTGATGGCTATTAGAAATTCCTATGCAACAACTGAAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACCGAGGTACGATGTTCTATCTTTCAGACCttt                                                                                      >  1:1132815/1‑143 (MQ=255)
                                               cTGTCACGTTTTGTGATGGCTATTAGAAATTCCTATGCAACAACTGAAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCttt                                                                                      <  1:401259/143‑1 (MQ=255)
                                               cTGTCACGTTTTGTGATGGCTATTAGAAATTCCTATGCAACAACTGAAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCttt                                                                                      <  2:993568/143‑1 (MQ=255)
                                                 gTCACGTTTTGTGATGGCTATTAGAAATTCCTATGCAACAACTGAAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTa                                                                                    <  2:350177/143‑1 (MQ=255)
                                                 gTCACGTTATGAAGTGGCTATTAGAAATTCCTATGCAAAAACTGAAAAAAAATGACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTa                                                                                    <  2:974250/143‑1 (MQ=255)
                                                  tCACGTTTTGTGATGGCTATTAGAAATTCCTATGCAACAACTGAAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTGAc                                                                                   >  2:1313516/1‑143 (MQ=255)
                                                   cACGTATTGTGATGGCTATGAGAAATTCCGAGGCAACAACTGAAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGCCAAGGTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACt                                                                                  <  1:1354330/143‑1 (MQ=255)
                                                     cGTTTTGTGATGGCTATTAGAAATTCCTATGCAACAACTGAAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTc                                                                                <  2:311825/143‑1 (MQ=255)
                                                       ttttGTGATGGCTATTAGAAATTCCTATGCAACAACTGAAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCAc                                                                              <  2:217780/143‑1 (MQ=255)
                                                        tttGTGATGGCTATTAGAAATTCCTATGCAACAACTGAAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACg                                                                             >  2:586115/1‑143 (MQ=255)
                                                        tttGTGATGGCTATTAGAAATTCCTATGCAACAACTGAAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACg                                                                             <  2:193396/143‑1 (MQ=255)
                                                         ttGTGATGGCTATTAGACATTCCTATGAAACAACAGAAAAAAAATTACAAAAAGTGCCTTCTGAACTGAACAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGt                                                                            <  2:1296583/143‑1 (MQ=255)
                                                           gtgATGGCTATTAGAAATTCCTATGCAACAACTGAAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGCAAAGTTAGTCGCGTAGGGAAAAGCGGTAAGACGTTCTATCTTTCAGACCTTTTACTTCACGTaa                                                                          >  2:235316/1‑143 (MQ=255)
                                                            tgATGGCTATTAGAAATTCCTATGCAACAACTGAAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAAt                                                                         <  2:829433/143‑1 (MQ=255)
                                                                 gCTATTAGAAATTCCTATGCAAAAAATGAAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGAt                                                                    <  2:118411/143‑1 (MQ=255)
                                                                   tATTAGAAATTCCTATGCAACAACTGAAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGAttt                                                                  <  2:231224/143‑1 (MQ=255)
                                                                     ttAGAAATTCCTATGCAACAACTGAAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTgg                                                                <  1:1329584/143‑1 (MQ=255)
                                                                              ccTATGCAACAACTGAAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATAtt                                                       <  1:1309944/143‑1 (MQ=255)
                                                                              ccTATGCAACAACTGAAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATAtt                                                       <  1:1309724/143‑1 (MQ=255)
                                                                               cgaTGCAACAACTGAAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATAttt                                                      <  2:1297857/141‑1 (MQ=255)
                                                                                   gCAACAACTGAAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGCAAAGTCAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCCGACCTTTTACTTCACGCAATCGGATTTCGCTCAATATTTTAgc                                                  >  2:846355/1‑143 (MQ=255)
                                                                                    caacgaCTGAAAAAAAATTACAAAAAGTGCTTGCTGAACTGAACAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAgcc                                                 <  2:322017/143‑1 (MQ=255)
                                                                                    caacaaCTGAAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCCCGTAATCGGATTTGGCTGCATATTTTAgcc                                                 >  1:1307942/1‑143 (MQ=255)
                                                                                     aacaacTGAAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAgccg                                                >  1:932385/1‑143 (MQ=255)
                                                                                        aacTGAAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGccc                                             >  2:593557/1‑143 (MQ=255)
                                                                                        aacTGAAAAAAAAGTACAAGAAGTGCTTTCTGAACTGAACAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGccc                                             <  2:1202584/143‑1 (MQ=255)
                                                                                           tGAAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATATTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATACTTTAGCCGCCCcag                                          >  2:750472/1‑143 (MQ=255)
                                                                                               aaaaaaaTTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCcagtca                                       <  1:715366/143‑1 (MQ=255)
                                                                                               aaaaaaaTTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCcagtca                                       >  1:825045/1‑143 (MQ=255)
                                                                                                aaaaaaTTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATTTTAAACACCCCCCACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCcagtcag                                      >  1:329596/1‑143 (MQ=255)
                                                                                                aaaaaaTTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCcagtcag                                      >  1:340133/1‑143 (MQ=255)
                                                                                                  aaaaTTACGAAAAGTTCTGTCCGAACTGAACAAAAAAGAGCAAAGTTAGACGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTCATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTa                                    <  2:84034/143‑1 (MQ=255)
                                                                                                   aaaTTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTaa                                   <  1:1260918/143‑1 (MQ=255)
                                                                                                      ttACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGa                                <  1:352675/143‑1 (MQ=255)
                                                                                                         ccataattGCGTTCTGAACTGAACAAAACACAGCAGAGTTAGTCGCGCAGGGTACAGAGGTATGATGTTCTATCTTTCAGACCTTTTTCTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTg                             <  2:315858/136‑1 (MQ=255)
                                                                                                         cAAAAACTGCTTTCTGAAGTGAACAAAAAAGCGCAAACTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTg                             <  2:1249399/143‑1 (MQ=255)
                                                                                                          aaaaaGTGCTTTCTGAACTGAACAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTgg                            <  1:1313045/143‑1 (MQ=255)
                                                                                                          aaaaaGTGCTTTCTGAACTGAACAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTgg                            <  2:1185903/143‑1 (MQ=255)
                                                                                                          aaaaaGTGCTTTCTGAACTGAACAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTgg                            >  1:1004179/1‑143 (MQ=255)
                                                                                                           aagaGTGCTTTCTGAACGGAACAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTggg                           <  2:994752/140‑1 (MQ=255)
                                                                                                            aaaGTGCTTTCTGAACTGAACAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTgggg                          >  1:1178103/1‑143 (MQ=255)
                                                                                                             aaGTGCTTTCTGAACTGAACAAAAAATAGCAAAGTTAGGCGCGTAGGGTACAGAGGTAAGATGTTCTATCTGTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCATTACTGACTGGGGc                         <  2:941392/143‑1 (MQ=255)
                                                                                                                    ttCTGAACTGAACAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGtttttt                  >  2:28832/1‑143 (MQ=255)
                                                                                                                         aactgaacAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTgg             <  2:149349/143‑1 (MQ=255)
                                                                                                                         aactgaacAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTCTTTATTgg             >  2:1044535/1‑143 (MQ=255)
                                                                                                                         aactgaacAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGATTTTTAGTgg             >  1:39741/1‑143 (MQ=255)
                                                                                                                         aaatgaaaAAAAAAGAGTAAAGGTAGTCGCGTAGGGTACAGAGGTAAAATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTgg             <  2:410957/140‑1 (MQ=255)
                                                                                                                           cggaacAAAAAAGAGATAAGTTAGTCGCGTAGGGGACAGAGGTAAGATGTTCTATCTTGCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGc           <  1:1201980/141‑1 (MQ=255)
                                                                                                                              aacAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAACCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGGAATGACCGGGGCGTTTTTTATTTGGCgaa        >  2:917571/1‑143 (MQ=255)
                                                                                                                               acAAAAAAGAGCAAAGTTCGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTGCACGTAATCGGATGTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCgaaa       <  2:1304266/143‑1 (MQ=255)
                                                                                                                                cAAAAAAGAGCAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTCTATTGGGCgaaac      >  2:736368/1‑142 (MQ=255)
                                                                                                                                    aaaGAGCAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGACCGAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGaaagaaaa  >  1:265831/1‑143 (MQ=255)
                                                                                                                                          |                                                                                                                                        
CTAAGATTCGCCTGAACATGCCTTTTTTCGTAAGTAAGCAACATAAGCTGTCACGTTTTGTGATGGCTATTAGAAATTCCTATGCAACAACTGAAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAA  >  NC_000913/4049749‑4050022

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: