Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,964,929 C→T E234E (GAG→GAA)  rhlB ← ATP‑dependent RNA helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,964,9290CT99.1% 403.7 / NA 116E234E (GAG→GAA) rhlBATP‑dependent RNA helicase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (52/63);  minor base G (1/0);  total (53/63)
Fisher's exact test for biased strand distribution p-value = 4.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 117 total aligned reads displayed.

GACACGATGACCATCTGCTGCCAGGTGGCCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCA‑TTTTTTCTTCGTTAGAAGGGTAGAAAAGCTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGC  >  NC_000913/3964782‑3965073
                                                                                                                                                    |                                                                                                                                                
gACACGATGACCATCTGCTGCCATGTGGCCCCAGATCTTTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCA‑TTTTTTCTTCGTTAGAAGGGTAGAAAAGttc                                                                                                                                                <  3:444393/150‑1 (MQ=255)
     gATGACCATCTGCTGCCAGGTGGCCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCA‑TTTTTTCTTCGTTAGAAGGGTAGAAAAGTTcttctt                                                                                                                                           <  2:138512/150‑1 (MQ=255)
     gATGACCATCTGCTGCCAGGTGGCCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCA‑TTTTTTCTTCGTTAGAAGGGTAGAAAAGTTcttctt                                                                                                                                           <  1:975839/150‑1 (MQ=255)
        gACCATCTGCTGCCAGGTGGCCCCAGATCTCTTCACAACGGTGTTTTGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCA‑TTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTaa                                                                                                                                        <  3:126232/150‑1 (MQ=255)
          ccATCTGCTGCCAGGTGGCCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCA‑TTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATa                                                                                                                                      >  1:339187/1‑150 (MQ=255)
           cATCTGCTGCCAGGTGGCCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCA‑TTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATAc                                                                                                                                     >  1:872220/1‑150 (MQ=255)
             tCTGCTGCCAGGTGGCCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCA‑TTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACgg                                                                                                                                   >  3:123404/1‑150 (MQ=255)
              ctgctgCCAGGTGGCCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCA‑TTTTTTATTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGt                                                                                                                                  <  1:351559/150‑1 (MQ=255)
               tgctgcCAGGTGGCCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCA‑TTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTg                                                                                                                                 <  3:339217/150‑1 (MQ=255)
                     cAGGTGGCCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCA‑TTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGt                                                                                                                           >  1:958274/1‑150 (MQ=255)
                          ggCCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCA‑TTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCt                                                                                                                      <  3:336044/150‑1 (MQ=255)
                            ccccAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCA‑TTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGt                                                                                                                    <  2:546240/150‑1 (MQ=255)
                             cccAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCA‑TTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGtt                                                                                                                   >  2:411377/1‑150 (MQ=255)
                             cccAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCA‑TTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGtt                                                                                                                   >  2:615755/1‑150 (MQ=255)
                              ccAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCA‑TTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTc                                                                                                                  >  2:955513/1‑150 (MQ=255)
                                aGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCA‑TTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCg                                                                                                                <  2:590109/150‑1 (MQ=255)
                                    ctctTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCA‑TTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTc                                                                                                            <  4:495209/150‑1 (MQ=255)
                                         cacaACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCA‑TTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACtt                                                                                                       >  1:59082/1‑150 (MQ=255)
                                          acaACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCA‑TTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTc                                                                                                      >  1:815769/1‑150 (MQ=255)
                                            aaCGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCA‑TTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCaa                                                                                                    <  1:245952/150‑1 (MQ=255)
                                            aaCGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCA‑TTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCaa                                                                                                    <  1:676711/150‑1 (MQ=255)
                                            aaCGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCA‑TTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCaa                                                                                                    >  4:61811/1‑150 (MQ=255)
                                            aaCGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTATTCTAGCAGCGTTTGCAGCAAACGCA‑TTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCaa                                                                                                    >  1:1128571/1‑150 (MQ=255)
                                             aCGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCA‑TTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAat                                                                                                   <  1:1583/150‑1 (MQ=255)
                                             aCGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCA‑TTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAat                                                                                                   >  2:789369/1‑150 (MQ=255)
                                                gtgtTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCA‑TTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAatata                                                                                                >  1:716335/1‑150 (MQ=255)
                                                   tttGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCA‑TTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCTGGTTCCACTTGAATATATTc                                                                                             <  2:477516/150‑1 (MQ=255)
                                                   tttGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCA‑TTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTc                                                                                             >  1:247887/1‑150 (MQ=255)
                                                   tttGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCA‑TTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTc                                                                                             <  3:463334/150‑1 (MQ=255)
                                                   tttGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCA‑TTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTc                                                                                             >  1:113201/1‑150 (MQ=255)
                                                   tttGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCA‑TTTTTTCTTCGTTAGAAGGGTAGAAAAGGTCTGCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTc                                                                                             >  1:48032/1‑150 (MQ=255)
                                                      ggTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCA‑TTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGc                                                                                          >  4:135125/1‑150 (MQ=255)
                                                         gtTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCA‑TTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCAtt                                                                                       <  1:659791/150‑1 (MQ=255)
                                                         gtTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCA‑TTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCAtt                                                                                       <  3:195875/150‑1 (MQ=255)
                                                             gCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCA‑TTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGcc                                                                                                                              >  3:51262/1‑107 (MQ=255)
                                                             gCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCA‑TTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGcc                                                                                                                              <  4:51262/107‑1 (MQ=255)
                                                               gAAAATAATCGCTCTGTCTGGCCGCTCTTCTTCGATCAGCGTTTGCAGCAAACGCA‑TTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCAt                                                                                 <  1:881375/150‑1 (MQ=255)
                                                               gAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCA‑TTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCAt                                                                                 >  1:453769/1‑150 (MQ=255)
                                                                       tCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCACTATATTCGGCATTGTTCATCTGCTCg                                                                          >  2:1125440/1‑150 (MQ=255)
                                                                         gCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCA‑TTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAAc                                                                       <  1:741395/150‑1 (MQ=255)
                                                                             tgtctgGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCA‑TTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCa                                                                   <  2:59082/150‑1 (MQ=255)
                                                                              gtctgGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCA‑TTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAg                                                                  <  2:727668/150‑1 (MQ=255)
                                                                              gtctgGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCA‑TTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAg                                                                  <  4:153291/150‑1 (MQ=255)
                                                                              gtctgGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCA‑TTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAg                                                                  <  2:958274/150‑1 (MQ=255)
                                                                              gtctgGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCA‑TTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAg                                                                  <  1:585384/150‑1 (MQ=255)
                                                                                 tgGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCA‑TTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGtt                                                                                    >  1:1047505/1‑129 (MQ=255)
                                                                                 tgGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCA‑TTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGtt                                                                                    <  2:1047505/129‑1 (MQ=255)
                                                                                 tgGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCA‑TTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTc                                                               <  1:500096/150‑1 (MQ=255)
                                                                                 tgGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCA‑TTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTc                                                               <  1:500298/150‑1 (MQ=255)
                                                                                  ggCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCA‑TTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCa                                                              >  1:772348/1‑150 (MQ=255)
                                                                                   gCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCA‑TTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCAc                                                             <  2:430603/150‑1 (MQ=255)
                                                                                        tcttctTCGATCAGCGTTTGCAGCAAACGCA‑TTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTAcc                                                        >  1:999395/1‑150 (MQ=255)
                                                                                        tcttctTCGATCAGCGTTTGCAGCAAACGCA‑TTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTAcc                                                        >  3:166055/1‑150 (MQ=255)
                                                                                        tcttctTCGATCAGCGTTTGCAGCAAACGCA‑TTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTAcc                                                        >  1:774446/1‑150 (MQ=255)
                                                                                        tcttctTCGATCAGCGTTTGCAGCAAACGCA‑TTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTAcc                                                        >  1:14415/1‑150 (MQ=255)
                                                                                        tcttctTCGATCAGCGTTTGCAGCAAACGCA‑TTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTAcc                                                        >  2:236017/1‑150 (MQ=255)
                                                                                        tcttctTCGATCAGCGTTTGCAGCAAACGCA‑TTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTAcc                                                        >  1:14839/1‑150 (MQ=255)
                                                                                            cttcGATCAGCGTTTGCAGCAAACGCA‑TTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGt                                                    <  4:245411/150‑1 (MQ=255)
                                                                                            cttcGATCAGCGTTTGCAGCAAACGCA‑TTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGt                                                    <  3:321394/150‑1 (MQ=255)
                                                                                            cttcGATCAGCGTTTGCAGCAAACGCA‑TTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGt                                                    >  1:1038855/1‑150 (MQ=255)
                                                                                            cttcGATCAGCGTTTGCAGCAAACGCA‑TTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGt                                                    <  2:872220/150‑1 (MQ=255)
                                                                                             ttcGATCAGCGTTTGCAGCAAACGCA‑TTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGc                                                                                                                               >  2:856777/1‑74 (MQ=255)
                                                                                                gATCAGCGTTTGCAGCAAACGCA‑TTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGa                                                >  3:336472/1‑150 (MQ=255)
                                                                                                  tCAGCGTTTGCAGCAAACGCA‑TTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGaaa                                              >  1:327105/1‑150 (MQ=255)
                                                                                                     gCGTTTGCAGCAAACGCA‑TTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCg                                           <  2:713895/150‑1 (MQ=255)
                                                                                                     gCGTTTGCAGCAAACGCA‑TTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCg                                           >  1:1011050/1‑150 (MQ=255)
                                                                                                      cGTTTGCAGCAAACGCA‑TTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGt                                          <  3:52609/150‑1 (MQ=255)
                                                                                                      cGTTTGCAGCAAACGCA‑TTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGt                                          <  1:824430/150‑1 (MQ=255)
                                                                                                      cGTTTGCAGCAAACGCA‑TTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGt                                          <  2:1112079/150‑1 (MQ=255)
                                                                                                       gTTTGCAGCAAACGCA‑TTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTg                                         >  2:458044/1‑150 (MQ=255)
                                                                                                       gTTTGCAGCAAACGCA‑TTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCAGTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTg                                         <  1:436275/150‑1 (MQ=255)
                                                                                                       gTTTGCAGCAAACGCA‑TTTTTTCTTCGTTAGAAGGGTAGAAAAGTTATTCTTTAATACGGTGGCACGTTTTCTAATCCGGCTCCACCTCAATATATTCGGCATTGTTCATCTGCTCGAACGACAGTTCAGGTCCCCGGTTCCAAAGCGTg                                         >  2:741565/1‑150 (MQ=255)
                                                                                                       gTTTGCAGCAAACGCA‑TTTTTTCATCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAGCGCCAGTGCACGTAACCGGTACGCAAGCGTg                                         >  2:741436/1‑150 (MQ=255)
                                                                                                          tGCAGCAAACGCA‑TTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCg                                      <  1:801516/150‑1 (MQ=255)
                                                                                                            cagcaAACGCA‑TTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGga                                    <  2:327105/150‑1 (MQ=255)
                                                                                                            cagcaAACGCA‑TTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGga                                    <  2:247887/150‑1 (MQ=255)
                                                                                                              gcaAACGCA‑TTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGgaga                                  >  2:252942/1‑150 (MQ=255)
                                                                                                                   cGCA‑TTTTTTCTTCGTTTGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGc                             <  1:1000162/150‑1 (MQ=255)
                                                                                                                   cGCA‑TTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGc                             <  2:952806/150‑1 (MQ=255)
                                                                                                                      a‑TTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATCTTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATg                          <  2:285684/150‑1 (MQ=255)
                                                                                                                        ttttttCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGCTCATCTGCTCGAACGCCCGTTCACGTACCCGGTACGAAAGCGTGGCGGAGCACAGCATGt                         >  1:346460/1‑150 (MQ=255)
                                                                                                                           tttCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGa                      <  1:955513/150‑1 (MQ=255)
                                                                                                                            ttcttcGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAg                     <  2:252070/150‑1 (MQ=255)
                                                                                                                               ttcGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCGTTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGgcg                  >  2:1111963/1‑150 (MQ=255)
                                                                                                                                  gTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTg               <  1:154525/150‑1 (MQ=255)
                                                                                                                                  gTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTg               <  1:522124/150‑1 (MQ=255)
                                                                                                                                   ttAGAAGGGTAGAAACGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTgg              <  2:48032/150‑1 (MQ=255)
                                                                                                                                    tAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGt             <  3:92609/150‑1 (MQ=255)
                                                                                                                                       aaGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTg          >  2:829278/1‑150 (MQ=255)
                                                                                                                                       aaGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTg          >  1:771735/1‑150 (MQ=255)
                                                                                                                                         gggTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTc                                                                                   <  2:667262/75‑1 (MQ=255)
                                                                                                                                         gggTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTc                                                                                   >  1:667262/1‑75 (MQ=255)
                                                                                                                                              gAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTAc                                                  >  3:80153/1‑103 (MQ=255)
                                                                                                                                              gAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTAc                                                  >  3:79953/1‑103 (MQ=255)
                                                                                                                                              gAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTAc                                                  <  4:79953/103‑1 (MQ=255)
                                                                                                                                              gAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTAc                                                  <  4:80153/103‑1 (MQ=255)
                                                                                                                                              gAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAggcgg   >  2:1028260/1‑150 (MQ=255)
                                                                                                                                               aaaaGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTa                                                   >  1:671034/1‑101 (MQ=255)
                                                                                                                                               aaaaGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTa                                                   <  2:671034/101‑1 (MQ=255)
                                                                                                                                               aaaaGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAggcggc  >  2:954358/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
GACACGATGACCATCTGCTGCCAGGTGGCCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCA‑TTTTTTCTTCGTTAGAAGGGTAGAAAAGCTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGC  >  NC_000913/3964782‑3965073

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