Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,803,227 A→G 6.2% N885N (AAT→AAC PP_3361 ← aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,803,2270AG 6.2% 605.2 / 18.7 193N885N (AAT→AACPP_3361aminotransferase
Reads supporting (aligned to +/- strand):  ref base A (89/92);  new base G (3/9);  total (92/101)
Fisher's exact test for biased strand distribution p-value = 1.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

Only 100 of 201 total aligned reads displayed.

CCAATCACGTCGAATTGCTGCTGAAGCTTGCGCAAACCCTCGAGCAGGTAATGGCCGACACTCAAGGCGTTTTCTTTCAACTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGCATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCC  >  NC_002947/3803078‑3803376
                                                                                                                                                     |                                                                                                                                                      
ccAATCACGTCGAATTGCTGCTGAAGCTTGCGCAAACCCTCGAGCAGGTAATGGCCGACACTCAAGGCGTTTTCTTTCAACTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTa                                                                                                                                                        <  1:1227468/150‑1 (MQ=255)
ccAATCACGTCGAATTGCTGCTGAAGCTTGCGCAAACCCTCGAGCAGGTAATGGCCGACACTCAAGGCGTTTTCTTTCAACTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTa                                                                                                                                                        <  1:3811331/150‑1 (MQ=255)
ccAATCACGTCGAATTGCTGCTGAAGCTTGCGCAAACCCTCGAGCAGGTAATGGCCGACACTCAAGGCGTTTTCTTTCAACTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTa                                                                                                                                                        >  2:29577/1‑150 (MQ=255)
ccAATCACGTCGAATTGCTGCTGAAGCTTGCGCAAACCCTCGAGCAGGTAATGGCCGACACTCAAGGCGTTTTCTTTCAACTGGTCGCGCTCGAATGCATAGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTa                                                                                                                                                        >  1:2199200/1‑150 (MQ=255)
  aaTCACGTCGAATTGCTGCTGAAGCTTGCGCAAACCCTCGAGCAGGTAATGGCCGACACTCAAGGCGTTTTCTTTCAACTGGTCGCGCTCGATTGCATCGAGTAACGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTAtt                                                                                                                                                      <  2:2860122/150‑1 (MQ=255)
   aTCACGTCGAATTGCTGCTGAAGCTTGCGCAAACCCTCGAGCAGGTAATGGCCGACACTCACGGCGTTTTCTTTCAACTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTg                                                                                                                                                     <  2:3371288/150‑1 (MQ=255)
    tCACGTCGAATTGCTGCTGAAGCTTGCGCAAACCCTCGAGCAGGTAATGGCCGACACTCAAGGCGTTTTCTTTCAACTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGa                                                                                                                                                    <  2:1960808/150‑1 (MQ=255)
      aCGTCGAATTGCTGCTGAAGCTTGCGCAAACCCTCGAGCAGGTAATGGCCGACACTCAAGGCGTTTTCTTTCAACTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGaaa                                                                                                                                                  >  1:3415916/1‑150 (MQ=255)
      aCGTCGAATTGCTGCTGAAGCTTGCGCAAACCCTCGAGCAGGTAATGGCCGACACTCAAGGCGTTTTCTTTCAACTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGaaa                                                                                                                                                  <  2:2029857/150‑1 (MQ=255)
        gTCGAATTGCTGCTGAAGCTTGCGCAAACCCTCGAGCAGGTAATGGCCGACACTCAAGGCGTTTTCTTTCAACTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATa                                                                                                                                                <  1:1792724/150‑1 (MQ=255)
            aaTTGCTGCTGAAGCTTGCGCAAACCCTCGAGCAGGTAATGGCCGACACTCAAGGCGTTTTCTTTCAACTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTcc                                                                                                                                            >  2:4924214/1‑150 (MQ=255)
             aTTGCTGCTGAAGCTTGCGCAAACCCTCGAGCAGGTAATGGCCGACACTCAAGGCGTTTTCTTTCAACTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCa                                                                                                                                           <  1:2669229/150‑1 (MQ=255)
              ttGCTGCTGAAGCTTGCGCAAACCCTCGAGCAGGTAATGGCCGACACTCAAGGCGTTTTCTTTCAACTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCAt                                                                                                                                          <  2:1410784/150‑1 (MQ=255)
                   gctgAAGCTTGCGCAAACCCTCGAGCAGGTAATGGCCGACACTCAAGGCGTTTTCTTTCAACTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCgt                                                                                                                                     <  1:4645093/150‑1 (MQ=255)
                     tgAAGCTTGCGCAAACCCTCGAGCAGGTAATGGCCGACACTCAAGGCGTTTTCTTTCAACTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCgttg                                                                                                                                   >  2:2677763/1‑150 (MQ=255)
                     tgAAGCTTGCACAAACCCTCGAGCAGGTAATGGCCGACACTCAAGGCGTTTTCTTTCAACTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTg                                                                                                                                                     <  1:2033339/132‑1 (MQ=255)
                      gAAGCTTGCGCAAACCCTCGAGCAGGTAATGGCCGACACTCAAGGCGTTTTCTTTCAACTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCgttgt                                                                                                                                  >  1:3395973/1‑150 (MQ=255)
                       aaGCTTGCGCAAACCCTCGAGCAGGTAATGGCCGACACTCAAGGCGTTTTCTTTCAACTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCgttgtt                                                                                                                                 >  2:3477964/1‑150 (MQ=255)
                       aaGCTTGCGCAAACCCTCGAGCAGGTAATGGCCGACACTCAAGGCGTTTTCTTTCAACTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCgttgtt                                                                                                                                 <  1:4266018/150‑1 (MQ=255)
                        aGCTTGCGCAAACCCTCGAGCAGGTAATGGCCGACACTCAAGGCGTTTTCTTTCAACTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGttgttg                                                                                                                                <  1:830300/150‑1 (MQ=255)
                        aGCTTGCGCAAACCCTCGAGCAGGTAATGGCCGACACTCAAGGCGTTTTCTTTCAACTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGttgttg                                                                                                                                <  1:830089/150‑1 (MQ=255)
                         gCTTGCGCAAACCCTCGAGCAGGTAATGGCCGACACTCAAGGCGTTTTCTTTCAACTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGa                                                                                                                               >  2:2102513/1‑150 (MQ=255)
                           ttGCGCAAACCCTCGAGCAGGTAATGGCCGACACTCAAGGCGTTTTCTTTCAACTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGaaa                                                                                                                             <  1:3241174/150‑1 (MQ=255)
                           ttGCGCAAACCCTCGAGCAGGTAATGGCCGACACTCAAGGCGTTTTCTTTCAACTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGaaa                                                                                                                             <  1:2990674/150‑1 (MQ=255)
                                 aaaCCCTCGAGCAGGTAATGGCCGACACTCAAGGCGTTTTCTTTCAACTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGCATCg                                                                                                                      <  1:1963925/150‑1 (MQ=255)
                                 aaaCCCTCGAGCAGGTAATGGCCGACACTCAAGGCGTTTTATTTCAACTGGTCGCGCTAGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGTTATTGAAATACTCCATGCCGTTGTTGAAAGCATCg                                                                                                                      <  2:3573730/150‑1 (MQ=255)
                                  aaCCCTCGAGCAGGTAATGGCCGACACTCAAGGCGTTTTCTTTCAACTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGCATCgg                                                                                                                     >  2:3268792/1‑150 (MQ=255)
                                   aCCCTCGAGCAGGTAATGGCCGACACTCAAGGCGTTTTCTTTCAACTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGCATCGGc                                                                                                                    >  1:2448392/1‑150 (MQ=255)
                                     ccTCGAGCAGGTAATGGCCGACACTCAAGGCGTTTTCTTTCAACTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGCATCGGCaa                                                                                                                  <  1:3742058/150‑1 (MQ=255)
                                     ccTCGAGCAGGTAATGGCCGACACTCAAGGCGTTTTCTTTCAACTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGCATCGGCaa                                                                                                                  >  2:1438711/1‑150 (MQ=255)
                                     ccTCGAGCAGGTAATGGCCGACACTCAAGGCGTTTTCTTTCAACTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGCATCGGCaa                                                                                                                  >  2:1748551/1‑150 (MQ=255)
                                     ccTCGAGCAGGTAATGGCCGACACTCAAGGCGTTTTCTTTCAACTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGCATCGGCaa                                                                                                                  <  2:89422/150‑1 (MQ=255)
                                      cTCGAGCAGGTAATGGCCGACACTCAAGGCGTTTTCTTTCAACTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGCATCGGCAAc                                                                                                                 >  2:3601631/1‑150 (MQ=255)
                                           gCAGGTAATGGCCGACACTCAAGGCGTTTTCTTTCAACTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGCATCGGCAACCTcgc                                                                                                            >  1:162416/1‑150 (MQ=255)
                                           gCAGGTAATGGCCGACACTCAAGGCGTTTTCTTTCAACTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCACGCCGTTGTTGAAAGCATCGGCAACCTcgc                                                                                                            >  2:800229/1‑150 (MQ=255)
                                             aGGTAATGGCCGACACTCAAGGCGTTTTCTTTCAACTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGCATCGGCAACCTcgcgc                                                                                                          >  2:4774573/1‑150 (MQ=255)
                                               gTAATGGCCGACACTCAAGGCGTTTTCTTTCAACTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGCATCGGCAACCTCGCGCGt                                                                                                        >  1:5078056/1‑150 (MQ=255)
                                               gTAATGGCCGACACTCAAGGCGTTTTCTTTCAACTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGAATTGAAATACTCCATGCCGTAGTTGAAAGCATCGGCAACCTCGCGCGt                                                                                                        >  1:5078205/1‑150 (MQ=255)
                                                tAATGGCCGACACTCAAGGCGTTTTCTTTCAACTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGCATCGGCAACCTCGCGCGTg                                                                                                       <  1:1006918/150‑1 (MQ=255)
                                                 aaTGGCCGACACTCAAGGCGTTTTCTTTCAACTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGCATCGGCAACCTCGCGCGTgg                                                                                                      <  2:4058533/150‑1 (MQ=255)
                                                   tGGCCGACACTCAAGGCGTTTTCTTTCAACTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGCATCGGCAACCTCGCGCGTGGTc                                                                                                    <  2:5022027/150‑1 (MQ=255)
                                                     gCCGACACTCAAGGCGTTTTCTTTCAACTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGCATCGGCAACCTCGCGCGTGGTCAc                                                                                                  >  2:1248112/1‑150 (MQ=255)
                                                     gCCGACACTCAAGGCGTTTTCTTTCAACTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGCATCGGCAACCTCGCGCGTGGTCAc                                                                                                  >  2:4167095/1‑150 (MQ=255)
                                                       cGACACTCAAGGCGTTTTCTTTCAACTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGCATCGGCAACCTCGCGCGTGGTCACTa                                                                                                <  2:614665/150‑1 (MQ=255)
                                                       cGACACTCAAGGCGTTTTCTTTCAACTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGCATCGGCAACCTCGCGCGTGGTCACTa                                                                                                >  1:1067629/1‑150 (MQ=255)
                                                            cTCAAGGCGTTTTCTTTCAACTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGCATCGGCAACCTCGCGCGTGGTCACTAccgcc                                                                                           <  1:1167891/150‑1 (MQ=255)
                                                             tCAAGGCGTTTTCTTTCAACTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGCATCGGCAACCTCGCGCGTGGTCACTAccgccg                                                                                          >  1:5007019/1‑150 (MQ=255)
                                                              cAAGGCGTTTTCTTTCAACTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGCATCGGCAACCTCGCGCGTGGTCACTACCGCCGa                                                                                         <  1:790962/150‑1 (MQ=255)
                                                               aaGGCGTTTTCTTTCAACTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGCATCGGCAACCTCGCGCGTGGTCACTACCGCCGAc                                                                                        >  1:303239/1‑150 (MQ=255)
                                                                    gTTTTCTTTCAACTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGCATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCgg                                                                                   >  2:1779560/1‑150 (MQ=255)
                                                                    gTTTTCTTTCAACTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGCATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCgg                                                                                   <  1:3965670/150‑1 (MQ=255)
                                                                     ttttCTTTCAACTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGCATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCggg                                                                                  <  2:3856406/150‑1 (MQ=255)
                                                                     ttttCTTTCAACTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGCATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCggg                                                                                  >  2:5095134/1‑150 (MQ=255)
                                                                       ttctttcAACTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGCATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTg                                                                                >  2:425489/1‑150 (MQ=255)
                                                                        tctttcAACTGGTCGCGCTCGAGTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGCATCGGCAACCTCGCGCGTGGTCACTACCGCCGACCTCGGGTgg                                                                               <  2:4781851/150‑1 (MQ=255)
                                                                            tcAACTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGCATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGt                                                                           >  2:1700108/1‑150 (MQ=255)
                                                                                cTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGCATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTAcc                                                                       >  1:2955152/1‑150 (MQ=255)
                                                                                  ggTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGCATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCaa                                                                     <  2:4685248/150‑1 (MQ=255)
                                                                                    tCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGCATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATg                                                                   >  1:1569631/1‑150 (MQ=255)
                                                                                      gcgcTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGCATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATggg                                                                 >  2:4579391/1‑150 (MQ=255)
                                                                                      gcgcTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGCATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATggg                                                                 <  1:2876843/150‑1 (MQ=255)
                                                                                         cTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGCATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGttt                                                              <  2:4694425/150‑1 (MQ=255)
                                                                                          tCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGCATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTg                                                             >  1:2766207/1‑150 (MQ=255)
                                                                                          tCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGCATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTg                                                             >  1:4782619/1‑150 (MQ=255)
                                                                                          tCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGCATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTg                                                             >  1:2576147/1‑150 (MQ=255)
                                                                                           cgatcacATTGAGCACTGCGAGGCCTGCGGCACAGGAAACCGGATTGCCACCGAAGGTGTTGAAATATTCCATGCCATTGTTAAAGCTGTCTGCCACTTCCCGCGTTGTGACAACCGCAGACATTGGATGACCGTTACCAATGGGTTTGc                                                            <  1:3069233/144‑1 (MQ=255)
                                                                                           cGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGCATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGc                                                            <  1:2701685/150‑1 (MQ=255)
                                                                                           cGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGCATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGc                                                            <  1:3209115/150‑1 (MQ=255)
                                                                                            gatcacATTGAGCACTGCGAGGCCTGCGGCACAGGAAACCGGATTGCCACCGAAGGTGTTGAAATATTCCATGCCATTGTTAAAGCTGTCTGCCACTTCCCGCGTTGTGACAACCGCAGACATTGGATGACCGTTACCAATGGGTTTGcc                                                           <  1:2633913/145‑1 (MQ=255)
                                                                                                gCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGCATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATg                                                       <  2:3159581/150‑1 (MQ=255)
                                                                                                  aTCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGCATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCt                                                     >  2:4414818/1‑150 (MQ=255)
                                                                                                   ttgAGCACTGCGAGGCCTGCGGCACAGGAAACCGGATTGCCACCGAAGGTGTTGAAATATTCCATGCCATTGTTAAAGCTGTCTGCCACTTCCCGCGTTGTGACAACCGCAGACATTGGATGACCGTTACCAATGGGTTTGCCCAtagtg                                                    >  2:1060468/3‑146 (MQ=255)
                                                                                                   tCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGCATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGtt                                                                          <  1:2402746/128‑1 (MQ=255)
                                                                                                   tCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGCATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTg                                                    >  1:1176475/1‑150 (MQ=255)
                                                                                                    cGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGCATc                                                                                                                       <  2:2213320/82‑1 (MQ=255)
                                                                                                        tACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGCATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGc                                               >  2:1103992/1‑150 (MQ=255)
                                                                                                          ccGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGCATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCAt                                             <  1:4405346/150‑1 (MQ=255)
                                                                                                           cGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGCATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATc                                            >  2:4110914/1‑150 (MQ=255)
                                                                                                              cAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGCATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCgg                                         <  1:2480007/150‑1 (MQ=255)
                                                                                                              cAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGCATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCgg                                         <  1:608532/150‑1 (MQ=255)
                                                                                                              cAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGCATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCgg                                         >  1:4580501/1‑150 (MQ=255)
                                                                                                                     aCCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGCATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACg                                  >  1:2518073/1‑150 (MQ=255)
                                                                                                                           cATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCCTGCCGTTGTTGAAAGCATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTgg                            <  1:4088603/150‑1 (MQ=255)
                                                                                                                           cATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGCATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTgg                            <  1:2677806/150‑1 (MQ=255)
                                                                                                                              gACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGCATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGtc                         >  1:3530070/1‑150 (MQ=255)
                                                                                                                                 acCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGCATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCg                      >  1:2212589/1‑150 (MQ=255)
                                                                                                                                 acCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGCATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCg                      >  2:3234878/1‑150 (MQ=255)
                                                                                                                                   cGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGCATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATg                                                                   >  1:3257519/1‑103 (MQ=255)
                                                                                                                                   cGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGCATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGaa                    <  1:1882206/150‑1 (MQ=255)
                                                                                                                                       ttgccaccGAAGGTGTTGAAATATTCCATGCCATTGTTAAAGCTGTCTGCCACTTCCCGCGTTGTGACAACCGCAGACATTGGATGACCGTTACCAATGGGTTTGCCCATAGTGACGATGTCAGGGACAACATCCTGCAACTCGAAAGcc                <  1:2170854/143‑1 (MQ=255)
                                                                                                                                       ttACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGCATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGcc                <  1:4694406/150‑1 (MQ=255)
                                                                                                                                        tACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGCATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGccc               >  1:1457912/1‑150 (MQ=255)
                                                                                                                                          ccGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGCATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCccac             <  2:3898604/150‑1 (MQ=255)
                                                                                                                                          ccGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGCATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCccac             <  1:2823616/150‑1 (MQ=255)
                                                                                                                                            ggCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGCATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCccacca           >  2:2445079/1‑150 (MQ=255)
                                                                                                                                            accGAAGGTGTTGAAATATTCCATGCCATTGTTAAAGCTGTCTGCCACTTCCCGCGTTGTGACAACCGCAGACATTGGATGACCGTTACCAATGGGTTTGCCCATAGTGACGATGTCAGGGACAACATCCTGCAACTCGAAAGCccacca           >  2:796280/3‑150 (MQ=255)
                                                                                                                                               gAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGCATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTg        >  2:4592070/1‑150 (MQ=255)
                                                                                                                                                  ggTGTTGAAATATTCCATGCCATTGTTAAAGCTGTCTGCCACTTCCCGCGTTGTGACAACCGCAGACATTGGATGACCGTTACCAATGGGTTTGCCCATAGTGACGATGTCAGGGACAACATCCTGCAACTCGAAAGCCCACCAGTGGCt     <  2:1732259/150‑1 (MQ=255)
                                                                                                                                                  ggTATTGAAATACTCCATGCCGTTGTTGAAAGCATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCt     <  1:1700108/150‑1 (MQ=255)
                                                                                                                                                     aTTGAAATACTCCATGCCGTTGTTGAAAGCATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGcc  >  1:717244/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
CCAATCACGTCGAATTGCTGCTGAAGCTTGCGCAAACCCTCGAGCAGGTAATGGCCGACACTCAAGGCGTTTTCTTTCAACTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGCATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCC  >  NC_002947/3803078‑3803376

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