Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,439,587 G→C 7.6% intergenic (‑80/+27) yhdN ← / ← rplQ DUF1992 family protein/50S ribosomal subunit protein L17

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,439,5870GC 7.6% 204.1 / 6.3 71intergenic (‑80/+27)yhdN/rplQDUF1992 family protein/50S ribosomal subunit protein L17
Reads supporting (aligned to +/- strand):  ref base G (36/29);  new base C (4/2);  total (40/31)
Fisher's exact test for biased strand distribution p-value = 6.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

TCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCGGCACGGCTCG  >  NC_000913/3439443‑3439735
                                                                                                                                                |                                                                                                                                                    
tCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCg                                                                                                                                                  >  2:581190/1‑149 (MQ=255)
tCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCg                                                                                                                                                  >  1:282253/1‑149 (MQ=255)
tCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCg                                                                                                                                                  >  2:188322/1‑149 (MQ=255)
tCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCg                                                                                                                                                  >  1:704172/1‑149 (MQ=255)
tCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCAATTGTAGATAAGTGGGGGTACTACGGGTATAAAAAAACCCGCCGGGGCg                                                                                                                                                  >  1:207330/1‑149 (MQ=255)
 taaccgcACCTTGCTCTGGCGCTTCTTCTAAATGGCGCTCTCCCCACTGGACAAGTAACCACATAAATACTCCAGGTGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCgg                                                                                                                                                 <  2:368454/143‑1 (MQ=255)
     cTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCCCGGCGGGttt                                                                                                                                             >  2:292247/1‑149 (MQ=255)
       cACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCCCGGCGGGttttt                                                                                                                                           <  2:55731/149‑1 (MQ=255)
       cACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCCCGGCGGGttttt                                                                                                                                           <  1:292247/149‑1 (MQ=255)
          cTTTCGTTCGCTCTTCTGCTATATGGCGACCGGCCCACTGGCCAAGTAACAACATAAATCCTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTa                                                                                                                                        <  2:148889/149‑1 (MQ=255)
          cTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTCACCACATAAATACTCCCGGGGATAAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGAAAACTACGGGTATAAAAAAACCAGCCGGGGAGGGtcttttga                                                                                                                                        >  2:362807/1‑147 (MQ=255)
          cTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTa                                                                                                                                        >  1:673436/1‑149 (MQ=255)
          cTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTa                                                                                                                                        >  2:641644/1‑149 (MQ=255)
          cTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTGTTTTa                                                                                                                                        >  1:97804/1‑149 (MQ=255)
          cTTTCGCTTGCGCTTCCGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTa                                                                                                                                        <  2:192488/149‑1 (MQ=255)
                cTTGCGCTTCTTCTATATGCCCCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCAGGGCGGGTTTTTTTACGTTGc                                                                                                                                  <  1:777916/149‑1 (MQ=255)
                              atatGGCGCTCTGCCCACAGGTAAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGCGCGGGTTTTTTTACGTTGCTTCAGATTACTCTg                                                                                                                    <  1:512372/149‑1 (MQ=255)
                               tatGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTgc                                                                                                                   >  1:75989/1‑149 (MQ=255)
                                atGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGCGGCGGGGTTTTTTACGTTGCGTCAGTTTACTCTgac                                                                                                                  >  1:124366/1‑147 (MQ=255)
                                     gcTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCagcagc                                                                                                             >  1:745635/1‑149 (MQ=255)
                                     gcTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGCGTTTTTACGTTGCTTCAGATTACTCTGCagcagc                                                                                                             >  1:440082/1‑149 (MQ=255)
                                         tGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGGGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCt                                                                                                         >  2:653849/1‑149 (MQ=255)
                                          gCCCACTGGTCAAGTAACCCCCTAAAGCCTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTg                                                                                                        <  1:658079/149‑1 (MQ=255)
                                              aCTGGTCCAGTAACCACAGAAATACTCCAGGGGATGAACAAAAAGAGTACCGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCttt                                                                                                    <  2:65632/149‑1 (MQ=255)
                                               cTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTtc                                                                                                   <  1:333537/149‑1 (MQ=255)
                                                   tCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTg                                                                                               <  1:716386/149‑1 (MQ=255)
                                                     aaGTCACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTACTGTAGATAAGTGGGGATACTACGGGTATAAAAACACCCGCCGGGGCGGGGTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTGTGCTTTCTCTGaa                                                                                             <  2:487494/149‑1 (MQ=255)
                                                     aaGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTGTGCAGCAGCTTCTGCTTTCTCTGaa                                                                                             <  1:447036/149‑1 (MQ=255)
                                                        tAACCACATAAATACTCCAGGGGATGGACAAAAAGAGTACAGCTATTGTAGATAGGTGGGGATACTACGGGTGTAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTGGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGa                                                                                          >  2:234835/1‑149 (MQ=255)
                                                        tAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGa                                                                                          >  1:233653/1‑149 (MQ=255)
                                                        catccAGATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTGTGAACGa                                                                                          >  2:496221/4‑149 (MQ=255)
                                                        caaCCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGa                                                                                          <  1:770362/148‑1 (MQ=255)
                                                                          aGGGGATGAACAAAAAGAGTACTGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTGTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTaa                                                                        <  2:457944/149‑1 (MQ=255)
                                                                          aGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAGCCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTaa                                                                        <  1:432129/149‑1 (MQ=255)
                                                                          aGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTaa                                                                        <  2:282253/149‑1 (MQ=255)
                                                                          aGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTaa                                                                        <  2:704172/149‑1 (MQ=255)
                                                                          aGGGGAGGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAACCCTGCCGGGGCTGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTaa                                                                        <  1:188322/149‑1 (MQ=255)
                                                                           ggggATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCCCTTTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAg                                                                       <  1:27064/149‑1 (MQ=255)
                                                                           ggggATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAg                                                                       >  2:266494/1‑149 (MQ=255)
                                                                                   aCAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTGTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAGGCCATCGgc                                                               >  1:722212/1‑149 (MQ=255)
                                                                                    cAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACAAAGAGTATGAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGgcg                                                              >  1:388477/1‑149 (MQ=255)
                                                                                     aaaaaGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTGTAAAAAAACCCGCCGGGGCGGGTTTTTTGACGTTGCGTCAGATTACTATGCACGAGCTTCTGCTTTCTGTGAAGGAACAACCAGCTCGATGGAAGCCAACGgcgc                                                             >  2:282483/1‑149 (MQ=255)
                                                                                     aaaaaGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGgcgc                                                             >  2:680365/1‑149 (MQ=255)
                                                                                     aaaaaGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGgcgc                                                             >  1:88486/1‑149 (MQ=255)
                                                                                     aaaaaGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGgcgc                                                             >  2:614966/1‑149 (MQ=255)
                                                                                     aaaaaGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGgcgc                                                             >  1:85311/1‑149 (MQ=255)
                                                                                     aaaaaGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGgcgc                                                             >  1:13788/1‑149 (MQ=255)
                                                                                      aaaaGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGCCGGGTTTTTTTACGTTGCTTCAGATCACTATGCAGCAGCATCTGCTTTCACTGAACGATCAACCAGCTCGATGTAAGCCATAGGcgcg                                                            >  1:560249/1‑149 (MQ=255)
                                                                                        aaGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGtt                                                          >  2:619625/1‑149 (MQ=255)
                                                                                        aaGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGtt                                                          <  1:436413/149‑1 (MQ=255)
                                                                                           agTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTc                                                       <  1:641644/149‑1 (MQ=255)
                                                                                                        tAGATAAGTGGGGATACTACGGGTATAGTGAAGCCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGGAAGCCATCGGCGCGTTGTCGCCTTCACGGCAg                                          >  2:7272/1‑149 (MQ=255)
                                                                                                        tAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGCAg                                          >  2:149530/1‑149 (MQ=255)
                                                                                                           aTAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCca                                       <  2:75989/149‑1 (MQ=255)
                                                                                                              aGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCacac                                    <  2:781807/149‑1 (MQ=255)
                                                                                                                tGGGGATACTACGGGTATAAAAACACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCGCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACtt                                  <  2:440082/149‑1 (MQ=255)
                                                                                                                tGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACtt                                  <  2:745635/149‑1 (MQ=255)
                                                                                                                tGGGGAGACTACGGGTATAAAAAAACCCGCCGGGGCGGGGTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCCCACtt                                  >  2:450120/1‑149 (MQ=255)
                                                                                                                  gggAGACTACGGGTATAAAAAACCCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCa                                <  2:301991/149‑1 (MQ=255)
                                                                                                                     agaatacGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGaa                             <  1:267896/145‑1 (MQ=255)
                                                                                                                       actacGGGTATAAACAAACCCGCCGGGGCGGGGTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCCCGGAAGCCACACTTCAGATTa                           >  1:428289/1‑149 (MQ=255)
                                                                                                                         tacGGGTATAAAAAAACCCGCCGGGACGGGTTGCACTTGGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGACATCGGCGCGTTGTCGCCTGCACGGCAGCCACACTTCAGAATACg                         >  1:194234/1‑149 (MQ=255)
                                                                                                                            gggTATAAAATAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGGCGCCTGCACGGAAGCCACACTTCAGAATACGAgt                      <  2:382035/149‑1 (MQ=255)
                                                                                                                              gTATAAAAAATCCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAgtgt                    <  1:275583/149‑1 (MQ=255)
                                                                                                                              gTATAAAAAACCCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAgtgt                    <  1:581190/149‑1 (MQ=255)
                                                                                                                                   aaaaaaCCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAacca               >  1:134406/1‑149 (MQ=255)
                                                                                                                                    aaaaaCCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAaccac              >  2:277780/1‑149 (MQ=255)
                                                                                                                                     aaaaCCCGCCCCGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGGACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGGCACACTTCAGAATACGGGTGTAgccacc             >  2:430159/1‑149 (MQ=255)
                                                                                                                                     aaaaCCCGCCCCGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAaccacc             >  2:582270/1‑149 (MQ=255)
                                                                                                                                      aaaCCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCCCGGAAGCCCCCCCTCAGAATACGAGTGTAACCACCg            >  1:320342/1‑149 (MQ=255)
                                                                                                                                        aCCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCGGc          >  2:533166/1‑149 (MQ=255)
                                                                                                                                                gggCGGGTTTTTTTACGATGCTTCAGATCACTCTGCAGCAGCTTCTGCTATCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCGGCACGGCTcg  <  2:207330/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
TCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCGGCACGGCTCG  >  NC_000913/3439443‑3439735

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