Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,489,858 G→A 5.9% A69V (GCA→GTA)  yjgR ← DUF853 family protein with NTPase fold

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,489,8580GA 5.9% 218.4 / 3.3 68A69V (GCA→GTA) yjgRDUF853 family protein with NTPase fold
Reads supporting (aligned to +/- strand):  ref base G (35/29);  new base A (2/2);  total (37/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.97e-01

AGATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGATATCCCACACCACCACCGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGGTTTTACCCGTCCCCGTTGCGCCGGTAATCAGCCCATGACGGTTA  >  NC_000913/4489720‑4489998
                                                                                                                                          |                                                                                                                                            
aGATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGATATCCCACACCACCACCGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGAc                                                                                                                                    >  1:241214/1‑149 (MQ=255)
      gAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGATATCCCACACCACCACCGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGt                                                                                                                              >  3:11952/1‑149 (MQ=255)
      gAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGATATCCCACACCACCACCGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGt                                                                                                                              <  3:655452/149‑1 (MQ=255)
        aaCCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGATATCCCACACCACCACCGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCa                                                                                                                            >  1:276266/1‑149 (MQ=255)
               gCCCGCACCGGATGGCCTTTCTCGCCAAAGATATCCCACACCACCACCGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGcc                                                                                                                     <  4:377359/149‑1 (MQ=255)
                 ccGCACCGGATGGCCTTTCTCGCCAAAGATATCCCACACCACCACCGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCtt                                                                                                                   >  4:30069/1‑149 (MQ=255)
                        ggATGGCCTTTCTCGCCAAAGATATCCCACACCACCACCGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTACCTGCGCGACGCCGGTCAGATCGCCTTTCACATc                                                                                                            >  2:542239/1‑149 (MQ=255)
                        ggATGGCCTTTCTCGCCAAAGATATCCCACACCACCACCGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTACCTGCGCGACGCCGGTCAGATCGCCTTTCACATc                                                                                                            >  2:413818/1‑149 (MQ=255)
                               cTTTCTCGCCAAAGATATCCCACACCACCACCGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATa                                                                                                     >  2:100311/1‑149 (MQ=255)
                                 ttCTCGCCAAAGATATCCCACACCACCACCGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATaaa                                                                                                   <  2:64468/149‑1 (MQ=255)
                                 ttCTCGCCAAAGATATCCCACACCACCACCGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATaaa                                                                                                   >  2:261395/1‑149 (MQ=255)
                                  tctcGCCAAAGATATCCCACACCACCACCGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAac                                                                                                  <  4:562471/149‑1 (MQ=255)
                                    tcGCCAAAGATATCCCACACCACCACCGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAacac                                                                                                >  2:503115/1‑149 (MQ=255)
                                     cGCCAAAGATATCCCACACCACCACCGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACAcc                                                                                               <  2:260490/149‑1 (MQ=255)
                                      gCCAAAGATATCCCACACCACCACCGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCg                                                                                              >  4:373535/1‑149 (MQ=255)
                                      gCCAAAGATATCCCACACCACCACCGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCg                                                                                              >  3:541365/1‑149 (MQ=255)
                                       ccAAAGATATCCCACACCACCACCGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCgg                                                                                             <  4:194092/149‑1 (MQ=255)
                                          aaGATATCCCACACCACCACCGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCAc                                                                                          <  3:337676/149‑1 (MQ=255)
                                              tatCCCACACCACCACCGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCg                                                                                      >  3:356179/1‑149 (MQ=255)
                                                 cccACACCACCACCGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGAtt                                                                                   >  4:230444/1‑149 (MQ=255)
                                                 cccACACCACCACCGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTCCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGAtt                                                                                   >  3:143677/1‑149 (MQ=255)
                                                  ccACACCACCACCGGCTTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGAttt                                                                                  >  3:656793/1‑149 (MQ=255)
                                                  ccACACCACCACCGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGAttt                                                                                  >  1:250401/1‑149 (MQ=255)
                                                    acacCACCACCGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTcc                                                                                >  4:160465/1‑149 (MQ=255)
                                                    acacCACCACCGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTcc                                                                                <  3:728764/149‑1 (MQ=255)
                                                    acacCACCACCGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTcc                                                                                <  1:261395/149‑1 (MQ=255)
                                                    acacCACCACCGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTACCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTcc                                                                                <  1:542239/149‑1 (MQ=255)
                                                    acacCACCACCGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTACCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTcc                                                                                <  1:413818/149‑1 (MQ=255)
                                                      accaccaccGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACCTCAGCCATAAACACCGGCACGCCGATTTCCGa                                                                              >  2:210718/1‑149 (MQ=255)
                                                      accaccaccGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTCCGa                                                                              >  3:34786/1‑149 (MQ=255)
                                                         accaccGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTCCGACaa                                                                           <  4:697120/149‑1 (MQ=255)
                                                           caccGGATTGGTATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTCCGACAATg                                                                         >  3:429600/1‑149 (MQ=255)
                                                           caccGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTCCGACAATg                                                                         >  3:734053/1‑149 (MQ=255)
                                                                     gCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAg                                                               >  2:130541/1‑149 (MQ=255)
                                                                         gCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGtttt                                                           <  2:276266/149‑1 (MQ=255)
                                                                             ttGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCa                                                       >  2:173733/1‑149 (MQ=255)
                                                                             ttGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCa                                                       <  1:513778/149‑1 (MQ=255)
                                                                                ccAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCg                                                    >  3:647939/1‑149 (MQ=255)
                                                                                  aGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTa                                                  >  3:213612/1‑149 (MQ=255)
                                                                                  aGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTa                                                  >  4:551281/1‑149 (MQ=255)
                                                                                  aGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTa                                                  >  4:127056/1‑149 (MQ=255)
                                                                                     cATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACg                                               <  1:100311/149‑1 (MQ=255)
                                                                                      aTTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACgg                                              >  1:362914/1‑149 (MQ=255)
                                                                                        tGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGGtt                                            <  4:138500/149‑1 (MQ=255)
                                                                                            gCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGAc                                                                                                                                    >  3:55138/1‑57 (MQ=255)
                                                                                            gCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGAc                                                                                                                                    <  4:55138/57‑1 (MQ=255)
                                                                                            gCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGGTTTTAc                                        <  3:127056/149‑1 (MQ=255)
                                                                                              cGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGGTTTTAccc                                      <  3:418185/149‑1 (MQ=255)
                                                                                              cGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGGTTTTAccc                                      >  3:51483/1‑149 (MQ=255)
                                                                                                 tatTTTCAACCCATGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCTGTCAGATCGCCTTTCACATCAGCCATTAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGGTTTCACCCGTc                                   <  2:422720/149‑1 (MQ=255)
                                                                                                       tAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGGTTTTACCCGTCCCCGtt                             <  4:212819/149‑1 (MQ=255)
                                                                                                       tAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGGTTTTACCCGTCCCCGtt                             >  3:219964/1‑149 (MQ=255)
                                                                                                        aaGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGGTTTTACCCGTCCCCGTTg                            <  2:678070/149‑1 (MQ=255)
                                                                                                        aaGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGGTTTTACCCGTCCCCGTTg                            <  3:230444/149‑1 (MQ=255)
                                                                                                           ccTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGATTTTACCCGTCCCCGTTgcgc                         <  4:734053/149‑1 (MQ=255)
                                                                                                             ttGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGGTTTTACCCGTCCCCGTTGCGCCg                       >  2:433955/1‑149 (MQ=255)
                                                                                                                 gAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGGTTTTACCCGTCCCCGTTGCGCCGGTaa                   <  2:693488/149‑1 (MQ=255)
                                                                                                                 gAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACAGCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGGTTTTACCCGTCCCCGTTGCGCCGGTaa                   <  3:551281/149‑1 (MQ=255)
                                                                                                                      gTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGAt                                                                                    <  3:612231/79‑1 (MQ=255)
                                                                                                                      gTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGAt                                                                                    >  4:612231/1‑79 (MQ=255)
                                                                                                                      gTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGGTTTTACCCGTCCCCGTTGCGCCGGTAATCAGc              >  4:326199/1‑149 (MQ=255)
                                                                                                                       tttttCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGGTTTTACCCGTCCCCGTTGCGCCGGTAATCAGcc             >  4:118098/1‑149 (MQ=255)
                                                                                                                          ttCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGGTTTTACCCGTCCCCGTTGCGCCGGTAATCAGCCCAt          <  4:693954/149‑1 (MQ=255)
                                                                                                                             cGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGGTTTTACCCGTCCCCGTTGCGCCGGTAATCAGCCCATGAc       <  2:233381/149‑1 (MQ=255)
                                                                                                                                cacCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGGTTTTACCCGTCCCCGTTGCGCCGGTAATCAGCCCATGACGGt    <  2:425229/149‑1 (MQ=255)
                                                                                                                                 acCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGGTTTTACCCGTCCCCGTTGCGCCGGTAATCAGCCCATGACGGtt   >  1:743016/1‑149 (MQ=255)
                                                                                                                                 acCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGGTTTTACCCGTCCCCGTTGCGCCGGTAATCAGCCCATGACGGtt   >  3:571284/1‑149 (MQ=255)
                                                                                                                                  ccGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGGTTTTACCCGTCCCCGTTGCGCCGGTAATCAGCCCATGACGGTTa  <  4:647939/149‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                            
AGATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGATATCCCACACCACCACCGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGGTTTTACCCGTCCCCGTTGCGCCGGTAATCAGCCCATGACGGTTA  >  NC_000913/4489720‑4489998

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑