Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,733,449 C→A 11.7% L1455M (CTG→ATG)  lhr → putative ATP‑dependent helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,733,4490CA11.7% 2224.4 / 193.6 673L1455M (CTG→ATG) lhrputative ATP‑dependent helicase
Reads supporting (aligned to +/- strand):  ref base C (318/275);  new base A (41/39);  total (359/314)
Fisher's exact test for biased strand distribution p-value = 7.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

Only 100 of 694 total aligned reads displayed.

ATTGCGCGATCTGGCGACACAAGCCACGCAAACGCGCCACTATACACCAGTGGCGCTCTCTGCCAACGATCCGGCTAATGTGTGGGGAAATCTTCTGCCCTGGCCTGCACATCCGGCAACGCTGGTTCCAACGCGTCGGGCGGGTGCGCTGGTGGTCGTTTCTGGCGGCAAATTGTTACTCTATCTGGCGCAAGGTGGCAAAAAAATGCTGGTCTGGCAGGAAAAAGAGGAATTACTCGCCCCAGAGGTTTTCCACGCGCTGACTACCGCACTGCGTCGCGAACCACGGCTGCGCTT  >  NC_000913/1733301‑1733597
                                                                                                                                                    |                                                                                                                                                    
 ttGCGCGATCTGGCGACACAAGCCACGCAAACGCGCCACTATACACCAGTGGCGCTCTCTGCCAACGATCCGGCTACTGTGTGGGGAAATCTTCTGCCCTGGCCTGCACATCCGGCAACGCTGGTTCCAACGCGTCGGGCGGGTGCGCt                                                                                                                                                     >  2:1999285/1‑149 (MQ=255)
     gcgATCTGGCGACACAAGCCACGCAAACGCGCCACTATACACCAGTGGCGCTCTCTGCCAACGATCCGGCTAATGTGTGGGGAAATCTTCTGCCCTGGCCTGCACATCCGGCAACGCTGGTTCCAATGCGTCGGGCGGGTGCGCtggtg                                                                                                                                                 >  4:5868892/1‑149 (MQ=255)
     gcgATCTGGCGACACAAGCCACGCAAACGCGCCACTATACACCAGTGGCGCTCTCTGCCAACGATCCGGCTAATGTGTGGGGAAATCTTCTGCCCTGGCCTGCACATCCGGCAACGCTGGTTCCAACGCGTCGGGCGGGTGCGCtggtg                                                                                                                                                 <  1:4966101/149‑1 (MQ=255)
     gcgATCTGGCGACACAAGCCACGCAAACGCGCCACTATACACCAGTGGCGCTCTCTGCCAACGATCCGGCTAATGTGTGGGGAAATCTTCTGCCCTGGCCTGCACATCCGGCAACGCTGGTTCCAACGCGTCGGGCGGGTGCGCtggtg                                                                                                                                                 >  1:2686270/1‑149 (MQ=255)
       gATCTGGCGACACAAGCCACGCAAACGCGCCACTATACACCAGTGGCGCTCTCTGCCAACGATCCGGCTAATGTGTGGGGAAATCTTCTGCCCTGGCCTGCACATCCGGCAACGCTGGTTCCAACGCGTCGGGCGGGTGCGCtggt                                                                                                                                                  <  2:4960569/146‑1 (MQ=255)
       gATCTGGCGACACAAGCCACGCAAACGCGCCACTATACACCAGTGGCGCTCTCTGCCAACGATCCGGCTAATGTGTGGGGAAATCTTCTGCCCTGGCCTGCACATCCGGCAACGCTGGTTCCAACGCGTCGGGCGGGTGCGCTggtggt                                                                                                                                               >  3:6242658/1‑149 (MQ=255)
          cTGGCGACACAAGCCACGCAAACGCGCCACTATACACCAGTGGCGCTCTCTGCCAACGATCCGGCTAATGTGTGGGGAAATCTTCTGCCCTGGCCTGCACATCCGGCAACGCTGGTTCCAACGCGTCGGGCGGGTGCGCTGGTGGTCGt                                                                                                                                            <  3:3784206/149‑1 (MQ=255)
          cTGGCGACACAAGCCACGCAAACGCGCCACTATACACCAGTGGCGCTCTCTGCCAACGATCCGGCTAATGTGTGGGGAAATCTTCTGCCCTGGCCTGCACATCCGGCAACGCTGGTTCCAACGCGTCGGGCGGGTGCGCTGGTGGTCGt                                                                                                                                            <  3:110214/149‑1 (MQ=255)
           tGGCGACACAAGCCACGCAAACGCGCCACTATACACCAGTGGCGCTCTCTGCCAACGATCCGGCTAATGTGTGGGGAAATCTTCTGCCCTGGCCTGCACATCCGGCAACGCTGGTTCCAACGCGTCGGGCGGGTGCGCTGGTGGTCGtt                                                                                                                                           >  4:2547470/1‑149 (MQ=255)
              agACCCAAGCCACGCAAACGCGCCACTAAATACGAGAGGCGCTGTCTGCCAACGATCCGGCTGAAGTGGGGTAAAATCTTGTGCGCTGGCCTGCCCGTCCGGCACCGCTGTTTCCAACGCGGCGGGAGGGTGCGATTGTGGTCGTATCt                                                                                                                                        <  4:1867858/148‑1 (MQ=255)
               gACACAAGCCACGCAAACGCGCCACTATACACCAGTGGCGCTCTCTGCCAACGATCCGGCTAATGTGTGGGGAAATCTTCTGCCCTGGCCTGCACATCCGGCAACGCTGGTTGCAACGCGTCGGGCGGGTGCGCTGGTGGTCGTTTCTg                                                                                                                                       >  3:4764041/1‑149 (MQ=255)
               gACACAAGCCACGCAAACGCGCCACTATACACCAGTGGCGCTCTCTGCCAACGATCCGGCTAATGTGTGGGGAAATCTTCTGCCCTGGCCTGCACATCCGGCAACGCTGGTTCCAACGCGTCGGGCGGGTGCGCTGGTGGTCGTTTCTg                                                                                                                                       >  3:4250699/1‑149 (MQ=255)
                acacAAGCCACGCAAACGCGCCACTATACTCCAGTGCCGCTCTCTGCCAACGATCCGGCTAATCTGTGGTGAAATCTTCTGCCCTGGCCTGCACATCCGGCAACGCTGGTTCCAACGCGGCGGGGGGCTGCGATGGTGGTCGTTTCggg                                                                                                                                      <  4:4428842/149‑1 (MQ=255)
                acacAAGCCACGCAAACGCGCCACTATACACCAGTGGCGCTCTCTGCCAACGATCCGGCTAATGTGTGGGGAAATCTTCTGCCCTGGCCTGCACATCCGGCAACGCTGGTTCCAACGCGTCGGGCGGGTGCGCTGGTGGTCGTTTCTgg                                                                                                                                      <  2:7184742/149‑1 (MQ=255)
                 cacaAGCCACGCAAACGCGCCACTATACACCTGTGGCGCTCTCTGCCACCGATCCGGCTAATGTGTGGGGAGATCTTCTGCCCTGGCCTGCACATCCGGCAACGCTGGTTCCAACGCGGCGGGCGGGTGCGCTGGTGGTCGTTTCTgga                                                                                                                                     <  2:3281431/149‑2 (MQ=255)
                 cacaAGCCACGCAAACGCGCCACTATACACCAGTGGCGCTCTCTGCCAACGATCCGGCTAATGTGTGGGGAAATCTTCTGCCCTGGCCTGCACATCCGGCAACGCTGGTTCCAACGCGTCGGGCGGGTGCGCTGGTGGTCGTTTCTggc                                                                                                                                     >  2:2302566/1‑149 (MQ=255)
                 cacaAGCCACGCAAACGCGCCACTATACACCAGTGGCGCTCTCTGCCAACGATCCGGCTAATGTGTGGGGAAATCTTCTGCCCTGGCCTGCACATCCGGCAACGCTGGTTCCAACGCGTCGGGCGGGTGCGCTGGTGGTCGTTTCTggc                                                                                                                                     >  1:6873488/1‑149 (MQ=255)
                   caAGCCACGCAAACGCGCCACTATACACCAGTGGCGCTCTCTGCCAACGATCCGGCTAATGTGTGGGGAAATCTTCTGCCCTGGCCTGCACATCCGGCAACGCTGGTTCCAACGCGTCGGGCGGGTGCGCTGGTGGTCGTTTCTggcgg                                                                                                                                   <  1:1416262/149‑1 (MQ=255)
                    aaGCCACGCAAACGCGCCACTATACACCAGTGGCGCTCTCTGCCAACGATCCGGCTAATGTGTGGGGAAATCTTCTGCCCTGGCCTGCACATCCGGCAACGCTGGTTCCAACGCGTCGGGCGGGTGCGCTGGTGGTCGTTTCTggcggc                                                                                                                                  >  2:3067206/1‑149 (MQ=255)
                       ccACGCAAACGCGCCACTATACACCAGTGGCGCTCTCTGCCAACGATCCGGCTAATGTGTGGGGAAATCTTCTGCCCTGGCCTGCACATCCGGCAACGCTGGTTCCAACGCGTCGGGCGGGTGCGCTGGTGGTCGTTTCTGGCGGCaaa                                                                                                                               >  2:2749917/1‑149 (MQ=255)
                           gCAAACGCGCCACTATACACCAGTGGCGCTCTCTGCCAACGATCCGGCTAATGTGTGGGGAAATCTTCTGCCCTGGCCTGCACATCCGGCAACGCTGGTTCCAACGCGTCGGGCGGGTGCGCTGGTGGTCGTTTCTGGCGGCAAATTGt                                                                                                                           <  3:7097673/149‑1 (MQ=255)
                            cAAACGCGCCACTATACACCAGTGGCGCTCTCTGCCAACGATCCGGCTAATGTGTGGGGAAATCTTCTGCCCTGGCCTGCACATCCGGCAACGCTGGTTCCAACGCGTCGGGCGGGTGCGCTGGTGGTCGTTTCTGGCGGCAAATTGtt                                                                                                                          <  1:718824/149‑1 (MQ=255)
                             aaaCGCGCCACTATACACCAGTGGCGCTCTCTGCCAACGATCCGGCTAATGTGTGGGGAAATCTTCTGCCCTGGCCTGCACATCCGGCAACGCTGGTTCCAACGCGTCGGGCGGGTGCGCtg                                                                                                                                                    >  2:1322035/1‑122 (MQ=255)
                             aaaCGCGCCACTATACACCAGTGGCGCTCTCTGCCAACGATCCGGCTAATGTGTGGGGAAATCTTCTGCCCTGGCCTGCACATCCGGCAACGCTGGTTCCAACGCGTCGGGCGGGTGCGCTGGTGGTCGTTTCTgg                                                                                                                                      <  3:6606023/136‑1 (MQ=255)
                             aaaCGCGCCACTATACACCAGTGGCGCTCTCTGCCAACGATCCGGCTAATGTGTGGGGAAATCTTCTGCCCTGGCCTGCACATCCGGCAACGCTGGTTCCAACGCGTCGGGCGGGTGCGCTGGTGGTCGTTTCTGGCGGCAAATTGTTa                                                                                                                         <  2:4407654/149‑1 (MQ=255)
                             aaaCGCGCCACTATACACCAGTGGCGCTCTCTGCCAACGATCCGGCTAATGTGTGGGGAAATCTTCTGCCCTGGCCTGCACATCCGGCAACGCTGGTTCCAACGCGTCGGGCGGGTGCGCTGGTGGTCGTTTCTGGCGGCAAATTGTTa                                                                                                                         <  4:819193/149‑1 (MQ=255)
                              aaCGCGCCACTGTACACCAGTGGCGCTCTCTGCCAACGATCCGGCTAATGTGTGGGGAAATCTTCTGCCCTGGCCTGCACATCCGGCAACGCTGGTTCCAACGCGTCGGGCGGGTGCGCTGGTGGTCGTTTCTGGCGGCAAATTGTTAc                                                                                                                        >  3:4041405/1‑149 (MQ=255)
                              aaCGCGCCACTATACACCAGTGGCGCTCTCTGCCAACGATCCGGCTAATGTGTGGGGAAATCTTCTGCCCTGGCCTGCACATCCGGCAACGCTGGTTCCAACGCGTCGGGCGGGTGCGCTGGTGGTCGTTTCTGGCGGCAAATTGTTAc                                                                                                                        >  2:1882272/1‑149 (MQ=255)
                                   gcCACTATACACCAGTGGCGCTCTCTGCCAACGATCCGGCTAATGTGTGGGGAAATCTTCTGCCCTGGCCTGCACATCCGGCAACGCTGGTTCCAACGCGTCGGGCGGGTGCGCTGGTGGTCGTTTCTGGCGGCACATTGTTACTCTAt                                                                                                                   >  4:1463910/1‑149 (MQ=255)
                                   gcCACTATACACCAGTGGCGCTCTCTGCCAACGATCCGGCTAATGTGTGGGGAAATCTTCTGCCCTGGCCTGCACATCCGGCAACGCTGGTTCCAACGCGTCGGGCGGGTGCGATGGTGGTCGTTTCTGGCGGCAAATTGTTACTCTAt                                                                                                                   >  1:594497/1‑149 (MQ=255)
                                     cACTATACACCAGTGGCGCTCTCTGCCAACGATCCGGCTAATGTGTGGGGAAATCTTCTGCCCTGGCCTGCACATCCGGCAACGCTGGTTCCAACGCGTCGGGCGGGTGCGCTGGTGGTCGTTTCTGGCGGCAAATTGTTACTCTATCt                                                                                                                 >  4:3467762/1‑149 (MQ=255)
                                            cacCAGTGGCGCTCTCTGCCAACGATCCGGCTAATGTGTGGGGAAATCTTCTGCCCTGGCCTGCACATCCGGCAACGCTGGTTCCAACGCGTCGGGCGGGTGCGCTGGTGGTCGTTTCTGGCGGCAAATTGt                                                                                                                           >  2:6949850/1‑132 (MQ=255)
                                            cacCAGTGGCGCTCTCTGCCAACGATCCGGCTAATGTGTGGGGAAATCTTCTGCCCTGGCCTGCACATCCGGCAACGCTGGTTCCAACGCGTCGGGCGGGTGCGCTGGTGGTCGTTTCTGGCGGCAAATTGTTACTCTATCTGGCGCaa                                                                                                          <  1:5725275/149‑1 (MQ=255)
                                             acCAGTGGCGCTCTCTGCCAACGATCCGGCTAATGTGTGGGGAAATCTTCTGCCCTGGCCTGCACATCCGGCAACGCTGGTTCCAACGCGTCGGGCGGGTGCGCTGGTGGTCGTTTCTGGCGGCAAATTGTTACTCTATCTGGCGCAAg                                                                                                         >  4:6991820/1‑149 (MQ=255)
                                             acCAGTGGCGCTCTCTGCCAACGATCCGGCTAATGTGTGGGGAAATCTTCTGCCCTGGCCTGCACATCCGGCAACGCTGGTTCCAACGCGTCGGGCGGGTGCGCTGGTGGTCGTTTCTGGCGGCAAATTGTTACTCTATCTGGCGCAAg                                                                                                         <  3:339621/149‑1 (MQ=255)
                                             acCAGTGGCGCTCTCTGCCAACGATCCGGCTAATGTGTGGGGAAATCTTCTGCCCTGGCCTGCACATCCGGCAACGCTGGTTCCAACGCGTCGGGCGGGTGCGCTGGTGGTCGTTTCTGGCGGCAAATTGTTACTCTATCTGGCGCAAg                                                                                                         >  4:3209201/1‑149 (MQ=255)
                                             acCAGTAGCGCTCTCTGCCAACGATCCGGCTAATGTGTGGGGAAATCTTCTGCCCTGGCCTGCACATCCGGCAACGCTGGTTCCAACGCGTCGGGCGGGTGCGCTGGTGGTCGTTTCTGGCGGCAAATTGTTACTCTATCTGGCGCAAg                                                                                                         >  2:5155491/1‑149 (MQ=255)
                                               cAGTGGCGCTCTCTGCCAACGATCCGGCTAATGTGTGGGGAAATCTTCTGCCCTGGCCTGCACATCCGGCAACGCTGGTTCCAACGCGTCGGGCGGGTGCGCTGGTGGTCGTTTCTGGCGGCAAATTGTTACTCTATCTGGCGCAAGGt                                                                                                       <  3:2218936/149‑1 (MQ=255)
                                                aGTGGCGCTCTCTGCCAACGATCCGGCTAATGTGTGGGGAAATCTTCTGCCCTGGCCTGCACATCCGGCAACGCTGGTTCCAACGCGTCGGGCGGGTGCGCTGGTGGTCGTTTCTGGCGGCAAATTGTTACTCTATCTGGCGCAAGGTg                                                                                                      <  1:5206365/149‑1 (MQ=255)
                                                  tGGCGCTCTCTGCCAACGATCCGGCTAATGTGTGGGGAAATCTTCTGCCCTGGCCTGCACATCCGGCAACGCTGGTTCCAACGCGTCGGGCGGGTGCGCTGGTGGTCGTTTCTGGCGGCAAATTGTTACTCTATCTGGCGCAAGGTGGc                                                                                                    >  1:34142/1‑149 (MQ=255)
                                                      gcTCTCTGCCAACGATCCGGCTAATGTGTGGGGAAATCTTCTGCCCTGGCCTGCACATCCGGCAACGCTGGTTCCAACGCGTCGGGCGGGTGCGATGGTGGTCGTTTCTGGCGGCAAATTGTTACTCTATCTGGCGCAAGGTGGCaaaa                                                                                                >  4:6442316/1‑149 (MQ=255)
                                                            tGCCAACGATCCGGCTAATGTGTGGGGAAATCTTCTGCCCTGGCCTGCACATCCGGCAACGCTGGTTCCAACGCGTCGGGCGGGTGCGATGGTGGTCGTTTCTGGCGGCAAATTGTTACTCTATCTGGCGCAAGGTGGCAAAAAAATGc                                                                                          <  3:2712053/149‑1 (MQ=255)
                                                            tGCCAACGATCCGGCTAATGTGTGGGGAAATCTTCTGCCCTGGCCTGCACATCCGGCAACGCTGGTTCCAACGCGTCGGGCGGATGCGCTGGTGGTCGTTTCTGGCGGCAAATTGTTACTCTATCTGGCGCAAGGTGGCAAAAAAATGc                                                                                          >  4:4032379/1‑149 (MQ=255)
                                                              ccAACGATCCGGCTAATGTGTGGGGAAATCTTCTGCCCTGGCCTGCACATCCGGCAACGCTGGTTCCAACGCGTCGGGCGGGTGCGCtggtg                                                                                                                                                 <  3:1964339/92‑1 (MQ=255)
                                                              ccAACGATCCGGCTAATGTGTGGGGAAATCTTCTGCCCTGGCCTGCACATCCGGCAACGCTGGTTCCAACGCGTCGGGCGGGTGCGCTGGTGGTCGTTTCTGGCGGCAAATTGTTACTCTATCTGGCGCAAGGTGGCAAAAAAATGCTg                                                                                        >  4:6802235/1‑149 (MQ=255)
                                                                aaCGATCCGGCTAATGTGTGGGGAAATCTTCTGCCCTGGCCTGCACATCCGGCAACGCTGGTTCCAACGCGTCGGGCGGGTGCGCTGGTGGTCGTTTCTGGCGGCAAATTGTTACTCTATCTGGCGCAAGGTGGCAAAAAAATGCTGGt                                                                                      <  4:2219114/149‑1 (MQ=255)
                                                                aaCGATCCGGCTAATGTGTGGGGAAATCTTCTGCCCTGGCCTGCACATCCGGCAACGCTGGTTCCAACGCGTCGGGCGGGTGCGCTGGTGGTCGTTTCTGGCGGCAAATTGTTACTCTATCTGGCGCAAGGTGGCAAAAAAATGCTGGt                                                                                      >  2:1228920/1‑149 (MQ=255)
                                                                aaCGATCCGGCTAATGTGTGGGGAAATCTTCTGCCCTGGCCTGCACATCCGGCAACGCTGGTTCCAACGCGTCGGGCGGGTGCGCTGGTGGTCGTTTCTGGCGGCAAATTGTTACTCTATCTGGCGCAAGGTGGCAAAAAAATGCTGGt                                                                                      >  4:2437605/1‑149 (MQ=255)
                                                                 aCGATCCGGCTAATGTGTGGGGAAATCTTCTGCCCTGGCCTGCACATCCGGCAACGCTGGTTCCAACGCGTCGGGCGGGTGCGCTGGTGGTCGTTTCTGGCGGCAAATTGTTACTCTATCTGGCGCAAGGTGGCAAAAAAATGCTGGTc                                                                                     <  2:4126143/149‑1 (MQ=255)
                                                                   gATCCGGCTAATGTGTGGGGAAATCTTCTGCCCTGGCCTGCACATCCGGCAACGCTGGTTCCAACGCGTCGGGCGGGTGCGCTGGTGGTCGTTTCTGGCGGCAAATTGTTACTCTATCTGGCGCAAGGTGGCAAAAAAATGCTGGTCTg                                                                                   >  4:3777218/1‑149 (MQ=255)
                                                                     tCCGGCTAATGTGTGGGGAAATCTTCTGCCCTGGCCTGCACATCCGGCAACGCTGGTTCCAACGCGTCGGGCGGGTGCGCtggt                                                                                                                                                  >  4:1172406/1‑84 (MQ=255)
                                                                      ccGGCTAATGTGTGGGGAAATCTTCTGCCCTGGCCTGCACATCCGGCAACGCTGGTTCCAACGCGTCGGGCGGGTGCGCTGGTGGTCGTTTCTGGCGGCAAATTGTTACTCTATCTGGCGCAAGGTGGCAAAAAAATGCTGGGCTGGCa                                                                                >  1:100831/1‑149 (MQ=255)
                                                                       cGGCTAATGTGTGGGGAAATCTTCTGCCCTGGCCTGCACATCCGGCAACGCTGGTTCCAATGCGTCGGGCGGGTGCGCTGGTGGTCGTTTCTGGCGGCAAATTGTTACTCTATCTGGCGCAAGGTGGCAAAAAAATGCTGGTCTGGCAg                                                                               <  3:4937855/149‑1 (MQ=255)
                                                                           tAATGTGTGGGGAAATCTTCTGCCCTGGCCTGCACATCCGGCAACGCTGGTTCCAACGCGTCGGGCGGGTGCGCTGGTGGTc                                                                                                                                              <  1:3228027/82‑1 (MQ=255)
                                                                               gtgtgGGGAAATCTTCTGCCCTGGCCTGCACATCCGGCAACGCTGGTTCCAACGCGTCGGGCGGGTGCGCTGGTGGTCGTTTCTGGCGGCAAATTGTTACTCTATCTGGCGCAAGGTGGCAAAAAAATGCTGGTCTGGCAGGAAAAaga                                                                       <  3:453612/149‑1 (MQ=255)
                                                                                tgtgGGGAAATCTTCTGCCCTGGCCTGCACATCCGGCAACGCTGGTTCCAACGCGTCGGGCGGGTGCGATGGTGGTCGTTTCTGGCGGCAAATTGTTACTCTATCTGGCGCAAGGTGGCAAAAAAATGCTGGTCTGGCAGGAAAAagag                                                                      <  3:3062690/149‑1 (MQ=255)
                                                                                 gtgGGGAAATCTTCTGCCCTGGCCTGCACATCCGGCAACGCTGGTTCCAACGCGTCGGGCGGGTGCGCTGGTGGTCGTTTCTGGCGGCAAATTGTTACTCTATCTGGCGCAAGGTGGCAAAAAAATGCTGGTCTGGCAGGAAAAAGAgg                                                                     <  4:5057518/149‑1 (MQ=255)
                                                                                    gggAAATCTTCTGCCCTGGCCTGCACATCCGGCAACGCTGGTTCCAACGCGTCGGGCGGGTGCGCTGGTGGTCGTTTCTGGCGGCAAATTGTTACTCTATCTGGCGCAAGGTGGCAAAAAAATGCTGGTCTGGCAGGAAAAAGAGGAAt                                                                  >  1:1925944/1‑149 (MQ=255)
                                                                                      gAAATCTTCTGCCCTGGCCTGCACATCCGGCAACGCTGGTTCCAACGCGTCGGGCGGGTGCGCTGGTGGTCGTTTCTGGCGGCAAATTGTTACTCTATCTGGCGCAAGGTGGCAAAAAAATGCTGGTCTGGCAGGAAAAAGAGGAATTa                                                                >  1:6550937/1‑149 (MQ=255)
                                                                                         aTCTTCTGCCCTGGCCTGCACATCCGGCAACGCTGGTTCCAACGCGTCGGGCGGGTGCGCTGGTGGTCGTTTCTgg                                                                                                                                      >  1:2435022/1‑76 (MQ=255)
                                                                                         aTCTTCTGCCCTGGCCTGCACATCCGGCAACGCTGGTTCCAACGCGTCGGGCGGGTGCGCTGGTGGTCGTTTCTGGCGGCAAATTGTTACTCTATCTGGCGCAAGGTGGCAAAAAAATGCTGGTCTGGCAGGAAAAAGAGGAATTACTc                                                             >  2:5390805/1‑149 (MQ=255)
                                                                                          tcttctGCCCTGGCCTGCACATCCGGCAACGCTGGTTCCAACGCGTCGGGCGGGTGCGCTGGTGGTCGTTTCTGGCGGCAAATTGTTACTCTATCTGGCGCAAGGTGGCAAAAAAATGCTGGTCTGGCAGGAAAAAGAGGAATTACTCg                                                            >  4:6209041/1‑149 (MQ=255)
                                                                                            ttctGCCCTGGCCTGCACATCCGGCAACGCTGGTTCCAACGCGTCGGGCGGGTGCGCTGGTGGTCGTTTCTGGCGGCAAATTGTTACTCTATCTGGCGCAAGGTGGCAAAAAAATGCTGGTCTGGCAGGAAAAAGAGGAATTACTCGcc                                                          >  4:1727184/1‑149 (MQ=255)
                                                                                                 cccTGGCCTGCACATCCGGCAACGCTGGTTCCAACGCGTCGGGCGGGTGCGCTGGTGGTCGTTTCTGGCGGCAAATTGTTACTCTATCTGGCGCAAGGt                                                                                                       >  4:474799/1‑99 (MQ=255)
                                                                                                 cccTGGCCTGCACATCCGGCAACGCTGGTTCCAACGCGTCGGGCGGGTGCGCTGGTGGTCGTTTCTGGCGGCAAATTGTTACTCTATCTGGCGCAAGGTGGCAAAAAAATGCTGGTCTGGCAGGATAAAGAGGAATTACTCGCCCCaga                                                     <  3:962487/149‑1 (MQ=255)
                                                                                                    tGGCCTGCACATCCGGCAACGCTGGTTCCAACGCGTCGGGCGGGTGCGCTGGTGGTCGTTTCTGGCGGCAAATTGTTa                                                                                                                         >  3:1568344/1‑78 (MQ=255)
                                                                                                    tGGCCTGCACATCCGGCAACGCTGGTTCCAACGCGTCGGGCGGGTGCGCTGGTGGTCGTTTCTGGCGGCAAATTGTTACTCTATCTGGCGCAAGGTGGCAAAAAAATGCTGGTCTGGCAGGAAAAAGAGGAATTACTCGCCCCAGAGGt                                                  >  2:5700254/1‑149 (MQ=255)
                                                                                                      gCCTGCACATCCGGCAACGCTGGTTCCAACGCGTCGGGCGGGTGCGATGGTGGTCGTTTCTGGCGGCAAATTGTTACTCTATCTGGCGCAAGGTGGCAAAAAAATGCTGGTCTGGCAGGAAAAAGAGGAATTACTCGCCCCAGAGGttt                                                >  3:129621/1‑149 (MQ=255)
                                                                                                       ccTGCACATCCGGCAACGCTGGTTCCAACGCGTCGGGCGGGTGCGCTGGTGGTCGTTTCTGGCGGCAAATTGTTACTCTATCTGGCGCAAGGTGGCAAAAAAATGCTGGTCTGGCAGGAAAAAGAGGAATTACTCGCCCCAGAGGtttt                                               <  3:6802235/149‑1 (MQ=255)
                                                                                                       ccTGCACATCCGGCAACGCTGGTTCCAACGCGTCGGGCGGGTGCGCTGGTGGTCGTTTCTGGCGGCAAATTGTTACTCTATCTGGCGCAAGGTGGCAAAAAAATGCTGGTCTGGCAGGAAAAAGAGGAATTACTCGCCCCAGAGGtttt                                               <  1:124793/149‑1 (MQ=255)
                                                                                                       ccTGCACATCCGGCAACGCTGGTTCCAACGCGTCGGGCGGGTGCGATGGTGGTCGTTTCTGGCGGCAAATTGTTACTCTATCTGGCGCAAGGTGGCAAAAAAATGCTGGTCTGGCAGGAAAAAGAGGaa                                                                   <  1:1752024/129‑1 (MQ=255)
                                                                                                        cTGCACATCCGGCAACGCTGGTTCCAACGCTTCGGGCGGGTGCGCTGGTGGTCGTTTCTGGCGGCAAATTGTTACTCTATCTGGCGCAAGGTGGCAAAAAAATGCTGGTCTGGCAGGAAAAAGAGGAATTACTCGCCCCAGAGGTTTTc                                              <  2:6324436/149‑1 (MQ=255)
                                                                                                        cTGCACATCCGGCAACGCTGGTTCCAACGCGTCGGGCGGGTGCGCTGGTGGTCGTTTCTGGCGGCAAATTGTTACTCTATCTGGCGCAAGGTGGCAAAAAAATGCTGGTCTGGCAGGAAAAaga                                                                       <  1:5371833/124‑1 (MQ=255)
                                                                                                        cTGCACATCCGGCAACGCTGGTTCCAACGCGTCGGGCGGGTGCGCTGGTGGTCGTTTCTGGCGGCAAATTGTTACTCTATCTGGCGCAAGGTGGCAAAAAAATGCTGGTCTGGCAGGAAAAAGAGGAATTACTCGCCCCAGAGGTTTTc                                              >  1:4738768/1‑149 (MQ=255)
                                                                                                          gCACATCCGGCAACGCTGGTTCCAACGCGTCGGGCGGGTGCGCTGGTGGTCGTTTCCGGCGGCAAATTGTTACTCTATCTGGCGCAAGGTGGCAAAAAAATGCTGGTCTGGCAGGAAAAaga                                                                       >  3:5167003/1‑122 (MQ=255)
                                                                                                          gCACATCCGGCAACGCTGGTTCCAACGCGTCGGGCGGGTGCGATGGTGGTCGTTTCTGGCGGCAAATTGTTACTCTATCTGGCGCAAGGTGGCAAAAAAATGCTGGTCTGGCAGGAAAAAGAGGAATTACTCGCCCCAGAGGTTTTCCa                                            >  2:6673876/1‑149 (MQ=255)
                                                                                                           cacaTCCGGCAACGCTGGTTCCAACGCGTCGGGCGGGTGCGCTGGTGGTCGTTTCTGGCGGCAAATTGTTACTCTATCTGGCGCAAGGTGGCAAAAAAATGCTGGTCTGGCAGGAAAAAGAGGAATTACTCGCCCCAGAGGTTTTCCAc                                           <  1:626763/149‑1 (MQ=255)
                                                                                                            acaTCCGGCAACGCTGGTTCCAACGCGTCGGGCGGGTGCGATGGTGGTCGTTTCTGGCGGCAAATTGTTACTCTATCTGGCGCAAGGt                                                                                                       <  1:2322088/88‑1 (MQ=255)
                                                                                                             caTCCGGCAACGCTGGTTCCAACGCGTCGGGCGGGTGCGCTGGTGGTCGTTTCTGGCGGCAAATTGTTACTCTATCTGGCGCAAGGTGGCAAAAAAATGCTGGTCTGGCAGGAAAAAGAGGAATTACTCGCCCCAGAGGTTTTCCAcgc                                         >  1:1143761/1‑149 (MQ=255)
                                                                                                                ccGGCAACGCTGGTTCCAACGCGTCGGGCGGGTGCGCTGGTGGTCGTTTCTGGCGGCAAATTGTTACTCTATCTGGCGCAAGGTGGCAAAAAAATGCTGGTCTGGCAGGAAAAAGAGGAATTACTCGCCCCAGAGGTTTTCCACGCGCt                                      <  1:2159520/149‑1 (MQ=255)
                                                                                                                ccGGCAACGCTGGTTCCAACGCGTCGGGCGGGTGCGCTGGTGGTCGTTTCTGGCGGCAAATTGTTACTCTATCTGGCGCAAGGTGGCAAAAAAATGCTGGTCTGGCAGGAAAAAGAGGAATTACTCGCCCCAGAGGTTTTCCACGCGCt                                      >  1:4178921/1‑149 (MQ=255)
                                                                                                                  ggCAACGCTGGTTCCAACGCGTCGGGCGGGTGCGCTGGTGGTCGTTTCTGGCGGCAAATTGTTACTCTATCTGGCGCAAGGTGGCAAAAAAATGCTGGTCTGGCAGGAAAAAGAGGAATTACTCGCCCCAGAGGTTTTCCACGCGCTGa                                    <  3:5558180/149‑1 (MQ=255)
                                                                                                                  ggCAACGCTGGTTCCAACGCGTCGGGCGGGTGCGCTGGTGGTCGTTTCTGGCGGCAAATTGTTACTCTATCTGGCGCAAGGTGGCAAAAAAATGCTGGTCTGGCAGGAAAAAGAGGAATTACTCGCCCCAGAGGTTTTCCACGCGCTGa                                    >  4:5532004/1‑149 (MQ=255)
                                                                                                                  ggCAACGCTGGTTCCAACGCGTCGGGCGGGTGCGATGGTGGTCGTTTCTGGCGGCAAATTGTTACTCTATCTGGCGCAAGGTGGCAAAAAAATGCTGGTCTGGCAGGAAAAAGAGGAATTACTCGCCCCAGAGGTTTTCCACGCGCTGa                                    <  1:769795/149‑1 (MQ=255)
                                                                                                                   gCAACGCTGGTTCCAACGCGTCGGGCGGGTGCGCTGGTGGTCGTTTCTGGCGGCAAATTGTTACTCTATCTGGCGCAAGGTGGCAAAAAAATGCTGGTCTGGCAGGAAAAAGAGGAATTACTCGCCCCAGAGGTTTTCCACGCGCTGAc                                   <  2:879392/149‑1 (MQ=255)
                                                                                                                          tGGTTCCAACGCGTCGGGCGGGTGCGCTGGTGGTCGTTTCTGGCGGCAAATTGTTACTCTATCTGGCGCAAGGTGGCAAAAAAATGCTGGTCTGGCAGGAAAAAGAGGAATTACTCGCCCCAGAGGTTTTCCACGCGCTGACTACCGCa                            >  1:4421448/1‑149 (MQ=255)
                                                                                                                             ttCCAACGCGTCGGGCGGGTGCGCTGGTGGTCGTTTCTGGCGGCAAATTGTTACTCTATCTGGCGCAAGGTGGCAAAAAAATGCTGGTCTGGCAGGAAAAAGAGGAATTACTCGCCCCAGAGGTTTTCCACGCGCTGACTACCGCACTg                         >  2:4606585/1‑149 (MQ=255)
                                                                                                                             ttCCAACGCGTCGGGCGGGTGCGATGGTGGTCGTTTCTGGCGGCAAATTGTTACTCTATCTGGCGCAAGGTGGCAAAAAAATGCTGGTCTGGCAGGAAAAAGAGGAATTACTCGCCCCAGAGGTTTTCCACGCGCTGACTACCGCACTg                         >  1:3715207/1‑149 (MQ=255)
                                                                                                                              tCCAACGCGTCGGGCGGGTGCGCTGGTGGTCGTTTCTGGCGGCAAATTGTTACTCTATCTGGCGCAAGGTGGCAAAAAAATTCTGGTCTGGCAGGAAAAAGAGGAATTACTCGCCCCAGAGGTTTTCCACGCGCTGACTACCGCACTGc                        >  2:1609073/1‑149 (MQ=255)
                                                                                                                                cAACGCGTCGGGCGGGTGCGCTGGTGGTCGTTTCTGGCGGCAAATTGTTACTCTATCTGGCGCAAGGTGGCAAAAAAATGCTGGTCTGGCAGGAAAAaga                                                                       <  3:2503013/100‑1 (MQ=255)
                                                                                                                                 aaCGCGTCGGGCGGGTGCGCTGGTGGTCGTTTCTGGCGGCAAATTGTTACTCTATCTGGCGCAAGGTGGCAAAAAAATGCTGGTCTGGCAGGAAAAAGAGGAATTACTCGCCCCAGAGGTTTTCCACGCGCTGACTACCGCACTGCGTc                     >  1:1514199/1‑149 (MQ=255)
                                                                                                                                 aaCGCGTCGGGCGGGTGCGCTGGTGGTCGTTTCTGGCGGCAAATTGTTACTCTATCTGGCGCAAGGTGGCAAAAAAATGCTGGTCTGGCAGGAAAAAGAGGAATTACTCGCCCCAGAGGTTTTCCACGCGCTGACTACCGCACTGCGTc                     <  2:244408/149‑1 (MQ=255)
                                                                                                                                 aaCGCGTCGGGCGGGTGCGCTGGTGGTCGTTTCTGGCGGCAAATTGTTACTCTATCTGGCGCAAGGTGGCAAAAAAATGCTGGTCTGGCAGGAAAAAGAGGAATTACTCGCCCCAGAGGTTTTCCACGCGCTGACTACCGCACTGCGTc                     >  4:6542503/1‑149 (MQ=255)
                                                                                                                                  aCGCGTCGGGCGGGTGCGCTGGTGGTCGTTTCTGGCGGCAAATTGTTACTCTATCTGGCGCAAGGTGGCAAAAAAATGCTGGTCTGGCAGGAAAAAGAGGAATTACTCGCCCCAGAGGTTTTCCACGCGCTGACTACCGCACTGCGTcg                    >  2:3819751/1‑149 (MQ=255)
                                                                                                                                  aCGCGTCGGGCGGGTGCGCTGGTGGTCGTTTCTGGCGGCAAATTGTTACTCTATCTGGCGCAAGGTGGCAAAAAAATGCTGGTCTGGCAGGAAAAAGAGGAATTACTCGCCCCAGAGGTTTTCCACGCGCTGACTACCGCACTGCGTcg                    >  2:6002166/1‑149 (MQ=255)
                                                                                                                                   cgcgTCGGGCGGGTGCGCTGGTGGTCGTTTCTGGCGGCAAATTGTTACTCTATCTGGCGCAAGGTGGCAAAAAAATGCTGGTCTGGCAGGAAAAAGAGGAATTACTCGCCCCAGAGGTTTTCCACGCGCTGACTACCGCACTGCGTcgc                   <  2:3533771/149‑1 (MQ=255)
                                                                                                                                       tCGGGCGGGTGCGATGGTGGTCGTTTCTGGCGGCAAATTGTTACTCTATCTGGCGCAAGGTGGCAAAAAAATGCTGGTCTGGCAGGAAAAAGAGGAATTACTCGCCCCAGAGGTTTTCCACGCGCTGACTACCGCACTGCGTCGCGAAc               <  2:106910/149‑1 (MQ=255)
                                                                                                                                           gcgggTGCGCTGGTGGTCGTTTCTGGCGGCAAATTGTTACTCTATCTGGCGCAAGGTGGCAAAAAAATGCTGGTCTGGCAGGAAAAAGAGGAATTACTCGCCCCAGAGGTTTTCCACGCGCTGACTACCGCACTGCGTCGCGAACCACg           >  3:597601/1‑149 (MQ=255)
                                                                                                                                                 gcgcTGGTGGTCGTTTCTGGCGGCAAATTGTTACTCTATCTGGCGCAAGGTGGCAAAAAAATGCTg                                                                                        >  2:1434714/1‑66 (MQ=255)
                                                                                                                                                   gcTGGTGGTCGTTTCTGGCGGCAAATTGTTACTCTATCTGGCGCAAGGTGGCAAAAAAATGCTGGTCTGGCAGGAAAAAGAGGAATTACTCGCCCCAGAGGTTTTCCACGCGCTGACTACCGCACTGCGTCGCGAACCACGGCTGCGCt   <  1:987878/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
ATTGCGCGATCTGGCGACACAAGCCACGCAAACGCGCCACTATACACCAGTGGCGCTCTCTGCCAACGATCCGGCTAATGTGTGGGGAAATCTTCTGCCCTGGCCTGCACATCCGGCAACGCTGGTTCCAACGCGTCGGGCGGGTGCGCTGGTGGTCGTTTCTGGCGGCAAATTGTTACTCTATCTGGCGCAAGGTGGCAAAAAAATGCTGGTCTGGCAGGAAAAAGAGGAATTACTCGCCCCAGAGGTTTTCCACGCGCTGACTACCGCACTGCGTCGCGAACCACGGCTGCGCTT  >  NC_000913/1733301‑1733597

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