Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,813,962 T→A 10.1% Q190L (CAA→CTA)  ligB ← DNA ligase, NAD(+)‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,813,9620TA10.1% 966.1 / 79.0 308Q190L (CAA→CTA) ligBDNA ligase, NAD(+)‑dependent
Reads supporting (aligned to +/- strand):  ref base T (148/129);  new base A (15/16);  total (163/145)
Fisher's exact test for biased strand distribution p-value = 7.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

Only 100 of 330 total aligned reads displayed.

AGCTCTTTTAAACGATCAGACATTAACTGCGGTCCATCCGGCCATGCCCAGACAAAAACGCCCAGAGAATTCAGCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGCGGCACAGCGGAAATTAAGCTAACTTTCTGCGTCCAGTCCTCGCCTTTCAGGCCATTACCGCGACTGATTGCTTT  >  CP009273/3813815‑3814110
                                                                                                                                                   |                                                                                                                                                    
aGCTCTTTTAAACGATCAGACATTAACTGCGGTCCATCCGGCCATGCCCAGACAAAAACGCCCAGAGAATTCAGCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATttgttgttg                                                                                                                                                     <  2:1236013/149‑1 (MQ=255)
  ctctTTTAAACGATCAGACATTAACTGCGGTCCATCCGGCCATGCCCAGACAAAAACGCCCAGAGAATTCAGCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGa                                                                                                                                                   <  1:3915878/149‑1 (MQ=255)
    ctTTTAAACGATCAGACATTAACTGCGGTCCATCCGGCCATGCCCAGACAAAAACGCCCAGAGAATTCAGCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGata                                                                                                                                                 <  1:1122730/149‑1 (MQ=255)
    ctTTTAAACGATCAGACATTAACTGCGGTCCATCCGGCCATGCCCAGACAAAAACGCCCAGAGAATTCAGCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGata                                                                                                                                                 <  2:4232942/149‑1 (MQ=255)
     ttttAAACGATCAGACATTAACTGCGGTCCATCCGGCCATGCCCAGACAAAAACGCCCAGAGAATTCAGCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTAGGatat                                                                                                                                                <  2:133511/149‑1 (MQ=255)
     ttttAAACGAACAGACATTAACTGCGGTCCATCCGGCAATGCCCAGACAAAAACGCCCAGAGAATTCAGCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGatat                                                                                                                                                <  2:2835281/149‑1 (MQ=255)
      tttAAACGATCAGACATTAACTGCGGTCCATCCGGCCATGCCCAGACAAAAACGCCCAGAGAATTCAGCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATg                                                                                                                                               >  2:892558/1‑149 (MQ=255)
       ttAAACGATCAGACATTAACTGCGGTCCATCCGGCCATGCCCAGACAAAAACGCCCAGAGAATTCAGCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGc                                                                                                                                              <  1:675411/149‑1 (MQ=255)
        tAAACGATCAGACATTAACTGCGGTCCATCCGGCCATGCCCAGACAAAAACGCCCAGAGAATTCAGCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGcc                                                                                                                                             >  1:3185682/1‑149 (MQ=255)
        tAAACGATCAGACATTAACTGCGGTCCATCCGGCCATGCCCAGACAAAAACGCCCAGAGAATTCAGCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGcc                                                                                                                                             >  2:1757358/1‑149 (MQ=255)
        tAAACGATCAGACATTAACTGCGGTCCATCCGGCCATGCCCAGACAAAAACGCCCAGAGAATTCAGCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTGTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGcc                                                                                                                                             <  1:4369300/149‑1 (MQ=255)
         aaaCGATCAGACATTAACTGCGGTCCATCCGGCCATGCCCAGACAAAAACGCCCAGAGAATTCAGCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGccc                                                                                                                                            <  1:2764559/149‑1 (MQ=255)
          aaCGATCAGACATTAACTGCGGTCCATCCGGCCAGGCCCAGACAAAAACGCCCAGAGAATTCAGCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGcccc                                                                                                                                           <  1:464252/149‑1 (MQ=255)
            cGATCAGACATTAACTGCGGTCCATCCGGCCATGCCCAGACAAAAACGCCCAGAGAATTCAGCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTc                                                                                                                                         <  2:658799/149‑1 (MQ=255)
                cAGACATTAACTGCGGTCCATCCGGCCATGCCCAGACAAAAACGCCCAGAGAATTCAGCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTAGGATATGCCCCTCgcgc                                                                                                                                     >  1:1803676/1‑149 (MQ=255)
                  gACATTAACTGCGGTCCATCCGGCCATGCCCAGACAAAAACGCCCAGAGAATTCAGCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTg                                                                                                                                   >  1:648993/1‑149 (MQ=255)
                     aTTAACTGCGGTCCATCCGGCCATGCCCAGACAAAAACGCCCAGAGAATTCAGCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGgag                                                                                                                                <  2:3824424/149‑1 (MQ=255)
                              ggTCCATCCGGCCATGCCCAGACAAACACGCCCAGAGAATTCAGCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTc                                                                                                                       <  2:4090201/149‑1 (MQ=255)
                               gTCCATCCGGCCATGCCCAGACAAACACGCCCAGAGAATTCAGCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTcc                                                                                                                      <  1:3695527/149‑1 (MQ=255)
                               gTCCATCCGGCCATGCCCAGACAAAAACGCCCAGAGAATTCAGCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCCTTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTcc                                                                                                                      >  2:4679618/1‑149 (MQ=255)
                                 ccATCCGGCCATGCCCAGACAAAAACGCCCAGAGAATTCAGCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTcccc                                                                                                                    >  2:3902872/1‑149 (MQ=255)
                                  cATCCGGCCATGCCCAGACAAAAACGCCCAGAGAATTCAGCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTccccc                                                                                                                   <  2:565513/149‑1 (MQ=255)
                                   aTCCGGCCATGCCCAGACAAAAACGCCCAGAGAATTCAGCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCt                                                                                                                  >  1:1545610/1‑149 (MQ=255)
                                   aTCCGGCCATGCCCAGACAAAAACGCCCAGAGAATTCAGCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCt                                                                                                                  >  2:4780577/1‑149 (MQ=255)
                                    tCCGGCCATGCCCAGACAAAAACGCCCAGAGAATTCAGCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTg                                                                                                                 >  1:2507242/1‑149 (MQ=255)
                                      cGGCCATGCCCAGACAAAAACGCCCAGAGAATTCAGCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGaa                                                                                                               >  2:3837291/1‑149 (MQ=255)
                                      cGGCCATGCCCAGACAAAAACGCCCAGAGAATTCAGCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATGTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGaa                                                                                                               >  2:870400/1‑149 (MQ=255)
                                             gCCCAGACAAAAACGCCCAGAGAATTCAGCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACt                                                                                                        <  2:1413152/149‑1 (MQ=255)
                                                 aGACAAAAACGCCCAGAGAATTCAGCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTAttgttg                                                                                                    <  1:1947542/149‑1 (MQ=255)
                                                  gACAAAAACGCCCAGAGAATTCAGCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTgg                                                                                                   <  2:814349/149‑1 (MQ=255)
                                                   aCAAAAACGCCCAGAGAATTCAGCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTACCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGc                                                                                                  >  1:2975349/1‑149 (MQ=255)
                                                    cAAAAACGCCCAGAGAATTCAGCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCt                                                                                                 >  1:1535865/1‑149 (MQ=255)
                                                      aaaaCGCCCAGAGAATTCAGCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTaa                                                                                               >  1:338553/1‑149 (MQ=255)
                                                      aaaaCGCCCAGAGAATTCAGCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAACAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTaa                                                                                               >  1:2622540/1‑149 (MQ=255)
                                                       aaaCGCCCAGAGAATTCAGCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTaaa                                                                                              >  1:663016/1‑149 (MQ=255)
                                                        aaCGCCCAGAGAATTCAGCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAg                                                                                             >  2:4500513/1‑149 (MQ=255)
                                                        aaCGCCCAGAGAATTCAGCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCACCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAg                                                                                             <  2:1097254/149‑1 (MQ=255)
                                                         aCGCCCAGAGAATTCAGCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAgg                                                                                            <  1:3833707/149‑1 (MQ=255)
                                                             ccAGAGAATTCAGCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCg                                                                                        >  1:547146/1‑149 (MQ=255)
                                                              cAGAGAATTCAGCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTAGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGc                                                                                       >  2:723081/1‑149 (MQ=255)
                                                               agagAATTCAGCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCt                                                                                      >  1:193936/1‑149 (MQ=255)
                                                                    aTTCAGCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACgg                                                                                 >  1:3980584/1‑149 (MQ=255)
                                                                    aTTCAGCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACgg                                                                                 >  1:2449021/1‑149 (MQ=255)
                                                                     ttCAGCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTAGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGt                                                                                >  1:698355/1‑149 (MQ=255)
                                                                      tCAGCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGtt                                                                               <  2:4588591/149‑1 (MQ=255)
                                                                       cAGCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGttt                                                                              <  2:377796/149‑1 (MQ=255)
                                                                        aGCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTg                                                                             >  2:1491791/1‑149 (MQ=255)
                                                                        aGCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTg                                                                             <  2:1664831/149‑1 (MQ=255)
                                                                         gCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGc                                                                            >  1:3734357/1‑149 (MQ=255)
                                                                         gCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGc                                                                            >  1:1130867/1‑149 (MQ=255)
                                                                             gtCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGCGGca                                                                        >  1:4321490/1‑149 (MQ=255)
                                                                              tCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGCGGcac                                                                       <  1:1275725/149‑1 (MQ=255)
                                                                               cGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTAGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGCGGcaca                                                                      <  1:583229/149‑1 (MQ=255)
                                                                                 cTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGCGGCACAGc                                                                    <  1:3602492/149‑1 (MQ=255)
                                                                                  tATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGCGGCACAGCg                                                                   >  2:3057029/1‑149 (MQ=255)
                                                                                  tATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGCGGCACAGCg                                                                   <  2:3723428/149‑1 (MQ=255)
                                                                                    tCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGCGGCACAGCGGa                                                                 <  1:837495/149‑1 (MQ=255)
                                                                                          tGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGCGGCACAGCGGAAATTaa                                                           <  1:2191738/149‑1 (MQ=255)
                                                                                           ggCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGCGGCACAGCGGAAATTAAg                                                          >  2:4081166/1‑149 (MQ=255)
                                                                                            gcgcATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGCGGCACAGCGGAAATTAAGc                                                         >  2:1393400/1‑149 (MQ=255)
                                                                                            gcgcATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTAGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGCGGCACAGCGGAAATTAAGc                                                         <  1:2521447/149‑1 (MQ=255)
                                                                                              gcATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGCGGCACAGCGGAAATTAAGCTa                                                       >  1:4031934/1‑149 (MQ=255)
                                                                                               catcatCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGCGGCACAGCGGAAATTAAGCTaa                                                      >  2:3482823/1‑149 (MQ=255)
                                                                                                atcatcAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGCGGCACAGCGGAAATTAAGCTAAc                                                     <  1:873694/149‑1 (MQ=255)
                                                                                                  catcaAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGCGGCACAGCGGAAATTAAGCTAACtt                                                   <  2:3615088/149‑1 (MQ=255)
                                                                                                  catcaAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGCGGCACAGCGGAAATTAAGCTAACtt                                                   <  2:4392919/149‑1 (MQ=255)
                                                                                                   atcaAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGCGGCACAGCGGAAATTAAGCTAACttt                                                  >  1:3479634/1‑149 (MQ=255)
                                                                                                   atcaAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGCGGCACAGCGGAAATTAAGCTAACttt                                                  <  1:2065887/149‑1 (MQ=255)
                                                                                                    tcaAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGCGGCACAGCGGAAATTAAGCTAACTTTc                                                 <  2:2465015/149‑1 (MQ=255)
                                                                                                     caAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGCGGCACAGCGGAAATTAAGCTAACTTTCt                                                <  1:3093334/149‑1 (MQ=255)
                                                                                                      aaGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGCGGCACAGCGGAAATTAAGCTAACTTTCTg                                               >  1:474075/1‑149 (MQ=255)
                                                                                                          cAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGCGGCACAGCGGAAATTAAGCTAACTTTCTGCGTc                                           <  1:2340182/149‑1 (MQ=255)
                                                                                                          cAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTAGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGCGGCACAGCGGAAATTAAGCTAACTTTCTGCGTc                                           >  2:2373233/1‑149 (MQ=255)
                                                                                                          cAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTAGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGCGGCACAGCGGAAATTAAGCTAACTTTCTGCGTc                                           <  1:2030576/149‑1 (MQ=255)
                                                                                                          cAGCAACTTTTGCGCGGGCATTTAGTCCCCCCATTTGTTGTAGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGCGGCACAGCGGAAATTAAGCTAACTTTCTGCGTc                                           <  2:4124112/149‑1 (MQ=255)
                                                                                                           aGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGCGGCACAGCGGAAATTAAGCTAACTTTCTGCGTcc                                          <  1:2427507/149‑1 (MQ=255)
                                                                                                            gCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGCGGCACAGCGGAAATTAAGCTAACTTTCTGCGTCCa                                         >  1:3746837/1‑149 (MQ=255)
                                                                                                            gCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGCGGCACAGCGGAAATTAAGCTAACTTTCTGCGTCCa                                         >  1:4209117/1‑149 (MQ=255)
                                                                                                               aCTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGCGGCACAGCGGAAATTAAGCTAACTTTCTGCGTCCAGTc                                      <  1:3871238/149‑1 (MQ=255)
                                                                                                               aCTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGCGGCACAGCGGAAATTAAGCTAACTTTCTGCGTCCAGTc                                      >  2:1627710/1‑149 (MQ=255)
                                                                                                               aCTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGCGGCACAGCGGAAATTAAGCTAACTTTCTGCGTCCAGTc                                      >  1:4352716/1‑149 (MQ=255)
                                                                                                                      cgcgGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGCGGCACAGCGGAAATTAAGCTAACTTTCTGCGTCCAGTCCTCGCCt                               <  1:1836728/149‑1 (MQ=255)
                                                                                                                       gcgGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGCGGCACAGCGGAAATTAAGCTAACTTTCTGCGTCCAGTCCTCGCCtt                              <  2:2507242/149‑1 (MQ=255)
                                                                                                                        cgGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGCGGCACAGCGGAAATTAAGCTAACTTTCTGCGTCCAGTCCTCGCCttt                             <  2:2083108/149‑1 (MQ=255)
                                                                                                                            cATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGCGGCACAGCGGAAATTAAGCTAACTTTCTGCGTCCAGTCCTCGCCTTTCAgg                         >  2:3129966/1‑149 (MQ=255)
                                                                                                                             aTTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGCGGCACAGCGGAAATTAAGCTAACTTTCTGCGTCCAGTCCTCGCCTTTCAGGc                        <  1:1606946/149‑1 (MQ=255)
                                                                                                                             aTTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGCGGCACAGCGGAAATTAAGCTAACTTTCTGCGTCCAGTCCTCGCCTTTCAGGc                        >  1:945238/1‑149 (MQ=255)
                                                                                                                                tATTCCCCCCATTTGTTGTAGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGCGGCACAGCGGAAATTAAGCTAACTTTCTGCGTCCAGTCCTCGCCTTTCAGGCCAt                     <  1:3334908/149‑1 (MQ=255)
                                                                                                                                 agtCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGCGGCACAGCGGAAATTAAGCTAACTTTCTGCGTCCAGTCCTCGCCTTTCAGGCCAtt                    <  2:2172435/147‑1 (MQ=255)
                                                                                                                                 aTTCCCCCCATTTGGTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGCGGCACAGCGGAAATTAAGCTAACTTTCTGCGTCCAGTCCTCGCCTTTCAGGCCAtt                    <  2:3354784/149‑1 (MQ=255)
                                                                                                                                   tCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGCGGCACAGCGGAAATTAAGCTAACTTTCTGCGTCCAGTCCTCGCCTTTCAGGCCATTAc                  <  1:4284333/149‑1 (MQ=255)
                                                                                                                                   tCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGCGGCACAGCGGAAATTAAGCTAACTTTCTGCGTCCAGTCCTCGCCTTTCAGGCCATTAc                  <  1:2274983/149‑1 (MQ=255)
                                                                                                                                    ccccccATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGCGGCACAGCGGAAATTAAGCTAACTTTCTGCGTCCAGTCCTCGCCTTTCAGGCCATTAcc                 <  1:3320687/149‑1 (MQ=255)
                                                                                                                                     cccccATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGCGGCACAGCGGAAATTAAGCTAACTTTCTGCGTCCAGTCCTCGCCTTTCAGGCCATTACcg                <  1:797858/149‑1 (MQ=255)
                                                                                                                                     cccccATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGCGGCACAGCGGAAATTAAGCTAACTTTCTGCGTCCAGTCCTCGCCTTTCAGGCCATTACcg                >  2:3276942/1‑149 (MQ=255)
                                                                                                                                     cccccATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCGAACGGTTTGCGGCACAGCGGAAATTAAGCTAACTTTCTGCGTCCAGTCCTCGCCTTTCAGGCCATTACcg                <  2:1812688/149‑1 (MQ=255)
                                                                                                                                           tttGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGCGGCACAGCGGAAATTAAGCTAACTTTCTGCGTCCAGTCCTCGCCTTTCAGGCCATTACCGCGACTg          >  1:2592295/1‑149 (MQ=255)
                                                                                                                                            ttgttgttgGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGCGGCACAGCGGAAATTAAGCTAACTTTCTGCGTCCAGTCCTCGCCTTTCAGGCCATTACCGCGACTGa         <  1:1014531/149‑1 (MQ=255)
                                                                                                                                              gttgttgGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGCGGCACAGCGGAAATTAAGCTAACTTTCTGCGTCCAGTCCTCGCCTTTCAGGCCATTACCGCGACTGAtt       <  2:2622540/149‑1 (MQ=255)
                                                                                                                                                 gttgGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGCGGCACAGCGGAAATTAAGCTAACTTTCTGCGTCCAGTCCTCGCCTTTCAGGCCATTACCGCGACTGATTGCt    >  1:1448199/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
AGCTCTTTTAAACGATCAGACATTAACTGCGGTCCATCCGGCCATGCCCAGACAAAAACGCCCAGAGAATTCAGCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGCGGCACAGCGGAAATTAAGCTAACTTTCTGCGTCCAGTCCTCGCCTTTCAGGCCATTACCGCGACTGATTGCTTT  >  CP009273/3813815‑3814110

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: