Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,000,498 T→A 6.9% R425R (CGT→CGA)  recQ → ATP‑dependent DNA helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,000,4980TA 6.9% 1076.5 / 52.0 332R425R (CGT→CGA) recQATP‑dependent DNA helicase
Reads supporting (aligned to +/- strand):  ref base T (143/166);  new base A (13/10);  total (156/176)
Fisher's exact test for biased strand distribution p-value = 3.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 342 total aligned reads displayed.

AAGCGCAAACTTGCCGTCGTCTGGTATTGCTGAACTATTTTGGCGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGTATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGTGATCCGCC  >  CP009273/4000350‑4000646
                                                                                                                                                    |                                                                                                                                                    
aaGCGCAAACTTGCCGTCGTCTGGTATTGCTGAACTATTTTGGCGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCgt                                                                                                                                                      <  1:2216159/149‑1 (MQ=255)
      aaaCTTGCCGTCGTCTGGTATTGCTGAACTATTTTGGCGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAAt                                                                                                                                                >  2:1161409/1‑149 (MQ=255)
       aaCTTGCCGTCGTCTGGTATTGCTGAACTATTTTGGCGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGACTc                                                                                                                                               >  2:2813231/1‑149 (MQ=255)
       aaCTTGCCGTCGTCTGGTATTGCTGAACTATTTTGGCGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATc                                                                                                                                               >  2:1785225/1‑149 (MQ=255)
       aaCTTGCCGTCGTCTGGTATTGCTGAACTATTTTGGCGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATc                                                                                                                                               >  2:420026/1‑149 (MQ=255)
        aCTTGCCGTCGTCTGGTATTGCTGAACTATTTTGGCGAAGGGCGTCAGGAGCCGTGCGGGACCTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCa                                                                                                                                              >  2:1502364/1‑149 (MQ=255)
            gCCGTCGTCTGGTATTGCTGAACTATTTTGGCGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCgg                                                                                                                                          <  1:1873496/149‑1 (MQ=255)
             ccGTCGTCTGGTATTGCTGAACTATTTTGGCGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGt                                                                                                                                         <  1:928498/149‑1 (MQ=255)
              cgtcgtCTGGTATTGCTGAACTATTTTGGCGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGtt                                                                                                                                        <  2:2422676/149‑1 (MQ=255)
              cgtcgtCTGGTATTGCTGAACTATTTTGGCGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGtt                                                                                                                                        <  1:3830848/149‑1 (MQ=255)
                tcgtcTGGTATTGCTGAACTATTTTGGCGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTg                                                                                                                                      >  1:4348262/1‑149 (MQ=255)
                   tcTGGTATTGCTGAACTATTTTGGCGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGa                                                                                                                                   <  2:3580320/149‑1 (MQ=255)
                    cTGGTATTGCTGAACTATTTTGGCGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGAt                                                                                                                                  <  2:1721184/149‑1 (MQ=255)
                     tGGTATTGCTGAACTATTTTGGCGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATg                                                                                                                                 <  2:82807/149‑1 (MQ=255)
                      ggTATTGCTGAACTATTTTGGCGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATgg                                                                                                                                >  2:4433676/1‑149 (MQ=255)
                         aTTGCTGAACTATTTTGGCGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGtt                                                                                                                             >  2:4611402/1‑149 (MQ=255)
                         aTTGCTGAACTATTTTGGCGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGtt                                                                                                                             <  1:3068931/149‑1 (MQ=255)
                           tGCTGAACTATTTTGGCGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGAGTGAATCAGCGGTTTGGGATGGGTTAt                                                                                                                           <  1:1405480/149‑1 (MQ=255)
                             cTGAACTATTTTGGCGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTAtgt                                                                                                                         <  1:867831/149‑1 (MQ=255)
                                 aCTATTTTGGCGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATgtggtg                                                                                                                     >  2:4328190/1‑149 (MQ=255)
                                    aTTTTGGCGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGaa                                                                                                                  >  2:1306465/1‑149 (MQ=255)
                                    aTTTTGGCGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGAGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGaa                                                                                                                  >  2:3645293/1‑149 (MQ=255)
                                     ttttGGCGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCCTTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAg                                                                                                                 <  2:3140766/149‑1 (MQ=255)
                                       ttGGCGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTg                                                                                                               <  1:2458358/149‑1 (MQ=255)
                                        tGGCGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGa                                                                                                              <  2:3398178/149‑1 (MQ=255)
                                        tGGCGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGAGATCTGCCTCGATCCGCCGAAACAGGACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGa                                                                                                              <  2:3253683/149‑1 (MQ=255)
                                         ggCGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGAt                                                                                                             <  1:3453144/149‑1 (MQ=255)
                                           cGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTc                                                                                                           <  2:3317024/149‑1 (MQ=255)
                                              aGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCgtg                                                                                                        >  1:342537/1‑149 (MQ=255)
                                               gggCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGCTGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCgtgg                                                                                                       <  2:1434032/149‑1 (MQ=255)
                                               gggCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCgtgg                                                                                                       >  1:3141655/1‑149 (MQ=255)
                                               gggCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCgtgg                                                                                                       <  2:2642883/149‑1 (MQ=255)
                                                      aGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCtaat                                                                                                >  2:260509/1‑149 (MQ=255)
                                                         aGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAAc                                                                                             >  2:3097312/1‑149 (MQ=255)
                                                          gCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAAcc                                                                                            >  1:565421/1‑149 (MQ=255)
                                                          gCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAAcc                                                                                            >  1:1706805/1‑149 (MQ=255)
                                                          gCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAAcc                                                                                            <  1:1672299/149‑1 (MQ=255)
                                                          gCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAAcc                                                                                            <  2:140095/149‑1 (MQ=255)
                                                          gCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGAGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAAcc                                                                                            >  1:4107089/1‑149 (MQ=255)
                                                             gTGCGGGAACTGCGATATCTGCCTCGATCCGCTGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGc                                                                                         <  1:1768216/149‑1 (MQ=255)
                                                             gTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGc                                                                                         <  2:730864/149‑1 (MQ=255)
                                                             gTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGc                                                                                         <  1:2571644/149‑1 (MQ=255)
                                                                 gggAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGTAt                                                                                     <  2:3858659/149‑1 (MQ=255)
                                                                  ggAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGTATc                                                                                    <  1:4743127/149‑1 (MQ=255)
                                                                    aaCTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGTATCcg                                                                                  <  2:3805778/149‑1 (MQ=255)
                                                                       tGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGTATCCGCGa                                                                               >  1:2092572/1‑149 (MQ=255)
                                                                       tGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGTATCCGCGa                                                                               <  1:3210169/149‑1 (MQ=255)
                                                                        gCGATATCTGCCTCGATCCGCCGAAAAAGTACTACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGTATCCGCGAc                                                                              <  1:1295568/149‑1 (MQ=255)
                                                                          gATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGTATCCGCGACTa                                                                            >  2:4524229/1‑149 (MQ=255)
                                                                          gATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGTATCCGCGACTa                                                                            >  2:549798/1‑149 (MQ=255)
                                                                             atCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGTATCCGCGACTATgg                                                                         <  1:4187738/149‑1 (MQ=255)
                                                                             atCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGAGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGTATCCGCGACTATgg                                                                         <  1:4617286/149‑1 (MQ=255)
                                                                              tCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGTATCCGCGACTATGGt                                                                        <  2:3853480/149‑1 (MQ=255)
                                                                              tCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGTATCCGCGACTATGGt                                                                        <  2:4587725/149‑1 (MQ=255)
                                                                                tGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGTATCCGCGACTATGGTCa                                                                      >  1:4554023/1‑149 (MQ=255)
                                                                                 gCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGTATCCGCGACTATGGTCAt                                                                     <  2:1748050/149‑1 (MQ=255)
                                                                                        tCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGTATCCGCGACTATGGTCATGACAAAc                                                              <  1:3477674/149‑1 (MQ=255)
                                                                                         ccgccgAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGTATCCGCGACTATGGTCATGACAAACt                                                             <  2:4759633/149‑1 (MQ=255)
                                                                                         ccgcagAAACAGTACGACGGTTCAAACGATGCTCAGATGGCCCTTTGCACCTTTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGTATCCGCGACTATGGTCATGACAAACt                                                             <  2:2013069/149‑1 (MQ=255)
                                                                                            ccgAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGTATCCGCGACTATGGTCATGACAAACTGaa                                                          <  1:2312257/149‑1 (MQ=255)
                                                                                            ccgAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGTATCCGCGACTATGGTCATGACAAACTGaa                                                          <  1:2588622/149‑1 (MQ=255)
                                                                                              gAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGTATCCGCGACTATGGTCATGACAAACTGAAAg                                                        <  1:362908/149‑1 (MQ=255)
                                                                                              gAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGTATCCGCGACTATGGTCATGACAAACTGAAAg                                                        >  2:209794/1‑149 (MQ=255)
                                                                                               aaaCAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGTATCCGCGACTATGGTCATGACAAACTGAAAGt                                                       >  2:1563013/1‑149 (MQ=255)
                                                                                                aaCAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGTATCCGCGACTATGGTCATGACAAACTGAAAGTc                                                      >  2:1489985/1‑149 (MQ=255)
                                                                                                 aCAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGTATCCGCGACTATGGTCATGACAAACTGAAAGTCt                                                     >  1:2969119/1‑149 (MQ=255)
                                                                                                   aGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGTATCCGCGACTATGGTCATCACAAACTGAAAGTCTAt                                                   <  2:4371428/149‑1 (MQ=255)
                                                                                                      acgacgGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGTATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGc                                                >  1:518458/1‑149 (MQ=255)
                                                                                                        gacgGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGTATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCAt                                              <  1:3901518/149‑1 (MQ=255)
                                                                                                         acgGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGTATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATg                                             <  1:4014238/149‑1 (MQ=255)
                                                                                                          cgGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGTATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATgg                                            <  1:3171933/149‑1 (MQ=255)
                                                                                                           ggTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGTATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATggg                                           <  1:237098/149‑1 (MQ=255)
                                                                                                            gTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGAGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGTATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGc                                          >  1:4121600/1‑149 (MQ=255)
                                                                                                               cAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGTATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGt                                       >  1:2204031/1‑149 (MQ=255)
                                                                                                                 aCCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGTATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGa                                     >  1:1080018/1‑149 (MQ=255)
                                                                                                                 aCCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGTATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGa                                     >  2:2911414/1‑149 (MQ=255)
                                                                                                                 aCCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGTATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGa                                     <  2:3983482/149‑1 (MQ=255)
                                                                                                                     aTGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTGATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGTATCCGCGACTCTGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATaaa                                 >  2:1605260/1‑149 (MQ=255)
                                                                                                                      tGCTCAGATTGCCCTTTCCACCATTGGTCGAGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGTATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATaaaa                                <  2:776894/149‑1 (MQ=255)
                                                                                                                       gCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGTATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAg                               <  2:2314008/149‑1 (MQ=255)
                                                                                                                         tCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGTATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGcc                             >  2:1975395/1‑149 (MQ=255)
                                                                                                                           aGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGTATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCAt                           >  1:1010440/1‑149 (MQ=255)
                                                                                                                            gATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGTATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATg                          >  2:3955562/1‑149 (MQ=255)
                                                                                                                                gCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATTGGTTATGTTGTGGAAGTGATTCTTGGTGCTAATAACCAGCGTATCCGCGACTATGGTCATGACAAACTGAAAGTATATGGCATGGGCCGTGATAAAAGCCATGAACa                      <  2:915035/149‑1 (MQ=255)
                                                                                                                                gCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGTATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACa                      >  2:240338/1‑149 (MQ=255)
                                                                                                                                 cccTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGTATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACAt                     <  1:4537125/149‑1 (MQ=255)
                                                                                                                                  ccTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGTATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACAtt                    <  1:3321524/149‑1 (MQ=255)
                                                                                                                                    tttCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGTATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTgg                  >  1:4404586/1‑149 (MQ=255)
                                                                                                                                    tttCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGTATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTgg                  >  2:493566/1‑149 (MQ=255)
                                                                                                                                     ttCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGTATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTggg                 <  1:2474866/149‑1 (MQ=255)
                                                                                                                                      tCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGTATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGt                <  1:1502364/149‑1 (MQ=255)
                                                                                                                                       ccaccaTTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGTATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTg               <  1:2146418/149‑1 (MQ=255)
                                                                                                                                        caccaTTTGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGTATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGa              <  1:4484239/149‑1 (MQ=255)
                                                                                                                                        caccaTTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGTATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGa              >  1:3640135/1‑149 (MQ=255)
                                                                                                                                           caTTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGTATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCg           >  1:3247997/1‑149 (MQ=255)
                                                                                                                                              tGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGTATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGTGa        >  1:4047840/1‑149 (MQ=255)
                                                                                                                                                 tCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGTATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGTGATcc     >  1:650482/1‑149 (MQ=255)
                                                                                                                                                   gtgtGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGTATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGTGATCCGc   >  1:3991112/1‑149 (MQ=255)
                                                                                                                                                    tgtgAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGTATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGTGATCCGcc  <  1:469734/149‑1 (MQ=255)
                                                                                                                                                    tgtgAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGTATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGTGATCCGcc  <  1:1048474/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AAGCGCAAACTTGCCGTCGTCTGGTATTGCTGAACTATTTTGGCGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGTATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGTGATCCGCC  >  CP009273/4000350‑4000646

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: ‑