Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,003,800 G→T 8.8% R127R (CGG→AGG)  rarD ← putative chloramphenical resistance permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,003,8000GT 8.8% 232.9 / 7.9 68R127R (CGG→AGG) rarDputative chloramphenical resistance permease
Reads supporting (aligned to +/- strand):  ref base G (25/37);  new base T (3/3);  total (28/40)
Fisher's exact test for biased strand distribution p-value = 6.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.31e-01

GCAATCTTCTTGCGTACCAGACCGTAGAAGGCAAAACTAAATGCCAGTCCCAGCGCGATAATAGGTAGCGAACCAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACG  >  NC_000913/4003664‑4003939
                                                                                                                                        |                                                                                                                                           
gCAATCTTCTTGCGTACCAGACCGTAGAAGGCAAAACTAAATGCCAGTCCCAGCGCGATAATAGGTAGCGAACCAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCcggcg                                                                                                                                          <  1:909066/140‑1 (MQ=255)
              tACCAGACCGTAGAAGGCAAAACTAAATGCCAGTCCCAGCGCGATAATAGGTAGCGAACCAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGACTCGCTCGCCGa                                                                                                                            >  1:233130/1‑140 (MQ=255)
              tACCAGACCGTAGAAGGCAAAACTAAATGCCAGTCCCAGCGCGATAATAGGTAGCGAACCAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGa                                                                                                                            >  1:232600/1‑140 (MQ=255)
                      cGTAGAAGGCAAAACTAAATGCCAGTCCCAGCGCGATAATAGGTAGCGAACCAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAatc                                                                                                                    <  1:178973/140‑1 (MQ=255)
                           aaGGCAAAACTAAATGCCAGTCCCAGCGCGATAATAGGTAGCGAACCAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATccc                                                                                                               <  1:1106181/140‑1 (MQ=255)
                           aaGGCAAAACTAAATGCCAGTCCCAGCGCGATAATAGGTAGCGAACCAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATccc                                                                                                               <  1:164757/140‑1 (MQ=255)
                              acAAAAATAAAAGACAGTCCCAGCGCGATAATATTTATATAAACAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAg                                                                                                            >  1:265430/2‑140 (MQ=255)
                                aaaaCTAAATGCCAGTCCCAGCGCGATAATAGGTAGCGAACCAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCTAGGAAAATCATCCCCAGca                                                                                                          >  1:173967/1‑140 (MQ=255)
                                aaaaCTAAATGCCAGTCCCAGCGCGATAATAGGTAGCGAACCAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCTAGGAAAATCATCCCCAGca                                                                                                          >  1:174010/1‑140 (MQ=255)
                                aaaaCTAAATGCCAGTCCCAGCGCGATAATAGGTAGCGAACCAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGca                                                                                                          <  2:955915/140‑1 (MQ=255)
                                aaaaCTAAATGCCAGTCCCAGCGCGATAATAGGTAGCGAACCAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGca                                                                                                          >  2:658078/1‑140 (MQ=255)
                                aaaaCTAAATGCCAGTCCCAGCGCGATAATAGGTAGCGAACCAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGca                                                                                                          >  1:399434/1‑140 (MQ=255)
                                     aaaaTGCCAGTCCCAGCGCGATAATAGGTAGCGAACCAAAAGTCCACAGCTGGACTAACACGCCTCATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCGCGCCGAGGAAACTCATCCCCAGCACAATg                                                                                                     <  2:673012/139‑1 (MQ=255)
                                        aTGCCAGTCCCAGCGCGATAATAGGTAGCGAACCAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTc                                                                                                  >  2:825822/1‑140 (MQ=255)
                                        aTGCCAGTCCCAGCGCGATAATAGGTAGCGAACAAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTc                                                                                                  >  1:670771/1‑140 (MQ=255)
                                                   aGCGCGATAATAGGTAGCGAACCAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGtt                                                                                       >  1:43264/1‑140 (MQ=255)
                                                     cgcgATAATAGGTAGCGAACCAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTtaa                                                                                     <  2:73087/140‑1 (MQ=255)
                                                                    cGAACCAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGTCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGc                                                                      >  1:958893/1‑140 (MQ=255)
                                                                    cGAACCAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGc                                                                      >  1:959024/1‑140 (MQ=255)
                                                                     gAACCAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCt                                                                     <  1:1221400/140‑1 (MQ=255)
                                                                          aacaGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCtt                                                                <  1:1192587/137‑1 (MQ=255)
                                                                          aaaaGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCtt                                                                <  1:313830/140‑1 (MQ=255)
                                                                          aaaaGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCtt                                                                <  1:261408/140‑1 (MQ=255)
                                                                            aaGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTcc                                                              <  1:708118/140‑1 (MQ=255)
                                                                            aaGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTcc                                                              <  1:292259/140‑1 (MQ=255)
                                                                            aaGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTcc                                                              <  1:693488/140‑1 (MQ=255)
                                                                            aaGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTcc                                                              <  1:129280/140‑1 (MQ=255)
                                                                            aaGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTcc                                                              <  1:1191118/140‑1 (MQ=255)
                                                                             aGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCa                                                             <  2:431448/140‑1 (MQ=255)
                                                                                     gCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATgg                                                     <  2:636385/140‑1 (MQ=255)
                                                                                     gCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATgg                                                     <  2:160940/140‑1 (MQ=255)
                                                                                     gCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATgg                                                     <  2:636493/140‑1 (MQ=255)
                                                                                        ggACTAACACGCCACATATCGCCAGTATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCTCCGAGGAAAATCATCCCCAGCTCACTGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGa                                                  <  2:332949/140‑1 (MQ=255)
                                                                                        ggACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGa                                                  <  2:1156456/140‑1 (MQ=255)
                                                                                               cacGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCCAGGCTCGCATCCAGCATATGGTGATTGTTCa                                           >  2:431886/1‑140 (MQ=255)
                                                                                               cacGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCa                                           >  2:432069/1‑140 (MQ=255)
                                                                                                   ccACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTTATAAAGTACCCAAGGCTCGCGTCCCGCATATGGAGATTGTTCACCGc                                       >  2:1205505/1‑140 (MQ=255)
                                                                                                   ccACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGc                                       >  1:1122534/1‑140 (MQ=255)
                                                                                                   ccACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGc                                       >  2:1163673/1‑140 (MQ=255)
                                                                                                   ccACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGGTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGc                                       >  2:1205265/1‑140 (MQ=255)
                                                                                                    cacaTATCGCCAGAATCACCGCCAGCCATTGCATCCTGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGcc                                      >  2:644745/1‑140 (MQ=255)
                                                                                                    cacaTATCGCCAGAATCACCGCCAGCCATTGCATCCTGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGcc                                      >  2:641959/1‑140 (MQ=255)
                                                                                                    cacaTATCGCCAGAATCACCGCCAGCCATTGCATCCTGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGcc                                      >  2:644882/1‑140 (MQ=255)
                                                                                                    cacaTATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGcc                                      <  1:697658/140‑1 (MQ=255)
                                                                                                    cacaTATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGcc                                      <  2:67309/140‑1 (MQ=255)
                                                                                                    cacaTATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGcc                                      >  1:194272/1‑140 (MQ=255)
                                                                                                    cacaTATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGcc                                      >  1:194339/1‑140 (MQ=255)
                                                                                                    cacaTATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGAGTGTTAATTAATTAAACAATTCTCGCTTCCAGCTTATGTTGATttttgtggtgt                                      <  1:27797/140‑8 (MQ=255)
                                                                                                         atCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGAt                                 <  1:1097472/140‑1 (MQ=255)
                                                                                                         atCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGAt                                 <  2:117292/140‑1 (MQ=255)
                                                                                                         atCGCCAGAATCACCGACAGACATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGAt                                 <  2:1045077/140‑1 (MQ=255)
                                                                                                          tCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATa                                <  2:171640/140‑1 (MQ=255)
                                                                                                          tCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATa                                <  2:172075/140‑1 (MQ=255)
                                                                                                            gCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATaaa                              <  1:692816/140‑1 (MQ=255)
                                                                                                             ccAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAAc                             >  1:439465/1‑140 (MQ=255)
                                                                                                              cAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACa                            <  1:771664/140‑1 (MQ=255)
                                                                                                              cAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACa                            <  1:143259/140‑1 (MQ=255)
                                                                                                              cAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACa                            <  2:202838/140‑1 (MQ=255)
                                                                                                                 aaTCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTa                         >  2:107959/1‑140 (MQ=255)
                                                                                                                 aaTCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTa                         >  2:107973/1‑140 (MQ=255)
                                                                                                                            gccaTTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGcca              <  1:69132/140‑1 (MQ=255)
                                                                                                                            gccaTTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGcca              <  1:68721/140‑1 (MQ=255)
                                                                                                                                 tGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAAt         <  2:62261/140‑1 (MQ=255)
                                                                                                                                 tGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAAt         >  2:1098905/1‑140 (MQ=255)
                                                                                                                                 tGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAAt         >  2:889548/1‑140 (MQ=255)
                                                                                                                                     tCCTGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGc     <  1:644882/140‑1 (MQ=255)
                                                                                                                                     tCCTGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGc     <  1:644745/140‑1 (MQ=255)
                                                                                                                                     tCCTGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGc     <  1:641959/140‑1 (MQ=255)
                                                                                                                                     tCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGc     <  2:431949/140‑1 (MQ=255)
                                                                                                                                       cggcggAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCAc   >  1:556834/1‑140 (MQ=255)
                                                                                                                                       cggcggAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCAc   >  1:555751/1‑140 (MQ=255)
                                                                                                                                        ggcggAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACg  >  2:1122524/1‑140 (MQ=255)
                                                                                                                                        |                                                                                                                                           
GCAATCTTCTTGCGTACCAGACCGTAGAAGGCAAAACTAAATGCCAGTCCCAGCGCGATAATAGGTAGCGAACCAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACG  >  NC_000913/4003664‑4003939

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