Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,912,748 T→G 7.2% intergenic (‑172/+20) htpX ← / ← prc putative endopeptidase/carboxy‑terminal protease for penicillin‑binding protein 3

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,912,7480TG 7.2% 284.8 / 14.2 83intergenic (‑172/+20)htpX/prcputative endopeptidase/carboxy‑terminal protease for penicillin‑binding protein 3
Reads supporting (aligned to +/- strand):  ref base T (34/43);  new base G (2/4);  total (36/47)
Fisher's exact test for biased strand distribution p-value = 6.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.87e-01

AAAACGCAATATGCGATACCCACATCGTATGGGTTACGCGACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATCCGGCTCCTGGTAATCTTTCGGTAGATCATCCAGTT  >  NC_000913/1912602‑1912895
                                                                                                                                                  |                                                                                                                                                   
aaaaCGCAATATGCGATACCCACATCGTATGGGTTACGCGACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAAttt                                                                                                                                                   >  1:729290/1‑149 (MQ=255)
  aaCGCAATATGCGATACCCACATCGTATGGGTTACGCGACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCt                                                                                                                                                 >  2:45963/1‑149 (MQ=255)
  aaCGCAATATGCGATACCCACATCGTATGGGTTACGCGACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCt                                                                                                                                                 <  1:807992/149‑1 (MQ=255)
      cAATATGCGATACCCACATCGTATGGGTTACGCGACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTtgtg                                                                                                                                             <  1:647064/149‑1 (MQ=255)
        atatGCGATACCCACATCGTATGGGTTACGCGACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGcc                                                                                                                                           <  3:789061/149‑1 (MQ=255)
            gCGATACCCACATCGTATGGGTTACGCGACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCtgat                                                                                                                                       <  3:640963/149‑1 (MQ=255)
                tACCCACATCGTATGGGTTACGCGACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCtgattgat                                                                                                                                   <  2:735356/149‑1 (MQ=255)
                  cccACATCGTATGGGTTACGCGACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGatat                                                                                                                                 <  3:216557/149‑1 (MQ=255)
                   ccACATCGTATGGGTTACGCGACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATAtt                                                                                                                                >  3:129571/1‑149 (MQ=255)
                    cacaTCGTATGGGTTACGCGACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAAGAAATTGTGCCtga                                                                                                                                        <  2:194927/140‑1 (MQ=255)
                    cacaTCGTATGGGTTACGCGACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAAGAAATTGTGCCtga                                                                                                                                        >  1:194927/1‑140 (MQ=255)
                     acaTCGTATGGGTTACGCGACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTAc                                                                                                                              >  1:589989/1‑149 (MQ=255)
                          gTATGGGTTACGCGACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGAc                                                                                                                         <  1:721671/149‑1 (MQ=255)
                                gTTACGCGACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGc                                                                                                                   >  2:679094/1‑149 (MQ=255)
                                  tACGCGACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCgg                                                                                                                 <  4:52272/149‑1 (MQ=255)
                                   aCGCGACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCggg                                                                                                                <  2:589989/149‑1 (MQ=255)
                                       gACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGcc                                                                                                                                           >  4:619929/1‑118 (MQ=255)
                                       gACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGcc                                                                                                                                           <  3:619929/118‑1 (MQ=255)
                                       gACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGgttgt                                                                                                            >  3:406861/1‑149 (MQ=255)
                                       gACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAAGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGgttgt                                                                                                            >  2:659187/1‑149 (MQ=255)
                                         cTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTc                                                                                                          <  3:171349/149‑1 (MQ=255)
                                         cTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTc                                                                                                          >  4:132081/1‑149 (MQ=255)
                                            ttttCAAGGCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCcgc                                                                                                       <  2:687378/149‑1 (MQ=255)
                                               tCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCggg                                                                                                    >  4:553267/1‑149 (MQ=255)
                                                  aGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCt                                                                                                 <  4:129571/149‑1 (MQ=255)
                                                   gTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTg                                                                                                >  1:8233/1‑149 (MQ=255)
                                                   gTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTg                                                                                                >  3:109790/1‑149 (MQ=255)
                                                        gATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCttt                                                                                           >  2:546677/1‑149 (MQ=255)
                                                           aGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTc                                                                                        <  1:546677/149‑1 (MQ=255)
                                                               cGTAGGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAAGAAATTGTGCCt                                                                                                                                          >  3:259793/1‑95 (MQ=255)
                                                               cGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAAGAAATTGTGCCt                                                                                                                                          <  4:259793/95‑1 (MQ=255)
                                                                 tAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCtga                                                                                                                                        <  3:187093/95‑1 (MQ=255)
                                                                 tAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCtga                                                                                                                                        >  4:187093/1‑95 (MQ=255)
                                                                    gTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGc                                                                               >  1:673355/1‑149 (MQ=255)
                                                                     tCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAAGAAATTGTGCCtga                                                                                                                                        <  3:700734/91‑1 (MQ=38)
                                                                        ccAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGa                                                                           >  4:106359/1‑149 (MQ=255)
                                                                           gAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCg                                                                        >  3:421950/1‑149 (MQ=255)
                                                                                  acaACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGa                                                                 >  4:169242/1‑149 (MQ=255)
                                                                                  acaACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGa                                                                 >  1:504016/1‑149 (MQ=255)
                                                                                  acaACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGa                                                                 >  3:699424/1‑149 (MQ=255)
                                                                                      cTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATAtt                                                             >  2:177009/1‑149 (MQ=255)
                                                                                      cTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATAtt                                                             <  3:434537/149‑1 (MQ=255)
                                                                                      cTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATAtt                                                             >  3:630742/1‑149 (MQ=255)
                                                                                       tttACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTc                                                            >  2:145907/1‑149 (MQ=255)
                                                                                       tttACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTc                                                            >  4:364810/1‑149 (MQ=255)
                                                                                        ttACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCa                                                           <  4:406861/149‑1 (MQ=255)
                                                                                          aCATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCAcc                                                         <  1:10985/149‑1 (MQ=255)
                                                                                            aTTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGt                                                       <  4:421950/149‑1 (MQ=255)
                                                                                                    gCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCa                                               <  3:553267/149‑1 (MQ=255)
                                                                                                      aTCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGa                                             <  2:309532/149‑1 (MQ=255)
                                                                                                       tCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGc                                                                                                                                            <  3:349940/53‑1 (MQ=255)
                                                                                                       tCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGc                                                                                                                                            >  4:349940/1‑53 (MQ=255)
                                                                                                         tGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATaa                                          <  3:132081/149‑1 (MQ=255)
                                                                                                         tGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATaa                                          <  4:109790/149‑1 (MQ=255)
                                                                                                                     tCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATCCGGCTcc                              >  3:440206/1‑149 (MQ=255)
                                                                                                                     tCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATCCGGCTcc                              >  1:326803/1‑149 (MQ=255)
                                                                                                                       ttGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATCCGGCTCCTg                            >  3:640349/1‑149 (MQ=255)
                                                                                                                         gTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCAc                                                          <  2:367106/117‑1 (MQ=255)
                                                                                                                         gTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCAc                                                          >  1:367106/1‑117 (MQ=255)
                                                                                                                          tCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATCCGGCTCCTGGTa                         <  1:679094/149‑1 (MQ=255)
                                                                                                                            gCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATCCGGCTCCTGGTAAt                       <  3:106359/149‑1 (MQ=255)
                                                                                                                             cTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATCCGGCTCCTGGTAATc                      <  1:659187/149‑1 (MQ=255)
                                                                                                                              tGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATCCGGCTCCTGGTAATCt                     <  2:673355/149‑1 (MQ=255)
                                                                                                                               gTTACAAAATCAGGCCAAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGCGCGATATTCACCGTCTCATCCAGATAAGGATcc                                    >  4:607914/1‑133 (MQ=255)
                                                                                                                               gTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATCCGGCTCCTGGTAATCtt                    >  4:134107/1‑149 (MQ=255)
                                                                                                                               gTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGCGCGATATTCACCGTCTCATCCAGATAAGGATcc                                    <  3:607914/133‑1 (MQ=255)
                                                                                                                                 tAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATCCGGCTCCTGGTAATCTTTc                  <  2:134949/149‑1 (MQ=255)
                                                                                                                                 tAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATCCGGCTCCTGGTAATCTTTc                  <  4:640349/149‑1 (MQ=255)
                                                                                                                                  aaaaaaTCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATCCGGCTCCTGGTAATCTTTCg                 <  2:326803/149‑1 (MQ=255)
                                                                                                                                  aaaaaaTCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATCCGGCTCCTGGTAATCTTTCg                 <  4:440206/149‑1 (MQ=255)
                                                                                                                                   aaaaaTCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATCCGGCTCCTGGTAATCTTTCgg                <  1:579303/149‑1 (MQ=255)
                                                                                                                                    tcaaTCAGGCACAAGAAATTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATaa                                          <  2:270101/120‑1 (MQ=39)
                                                                                                                                    aaaaTCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATCCGGCTCCTGGTAATCTTTCGGt               >  1:803590/1‑149 (MQ=255)
                                                                                                                                     aaaTCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATCCGGCTCCTGGTAATCTTTCGGTa              <  2:803590/149‑1 (MQ=255)
                                                                                                                                        tCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGCGCGGGTTGTTCCGCGGGTCTGGCTTTTTCaa                                                                                      >  3:207758/1‑74 (MQ=255)
                                                                                                                                        tCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCaa                                                                                      <  4:207758/74‑1 (MQ=255)
                                                                                                                                        tCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATCCGGCTCCTGGTAATCTTTCGGTAGat           >  4:634603/1‑149 (MQ=255)
                                                                                                                                        tCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATCCGGCTCCTGGTAATCTTTCGGTAGat           <  1:145907/149‑1 (MQ=255)
                                                                                                                                        tCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATCCGGCTCCTGGTAATCTTTCGGTAGat           <  3:364810/149‑1 (MQ=255)
                                                                                                                                        tCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATCCGGCTCCTGGTAATCTTTCGGTAGat           <  2:504016/149‑1 (MQ=255)
                                                                                                                                              acaATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGt                                                                     >  2:451896/1‑85 (MQ=255)
                                                                                                                                              acaATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGt                                                                     <  1:451896/85‑1 (MQ=255)
                                                                                                                                              acaATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATCCGGCTCCTGGTAATCTTTCGGTAGATCATCCa     <  1:45963/149‑1 (MQ=255)
                                                                                                                                                 aTTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATCCGGCTCCTGGTAATCTTTCGGTAGATCATCCAGtt  <  4:699424/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
AAAACGCAATATGCGATACCCACATCGTATGGGTTACGCGACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATCCGGCTCCTGGTAATCTTTCGGTAGATCATCCAGTT  >  NC_000913/1912602‑1912895

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