Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,533,148 C→T 5.7% L33L (CTG→TTG)  pptA → 4‑oxalocrotonate tautomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,533,1480CT 5.7% 222.5 / 2.4 70L33L (CTG→TTG) pptA4‑oxalocrotonate tautomerase
Reads supporting (aligned to +/- strand):  ref base C (32/34);  new base T (2/2);  total (34/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.35e-01

ACAATGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGCGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTCATCTGAACAGTAAAGACAGTTCGATAAGCATTGCTCTACAGCAGATTCAACCAGAATCTTGGCAAGCTATCTGGGATGCCGAAATCGCGCCCCAAATGGAGGCTTTGATAAAGAAACCTGGTTATAGCATGAATGCTTAACCCGCC  >  NC_000913/1533006‑1533291
                                                                                                                                              |                                                                                                                                               
ccaATGGGTGTTAATTTTACTGACATCTTTCAAAAGGAGCGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTCATCTGAAc                                                                                                                                            <  2:313438/147‑1 (MQ=255)
acaATGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGCGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTCATCTGaa                                                                                                                                             >  4:232500/1‑147 (MQ=255)
 caATGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGCGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTCATCTGAACAg                                                                                                                                          <  1:149319/149‑1 (MQ=255)
    tgtgtgTTAATTTTACTGACATCTTTCAAAAGGAGCGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTCATCTGAACAGTa                                                                                                                                        >  7:152704/1‑148 (MQ=255)
      tgtgtTAATTTTACTGACATCTTTCAAAAGGAGCGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTCATCTGAACAGTaaa                                                                                                                                      <  3:297412/148‑1 (MQ=255)
        tgtTAATTTTACTGACATCTTTCAAAAGGAGCGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTCATCTGAACAGTa                                                                                                                                        <  2:361903/144‑1 (MQ=255)
        tgtTAATTTTACTGACATCTTTCAAAAGGAGCGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTCATCTGAACAGTa                                                                                                                                        >  1:361903/1‑144 (MQ=255)
        tgtTAATTTTACTGACATCTTTCAAAAGGAGCGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTCATCTGAACAGTAAAGa                                                                                                                                    >  2:362864/1‑148 (MQ=255)
          ttAATTTTACTGACATCTTTCAAAAGGAGCGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTCATCTGAACAGTAAAGACAg                                                                                                                                 >  5:169841/1‑149 (MQ=255)
            aaTTTTACTGACATCTTTCAAAAGGAGCGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTCATTTGAACAGTAAAGACAGtt                                                                                                                               >  3:174426/1‑149 (MQ=255)
             aTTTTACTGACATCTTTCAAAAGGAGCGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTCATCTGAACAGTAAAGACAGTTc                                                                                                                              >  5:157475/1‑149 (MQ=255)
             aTTTTACTGACATCTTTCAAAAGGAGCGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTCATCTGAACAGTAAAGACAGTTc                                                                                                                              <  3:223372/149‑1 (MQ=255)
              ttttACTGACATCTTTCAAAAGGAGCGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTCATCTGAACAGTAAAGACAGTTc                                                                                                                              >  5:279743/1‑148 (MQ=255)
                 tACTGACATCTTTCAAAAGGAGCGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTCATCTGAACAGTAAAGACAGTTCGAt                                                                                                                           <  6:169841/148‑1 (MQ=255)
                   cTGACATCTTTCAAAAGGAGCGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTCATCTGAACAGTAAAGACAGTTCGATaa                                                                                                                         <  5:279949/148‑1 (MQ=255)
                        aTCTTTCAAAAGGAGCGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTCATCTGAACAGTAAAGACAGTTCGATAAGCATTg                                                                                                                   >  8:47442/1‑149 (MQ=255)
                            ttCAAAAGGAGCGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTCATCTGAACAGTAAAGACAGTTCGATAAGCATTGctc                                                                                                                >  1:207814/1‑148 (MQ=255)
                             tCAAAAGGAGCGTATTCATGGCGGACATCGAGATTAAATGTTTTCGTCGTGAACTGGAAGAACAACAAAAAGCAGCACTTGCTGCAGATATTTCCGACGTTATTATTCGTCATCTGAACAGTAAAGAGAGTTCGATAAGCATTGct                                                                                                                 <  8:118221/146‑1 (MQ=255)
                               aaaaGGAGCGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTCATCTGAACAGTAAAGACAGTTCGATAAGCATTGCTCTa                                                                                                              >  3:154879/1‑147 (MQ=255)
                               aaaaGGAGCGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTCATCTGAACAGTAAAGACAGTTCGATAAGCATTGCTCTAca                                                                                                            <  7:243958/149‑1 (MQ=255)
                                  aGGAGCGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTTAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTCATCTGAACAGTAAAGACAGTTCGATAAGCATTGCTCTAcagca                                                                                                         <  4:21287/149‑1 (MQ=255)
                                        gTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTCATCTGAACAGTAAAGACAGTTCGATAAGCATTGCTCTACAGCAGATTc                                                                                                    >  8:302012/1‑148 (MQ=255)
                                        gTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTCATCTGAACAGTAAAGACAGTTCGATAAGCATTGCTCTACAGCAGATTc                                                                                                    >  7:119357/1‑148 (MQ=255)
                                           aTCATGCCGCTCATAGACATTAAATGTTTTCGGGTTTAACTGGACTTACAACAAATAGCAGCACTTGCTGCTGATATTACCGACGTTATTATTCGTCATCTGAACAGTAAAGACAGTTCGATAAGCATTGCTCTACAGCAGATTCAAc                                                                                                 <  6:35403/148‑1 (MQ=255)
                                           aTCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAATCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTCATCTGAACAGTAAAGACAGTTCGATAAGCATTGCTCTACAGCAGATTCAAc                                                                                                 <  4:100509/148‑1 (MQ=255)
                                                     acaTCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTCATCTGAACAGTAAAGACAGTTCGATAAGCATTGCTCTACAGCAGATTCAACCAGAATCTTgg                                                                                      >  2:224251/1‑149 (MQ=255)
                                                      caTCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTCATCTGAACAGTAAAGACAGTTCGATAAGCATTGCTCTACAGCAGATTCAACCAGAATCTTTGc                                                                                     <  8:19608/149‑1 (MQ=255)
                                                          gACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTAGTCGTCCTCTGAACAGTAAAGACAGTTCGATAAGCATTGCTCTACAGCAGATTCAACCAGAATCTTGGCAAGc                                                                                 >  5:265723/1‑149 (MQ=255)
                                                            cATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTCATCTGAACAGTAAAGACAGTTCGATAAGCATTGCTCTACAGCAGATTCAACCAGAATCTTGGCAAGCTa                                                                               >  4:279445/1‑149 (MQ=255)
                                                             aTTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTCATCTGAACAGTAAAGACAGTTCGATAAGCATTGCTCTACAGCAGATTCAACCAGAATCTTGGCAAGCTAt                                                                              >  2:90749/1‑149 (MQ=255)
                                                                aaaTGTTTTCCGCGTGAACTGGACGAAAAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTCATCTGAACAGTAAAGACAGTTCGATAAGCATTGCTCTACAGCAGATTCAACCAGAATCTTGGCAAGCTATCTg                                                                           >  7:325260/1‑149 (MQ=255)
                                                                 aaTGTTTTCCGCGTTAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTCATCTGAACAGTAAAGACAGTTCGATAAGCATTGCTCTACAGCAGATTCAACCAGAATCTTGGCAAGCTATCTg                                                                           <  8:152704/148‑1 (MQ=255)
                                                                 aaTGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTCATCTGAACAGTAAAGACAGTTCGATAAGCATTGCTCTACAGCAGATTCAACCAGAATCTTGGCAAGCTATCTg                                                                           <  6:157475/148‑1 (MQ=255)
                                                                           gcgTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTCATCTGAACAGTAAAGACAGTTCGATAAGCATTGCTCTACAGCAGATTCAACCAGAATCTTGGCAAGCTATCTGGGATGCCGaaa                                                                <  5:345165/149‑1 (MQ=255)
                                                                                aaCTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTCATCTGAACAGTAAAGACAGTTCGATAAGCATTGCTCTACAGCAGATTCAACCAGAATCTTGGCAAGCTATCTGGGATGCCGAAATcgcg                                                           <  4:154879/149‑1 (MQ=255)
                                                                                aaCTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTCATCTGAACAGTAAAGACAGTTCGATAAGCATTGCTCTACAGCAGATTCAACCAGAATCTTGGCAAGCTATCTGGGATGCCGAAATcgc                                                            <  3:232500/148‑1 (MQ=255)
                                                                                     gACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTCATCTGAACAGTAAAGACAGTACGAGAAGCATTGCTCTACAGCAGATTCAACCAGAATCTTGGCAAGCTATCTGGGATGCCGAAATCGCGcccc                                                       >  5:53340/1‑148 (MQ=255)
                                                                                      aCGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTCATCTGAACAGTAAAGACAGTTCGATAAGCATTGCTCTACAGCAGATTCAACCAGAATCTTGGCAAGCTATCTGGGATGCCGAAATCGCGCCCCa                                                      >  6:203001/1‑148 (MQ=255)
                                                                                      aCGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTCATCTGAACAGTAAAGACAGTTCGATAAGCATTGCTCTACAGCAGATTCAACCAGAATCTTGGCAAGCTATCTGGGATGCCGAAATCGCGCCCCa                                                      >  5:229802/1‑148 (MQ=255)
                                                                                        gAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTCATCTGAACAGTAAAGACAGTTCGATAAGCATTGCTCTACAGCAGATTCAACCAGAATCTTGGCAAGCTATCTGGGATGCCGAAATCGCGCCCCAAAt                                                   >  7:130202/1‑149 (MQ=255)
                                                                                            aacaaAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTCATCTGAACAGTAAAGACAGTTCGATAAGCATTGCTCTACAGCAGATTCAACCAGAATCTTGGCAAGCTATCTGGGATGCCGAAATCGCGCCCCAAATGGAg                                               <  3:362098/149‑1 (MQ=255)
                                                                                              caaAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTCATCTGAACAGTAAAGACAGTTCGATAAGCATTGCTCTACAGCAGATTCAACCAGAATCTTGGCAAGCTATCTGGGATGCCGAAATCGCGACCCAAATGGAGGc                                             >  4:371035/1‑149 (MQ=255)
                                                                                              caaAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTCATCTGAACAGTAAAGACAGTTCGATAAGCATTGCTCTACAGCAGATTAAACCAGAATCTTGGCAAGCTATCTGGGATGCCGAAATCGCGCCCCAAATGGAgga                                             >  6:366006/1‑148 (MQ=255)
                                                                                               aaaaaGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTCATTTGAACAGTAAAGACAGTTCGATAAGCATTGCTCTACAGCAGATTCAACCAGAATCTTGGCAAGCTATCTGGGATGCCGAAATCGCGCCCCAAATGGAGGCt                                            >  6:123076/1‑149 (MQ=255)
                                                                                                 aaaGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTCATCTGAACAGTAAAGACAGTTCGATAAGCATTGCTCTACAGCAGATTCAACCAGAATCTTGGCAAGCTATCTGGGATGCCGAAATCGCGCCCCAAATGGAGGCttt                                          <  1:362864/149‑1 (MQ=255)
                                                                                                  aaGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTCATCTGAACAGTAAAGACAGTTCGATAAGCATTGCTCTACAGCAGATTCAACCAGAATCTTGGCAAGCTATCTGGGATGCCGAAATCGCGCCCCAAATGGAGGCTTTg                                         <  7:302012/149‑1 (MQ=255)
                                                                                                   agcagcACTTGCTGCAGATATTACCGACGTTATTATTCGTCATCTGAACAGTAAAGACAGTTCTATAAGCATTGCTCTACAGCATATTCAACCAGAATCTTGGCAAGCTATCTGGGATGCCGAAATCGCGCCCCAAATGGAGGCTTTGa                                        <  2:328123/149‑1 (MQ=255)
                                                                                                        caCTTGCTGCAGATATTACCGACGTTATTATTCGTCATCTGAACAGTAAAGACAGTTCGATAAGCATTGCTCTACAGCAGATTCAACCAGAATCTTGGCAAGCTATCTGGGATGCCGAAATCGCGCCCCAAATGGAGGCTTTGATa                                      >  5:9021/1‑146 (MQ=255)
                                                                                                              ctgcAGATATTACCGACGTTATTATTCGTCATTTGAACAGTAAAGACAGTTCGATAAGCATTGCTCTACAGCAGATTCAACCAGAATCTTGGCAAGCTATCTGGGATGCCGAAATCGCGCCCCAAATGGAGGCTTTGATAAAGAAAcc                              <  4:174426/148‑1 (MQ=255)
                                                                                                              ctgcAGATATTACCGACGTTATTATTCGTCATCTGAACAGTAAAGACAGTTCGATAAGCATTGCTCTACAGCAGATTCAACCAGAATCTTGGCAAGCTATCTGGGATGCCGAAATCGCGCCCCAAATGGAGGCTTTGATAAAGAAACCt                             <  8:130202/149‑1 (MQ=255)
                                                                                                                 cAGATATTACCGACGTTATTATTCGTCATCTGAACAGTAAAGACCGTTCGATAAGCAATGCTCTACAGCAGATTCAACCAGAATCTTGGCAAGCTATCTGGGATGCCGAAATCGCGCCCCAAATGGAGGCTTTGATAAAGAAACCt                             <  1:20938/146‑1 (MQ=255)
                                                                                                                 cAGATATTACCGACGTTATTATTCGTCATCTGAACAGTAAAGACAGTTCGATAAGCATTGCTCTACAGCAGATTCAACCAGAATCTTGGCAAGCTATCTGGGATGCCGAAATCGCGCCCCAAATGGAGGCTTTGATAAAGAAACCTgg                           <  7:47442/148‑1 (MQ=255)
                                                                                                                    atatTACCGACGTTATTATTCGTCATCTGAACAGTAAAGACAGTTCGATAAGCATTGCTCTACAGCAGATTCAACCAGAATCTTGGCAAGCTATCTGGGATGCCGAAATCGCGCCCCAAATGGAGGCTTTGATaaa                                    <  8:372593/136‑1 (MQ=255)
                                                                                                                    atatTACCGACGTTATTATTCGTCATCTGAACAGTAAAGACAGTTCGATAAGCATTGCTCTACAGCAGATTCAACCAGAATCTTGGCAAGCTATCTGGGATGCCGAAATCGCGCCCCAAATGGAGGCTTTGATaaa                                    >  7:372593/1‑136 (MQ=255)
                                                                                                                    atatTACCGACGTTATTATTCGTCATCTGAACAGTAAAGACAGTTCGATAAGCATTGCTCTACAGCAGATTCAACCAGAATCTTGGCAAGCTATCTGGGATGCCGAAATCGCGCCCCAAATGGAGGCTTTGATAAAGAAACCTGGTtat                       <  1:90749/149‑1 (MQ=255)
                                                                                                                     tatTACCGACGTTATTATTCGTCATCTGAACAGTAAAGACAGTTCGATAAGCATTGCTCTACAGCAGATTCAACCAGAATCTTGGCAAGCTATCTGGGATGCCGAAATCGCGCCCCAAATGGAGGCTTTGATAAAGAAACCTGGTtata                      <  5:203001/149‑1 (MQ=255)
                                                                                                                     tatTACCGACGTTATTATTCGTCATCTGAACAGTAAAGACAGTTCGATAAGCATTGCTCTACAGCAGATTCAACCAGAATCTTGGCAAGCTATCTGGGATGCCGAAATCGCGCCCCAAATGGAGGCTTTGATAAAGAAACCTGGTtata                      <  2:207814/149‑1 (MQ=255)
                                                                                                                      atTACCGACGTTATTATTCGTCATCTGAACAGTAAAGACAGTTCGATAAGCATTGCTCTACAGCAGATTCAACCAGAATCTTGGCAAGCTATCTGGGATGCCGAAATCGCGCCCCAAATGGAGGCTTTGa                                        >  2:58190/1‑130 (MQ=255)
                                                                                                                      atTACCGACGTTATTATTCGTCATCTGAACAGTAAAGACAGTTCGATAAGCATTGCTCTACAGCAGATTCAACCAGAATCTTGGCAAGCTATCTGGGATGCCGAAATCGCGCCCCAAATGGAGGCTTTGa                                        <  1:58190/130‑1 (MQ=255)
                                                                                                                       ttACCGACGTTATTATTCGTCATCTGAACAGTAAAGACAGTTCGATAAGCATTGCTCTACAGCAGATTCAACCAGAATCTTGGCAAGCTATCTGGGATGCCGAAATCGCGCCCCAAATGGAGGCTTTGATAAAGAAACCTGGTTATAg                     >  2:216668/1‑148 (MQ=255)
                                                                                                                       ttACCGACGTTATTATTCGTCATCTGAACAGTAAAGACAGTTAGATAAGCATTGCTCTACAGCAGATTCAACCAGAATCTTGGCAAGCTATCTGGGATGCCGAAATCGCGCCCCAAATGGAGGCTTTGATAAAGAAACCTGGTTATAGc                    <  6:279743/149‑1 (MQ=255)
                                                                                                                            gACGTTATTATTCGTCATCTGAACAGTAAAGACAGTTCGATAAGCATTGCTCTACAGCAGATTCAACCAGAATCTTGGCAAGCTATCTGGGATGCCGAAATCGCGCCCCAAATGGAGGCGTTGATATAGAATCCTGGTTATAGCATGaa               <  3:371035/149‑1 (MQ=255)
                                                                                                                             aCGTTATTATTCGTCATCTGAACAGTAAAGACAGTTCGATAAGCATTGCTCTACAGCAGATTCAACCAGAATCTTGGCAAGCTATCTGGGATGCCGAAATCGCGCCCCAAATGGAGGCTTTGATAAAGAAACCTGGTTATAGCATGAAt              >  4:341677/1‑149 (MQ=255)
                                                                                                                             aCGTTATTATTCGTCATCTGAACAGTAAAGACAGTTCGATAAGCATTGCTCTACAGCAGATTCAACCAGAATCTTGGCAAGCTATCTGGGATGCCGAAATCGCGCCCCAAATGGAGGCTTTGATAAAGAAACCTGGTTATAGCATGAAt              >  4:347159/1‑149 (MQ=255)
                                                                                                                                      ttCGTCATCTGAACAGTAAAGACAGTTCGATAAGCATTGCTCTACAGCAGATTCAACCAGAATCTTGGCAAGCTATCTGGGATGCCGAAATCGCGCCCCAAATGGAGGCTTTGATAAAGAAACCTGGTTATAGCATGAATGGTTAAccc     >  2:232277/1‑149 (MQ=255)
                                                                                                                                         gTCATTTGAACAGTAAAGACAGTTCGATAAGCATTGCTCTACAGCAGATTCAACCAGAATCTTGGCAAGCTATCTGGGATGCCGAAATCGCGCCCCAAATGGAGGCTTTGATAAAGAAACCTGGTTATAGCATGAATGCTTAACccgcc  <  5:123076/149‑1 (MQ=255)
                                                                                                                                         gTCATCTGAACAGTAAAGACAGTTCGATAAGCATTGCTCTACAGCAGATTCAACCAGAATCTTGGCAAGCTATCTGGGATGCCGAAATCGCGCCCCAAATGGAGGCTTTGATAAAGAAACCTGGTTATAGCATGAATGCTTAACccgcc  >  6:255507/1‑149 (MQ=255)
                                                                                                                                         gTCATCTGAACAGTAAAGACAGTTCGATAAGCATTGCTCTACAGCAGATTCAACCAGAATCTTGGCAAGCTATCTGGGATGCCGAAATCGCGCCCCAAATGGAGGCTTTGATAAAGAAACCTGGTTATAGCATGAATGCTTAACccgc   <  6:9021/148‑1 (MQ=255)
                                                                                                                                          tCATCTGAACAGTAAAGACAGTTCGATAAGCATTGCTCTACAGCAGATTCAACCAGAATCTTGGCAAGCTATCTGGGATGCCGAAATCGCGCCCCAAATGGAGGCtt                                           >  2:267098/1‑107 (MQ=255)
                                                                                                                                          tCATCTGAACAGTAAAGACAGTTCGATAAGCATTGCTCTACAGCAGATTCAACCAGAATCTTGGCAAGCTATCTGGGATGCCGAAATCGCGCCCCAAATGGAGGCtt                                           <  1:267098/107‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                               
ACAATGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGCGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTCATCTGAACAGTAAAGACAGTTCGATAAGCATTGCTCTACAGCAGATTCAACCAGAATCTTGGCAAGCTATCTGGGATGCCGAAATCGCGCCCCAAATGGAGGCTTTGATAAAGAAACCTGGTTATAGCATGAATGCTTAACCCGCC  >  NC_000913/1533006‑1533291

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