Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,439,586 G→C 7.2% intergenic (‑79/+28) 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,5860GC 7.2% 210.3 / 8.4 71intergenic (‑79/+28)yhdN/rplQDUF1992 family protein/50S ribosomal subunit protein L17
Reads supporting (aligned to +/- strand):  ref base G (37/28);  new base C (3/2);  total (40/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.50e-01

TCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCGGC  >  NC_000913/3439443‑3439727
                                                                                                                                               |                                                                                                                                             
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)
                                                                                                                                               |                                                                                                                                             
TCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCGGC  >  NC_000913/3439443‑3439727

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: