Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,999,888 A→G 6.2% K222R (AAA→AGA)  recQ → ATP‑dependent DNA helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,999,8880AG 6.2% 1051.5 / 43.1 305K222R (AAA→AGA) recQATP‑dependent DNA helicase
Reads supporting (aligned to +/- strand):  ref base A (152/134);  new base G (9/10);  total (161/144)
Fisher's exact test for biased strand distribution p-value = 6.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.18e-01

Only 100 of 311 total aligned reads displayed.

CCGACGCTGCCGTTTATGGCGCTGACCGCCACAGCCGACGACACCACGCGCCAGGATATCGTGCGCCTGCTGGGGCTGAACGATCCGCTGATTCAAATCAGCAGTTTTGACCGTCCGAATATTCGCTACATGCTGATGGAGAAGTTCAAACCGCTCGATCAGTTGATGCGCTACGTGCAGGAACAGCGCGGTAAGTCAGGCATTATCTACTGCAACAGCCGCGCGAAAGTAGAAGACACCGCTGCGCGCCTGCAAAGCAAGGGAATTAGCGCGGCGGCCTATCATGCCGGGCTGGAA  >  CP009273/3999740‑4000036
                                                                                                                                                    |                                                                                                                                                    
 cGACGCTGCCGTTTATGGCGCTGACCGCCACAGCCGACGACACCACGCGCCAGGATATCGTGCGCCTGCTGGGGCTGAACGATCCGCTGATTCAAATCAGCAGTTTTGACCGTCCGAATATTCGCTACATGCTGATGGAGAAGTTCaaa                                                                                                                                                     <  2:3737895/149‑1 (MQ=255)
  gACGCTGCCGTTTATGGCGCTGACCGCCACAGCCGACGACACCACGCGCCAGGATATCGTGCGCCTGCTGGGGCTGAACGATCCGCTGATTCAAATCAGCAGTTTTGACCGTCCGAATATTCGCTACATGCTGATGGAGAAGTTCAAAc                                                                                                                                                    >  1:5171948/1‑149 (MQ=255)
   aCGCTGCCGTTTATGGCGCTGACCGCCACAGCCGACGACACCACGCGCCAGGATATCGTGCGCCTGCTGGGGCTGAACGATCCGCTGATTCAAATCAGCAGTTTTGACCGTCCGAATATTCGCTACATGCTGATGGAGAAGTTCAGAcc                                                                                                                                                   <  2:4193094/149‑1 (MQ=255)
       tGCCGTTTATGGCGCTGACCGCCACAGCCGACGACACCACGCGCCAGGATATCGTGCGCCTGCTGGGGCTGAACGATCCGCTGATTCAAATCAGCAGTTTTGACCGTCCGAATATTCGCTACATGCTGATGGAGAAGTTCAAACCGCTc                                                                                                                                               <  1:1362485/149‑1 (MQ=255)
       tGCCGTTTATGGCGCTGACCGCCACAGCCGACGACACCACGCGCCAGGATATCGTGCGCCTGCTGGGGCTGAACGATCCGCTGATTCAAATCAGCAGTTTTGACCGTCCGAATATTCGCTACATGCTGATGGAGAAGTTCAAACCGCTc                                                                                                                                               >  2:67180/1‑149 (MQ=255)
         ccGTTTATGGCGCTGACCGCCACAGCCGACGACACCACGCGCCAGGATATCGTGCGCCTGCTGGGGCTGAACGATCCGCTGATTCAAATCAGCAGTTTTGACCGTCCGAATATTCGCTACATGCTGATGGAGAAGTTCAAACCGCTCGa                                                                                                                                             <  2:2354394/149‑1 (MQ=255)
          cGTTTATGGCGCTGACCGCCACAGCCGACGACACCACGCGCCAGGATATCGTGCGCCTGCTGGGGCTGAACGATCCGCTGATTCAAATCAGCAGTTTTGACCGTCCGAATATTCGCTACATGCTGATGGAGAAGTTCAAACCGCTCGAt                                                                                                                                            >  1:1777152/1‑149 (MQ=255)
              tATGGCGCTGACCGCCACAGCCGACGACACCACGCGCCAGGATATCGTGCGCCTGCTGGGGCTGAACGATCCGCTGATTCAAATCAGCAGTTTTGACCGTCCGAATATTCGCTACATGCTGATGGAGAAGTTCAGACCGCTCGATCAGt                                                                                                                                        >  1:1337280/1‑149 (MQ=255)
                tGGCGCTGACCGCCACAGCCGACGACACCACGCGCCAGGATATCGTGCGCCTGCTGGGGCTGAACGATCCGCTGATTCAAATCAGCAGTTTTGACCGTCCGAATATTCGCTACATGCTGATGGAGAAGTTCAAACCGCTCGATCAGTTg                                                                                                                                      >  2:4554874/1‑149 (MQ=255)
                   cgcTGACCGCCACAGCCGACGACACCACGCGCCAGGATATCGTGCGCCTGCTGGGGCTGAACGATCCGCTGATTCAAATCAGCAGTTTTGACCGTCCGAATATTCGCTACATGCTGATGGAGAAGTTCAAACCGCTCGATCAGTTGATg                                                                                                                                   <  2:4695905/149‑1 (MQ=255)
                      tGACCGCCACAGCCGACGACACCACGCGCCAGGATATCGTGCGCCTGCTGGGGCTGAACGATCCGCTGATTCAAATCAGCAGTTTTGACCGTCCGAATATTCGCTACATGCTGATGGAGAAGTTCAAACCGCTCGATCAGTTGATgcgc                                                                                                                                >  1:1458804/1‑149 (MQ=255)
                      tGACCGCCACAGCCGACGACACCACGCGCCAGGATATCGTGCGCCTGCTGGGGCTGAACGATCCGCTGATTCAAATCAGCAGTTTTGACCGTCCGAATATTCGCTACATGCTGATGGAGAAGTTCAAACCGCTCGATCAGTTGATgcgc                                                                                                                                >  2:454826/1‑149 (MQ=255)
                         ccGCCACAGCCGACGACACCACGCGCCAGGATATCGTGCGCCTGCTGGGGCTGAACGATCCGCTGATTCAAATCAGCAGTTTTGACCGTCCGAATATTCGCTACATGCTGATGGAGAAGTTCAAACCGCTCGATCAGTTGATGCGCTAc                                                                                                                             >  2:1557508/1‑149 (MQ=255)
                         ccGCCACAGCCGACGACACCACGCGCCAGGATATCGTGCGCCTGCTGGGGCTGAACGATCCGCTGATTCAAATCAGCAGTTTTGACCGTCCGAATATTCGCTACATGCTGATGGAGAAGTTCAAACCGCTCGATCAGTTGATGCGCTAc                                                                                                                             >  2:814699/1‑149 (MQ=255)
                         ccGCCACAGCCGACGACACCACGCGCCAGGATATCGTGCGCCTGCTGGGGCTGAACGATCCGCTGATTCAAATCAGCAGTTTTGACCGTCCGAATATTCGCTACATGCTGATGGAGAAGTTCAAACCGCTCGATCAGTTGATGCGCTAc                                                                                                                             >  1:788062/1‑149 (MQ=255)
                           gCCACAGCCGACGACACCACGCGCCAGGATATCGTGCGCCTGCTGGGGCTGAACGATCCGCTGATTCAAATCAGCAGTTTTGACCGTCCGAATATTCGCTACATGCTGATGGAGAAGTTCAAACCGCTCGATCAGTTGATGCGCTACGt                                                                                                                           <  2:1668625/149‑1 (MQ=255)
                             cacaGCCGACGACACCACGCGCCAGGATATCGTGCGCCTGCTGGGGCTGAACGATCCGCTGATTCAAATCAGCAGTTTTGACCGTCCGAATATTCGCTACATGCTGATGGAGAAGTTCAAACCGCTCGATCAGTTGATGCGCTACGTGc                                                                                                                         <  2:4476241/149‑1 (MQ=255)
                               caGCCGACGACACCACGCGCCAGGATATCGTGCGCCTGCTGGGGCTGAACGATCCGCTGATTCAAATCAGCAGTTTTGACCGTCCGAATATTCGCTACATGCTGATGGAGAAGTTCAAACCGCTCGATCAGTTGATGCGCTACGTGCAg                                                                                                                       <  2:1455793/149‑1 (MQ=255)
                                 gCCGACGACACCACGCGCCAGGATATCGTGCGCCTGCTGGGGCTGAACGATCCGCTGATTCAAATCAGCAGTTTTGACCGTCCGAATATTCGCTACATGCTGATGGAGAAGTTCAAACCGCTCGATCAGTTGATGCGCTACGTGCAGGa                                                                                                                     >  2:3629392/1‑149 (MQ=255)
                                   cgacgaCACCACGCGCCAGGATATCGTGCGCCTGCTGGGGCTGAACGATCCGCTGATTCAAATCAGCAGTTTTGACCGTCCGAATATTCGCTACATGCTGATGGAGAAGTTCAGACCGCTCGATCAGTTGAt                                                                                                                                    >  2:1817196/1‑132 (MQ=255)
                                   cgacgaCACCACGCGCCAGGATATCGTGCGCCTGCTGGGGCTGAACGATCCGCTGATTCAAATCAGCAGTTTTGACCGTCCGAATATTCGCTACATGCTGATGGAGAAGTTCAGACCGCTCGATCAGTTGAt                                                                                                                                    <  1:1817196/132‑1 (MQ=255)
                                   cgacgaCACCACGCGCCAGGATATCGTGCGCCTGCGGGGGCTGAACGATCCGCTGATTCAAATCAGCAGTTTTGACCGTCCGAATATTCGCTACAGGCTGATGGAGAAGTTCAAACCGCTCGATCAGTTGATGCGCTACGTGCAGGAAc                                                                                                                   <  1:3470855/149‑1 (MQ=255)
                                    gacgacACCACGCGCCAGGATATCGTGCGCCTGCTGGGGCTGAACGATCCGCTGATTCAAATCAGCAGTTTTGACCGTCCGAATATTCGCTACATGCTGATGGAGAAGTTCAAACCGCTCGATCAGTTGATGCGCTACGTGCAGGAACa                                                                                                                  <  2:4321390/149‑1 (MQ=255)
                                     acgacACCACGCGCCAGGATATCGTGCGCCTGCTGGGGCTGAACGATCCGCTGATTCAAATCAGCAGTTTTGACCGTCCGAATATTCGCTACATGCTGATGGAGAAGTTCAAACCGCTCGATCAGTTGATGCGCTACGTGCAGGAACAg                                                                                                                 <  2:3996924/149‑1 (MQ=255)
                                        acacCACGCGCCAGGAGATCGTGCGCCTGCTGGGGCTGAACGATCCGCTGATTCAAATCAGCAGTTTTGACCGTCCGAATATTCGCTACATGCTGATGGAGAAGTTCAAACCGCTCGATCAGTTGATGCGCTACGTGCAGGAACAgcgc                                                                                                              <  2:3540788/149‑1 (MQ=255)
                                         caccacGCGCCAGGATATCGTGCGCCTGCTGGGGCTGAACGATCCGCTGATTCAAATCAGCAGTTTTGACCGTCCGAATATTCGCTACATGCTGATGGAGAAGTTCAAACCGCTCGATCAGTTGATGCGCTACGTGCAGGAACcgcgcg                                                                                                             >  1:3544880/1‑149 (MQ=255)
                                         caccacGCGCCAGGATATCGTGCGCCTGCTGGGGCTGAACGATCCGCTGATTCAAATCAGCAGTTTTGACCGTCCGAATATTCGCTACATGCTGATGGAGAAGTTCAAACCGCTCGATCAGTTGATGCGCTACGTGCAGGAACAGcgcg                                                                                                             >  2:2701741/1‑149 (MQ=255)
                                            cacGCGCCAGGATATCGTGCGCCTGCTGGGGCTGAACGATCCGCTGATTCAAATCAGCAGTTTTGACCGTCCGAATATTCGCTACATGCTGATGGAGAAGTTCAAACCGCTCGATCAGTTGATGCGCTACGTGCAGGAACAGCGCGGTa                                                                                                          <  2:508960/149‑1 (MQ=255)
                                            cacGCGCCAGGATATCGTGCGCCTGCTGGGGCTGAACGATCCGCTGATTCAAATCAGCAGTTTTGACCGTCCGAATATTCGCTACATGCTGATGGAGAAGTTCAAACCGCTCGATCAGTTGATGCGCTACGTGCAGGAACAGCGCGGTa                                                                                                          <  1:3869886/149‑1 (MQ=255)
                                             acGCGCCAGGATATCGTGCGCCTGCTGGGGCTGAACGATCCGCTGATTCAAATCAGCAGTTTTGACCGTCCGAATATTCGCTACATGCTGATGGAGAAGTTCAAACCGCTCGATCAGTTGATGCGCTACGTGCAGGAACAGCGCGGTaa                                                                                                         <  2:882011/149‑1 (MQ=255)
                                             acGCGCCAGGATATCGTGCGCCTGCTGGGGCTGAACGATCCGCTGATTCAAATCAGCAGTTTTGACCGTCCGAATATTCGCTACATGCTGATGGAGAAGTTCAAACCGCTCGATCAGTTGATGCGCTACGTGCAGGAACAGCGCGGTaa                                                                                                         >  1:2394906/1‑149 (MQ=255)
                                                 gcCAGGATATCGTGCGCCTGCTGGGGCTGAACGATCCGCTGATTCAAATCAGCAGTTTTGACCGTCCGAATATTCGCTACATGCTGATGGAGAAGTTCAAACCGCTCGATCAGTTGATGCGCTACGTGCAGGAACAGCGCGGTAAGTCa                                                                                                     >  1:3014575/1‑149 (MQ=255)
                                                  ccAGGATATCGTGCGCCTGCTGGGGCTGAACGATCCGCTGATTCAAATCAGCAGTTTTGACCGTCCGAATATTCGCTACATGCTGATGGAGAAGTTCAAACCGCTCGATCAGTTGATGCGCTACGTGCAGGAACAGCGCGGTAAGTCAg                                                                                                    >  2:2854129/1‑149 (MQ=255)
                                                   cAGGATATCGTGCGCCTGCTGGGGCTGAACGATCCGCTGATTCAAATCAGCAGTTTTGACCGTCCGAATATTCGCTACATGCTGATGGAGAAGTTCAAACCGCTCGATCAGTTGATGCGCTACGTGCAGGAACAGCGCGGTAAGTCAgg                                                                                                   >  2:2909652/1‑149 (MQ=255)
                                                      gATATCGTGCGCCTGCTGGGGCTGAACGATCCGCTGATTCAAATCAGCAGTTTTGACCGTCCGAATATTCGCTACATGCTGATGGAGAAGTTCAAACCGCTCGATCAGTTGATGCGCTACGTGCAGGAACAGCGCGGTAAGTCAGGCAt                                                                                                <  2:2590558/149‑1 (MQ=255)
                                                       atatCGTGCGCCTGCTGGGGCTGAACGATCCGCTGATTCAAATCAGCAGTTTTGACCGTCCGAATATTCGCTACATGCTGATGGAGAAGTTCAAACCGCTCGATCAGTTGATGCGCTACGTGCAGGAACAGCGCGGTAAGTCAGGCAtt                                                                                               <  1:4432617/149‑1 (MQ=255)
                                                        tatCGTGCGCCTGCTGGGGCTGAACGATCCGCTGATTCAAATCAGCAGTTTTGACCGTCCGAATATTCGCTACATGCTGATGGAGAAGTTCAGACCGCTCGATCAGTTGATGCGCTACGTGCAGGAACAGCGCGGTAAGTCAGGCATTa                                                                                              <  2:329839/149‑1 (MQ=255)
                                                        tatCGTGCGCCTGCTGGGGCTGAACGATCCGCTGATTCAAATCAGCAGTTTTGACCGTCCGAATATTCGCTACATGCTGATGGAGAAGTTCAAACCGCTCGATCAGTTGATGCGCTACGTGCAGGAACAGCGCGGTAAGTCAGGCATTa                                                                                              >  2:1606990/1‑149 (MQ=255)
                                                        tatCGTGCGCCTGCTGGGGCTGAACGATCCGCTGATTCAAATCAGCAGTTTTGACCGTCCGAATATTCGCGACATGCTGATGGAGAAGTTCAAACCGCTCGATCAGTTGATGCGCTACGTGCAGGAACAGCGCGGTAAGTCAGGCATTa                                                                                              <  2:2383905/149‑1 (MQ=255)
                                                         atCGTGCGCCTGCTGGGGCTGAACGATCCGCTGATTCAAATCAGCAGTTTTGACCGTCCGAATATTCGCTACATGCTGATGGAGAAGTTCAAACCGCTCGATCAGTTGATGCGCTACGTGCAGGAACAGCGCGGTAAGTCAGGCATTAt                                                                                             <  1:550041/149‑1 (MQ=255)
                                                         atCGTGCGCCTGCTGGGGCTGAACGATCCGCTGATTCAAATCAGCAGTTTTGACCGTCCGAATATTCGCTACATGCTGATGGAGAAGTTCAAACCGCTCGATCAGTTGATGCGCTACGTGCAGGAACAGCGCGGTAAGTCAGGCATTAt                                                                                             >  1:4331151/1‑149 (MQ=255)
                                                          tCGTGCGCCTGCTGGGGCTGAACGATCCGCTGATTCAAATCAGCAGTTTTGACCGTCCGAATATTCGCTACATGCTGATGGAGAAGTTCAAACCGCTCGATCAGTTGATGCGCTACGTGCAGGAACAGCGCGGTAAGTCAGGCATTATc                                                                                            <  1:3286464/149‑1 (MQ=255)
                                                            gTGCGCCTGCTGGGGCTGAACGATCCGCTGATTCAAATCAGCAGTTTTGACCGTCCGAATATTCGCTACATGCTGATGGAGAAGTTCAAACCGCTCGATCAGTTGATGCGCTACGTGCAGGAACAGCGCGGTAAGTCAGGCATTATCTa                                                                                          <  2:3275652/149‑1 (MQ=255)
                                                              gcgcCTGCTGGGGCTGAACGATCCGCTGATTCAAATCAGCAGTTTTGACCGTCCGAATATTCGCTACATGCTGATGGAGAAGTTCAAACCGCTCGATCAGTTGATGCGCTACGTGCAGGAACAGCGCGGTAAGTCAGGCATTATCTACt                                                                                        <  2:3708869/149‑1 (MQ=255)
                                                                  ctgctgGGGCTGAACGATCCGCTGATTCAAATCAGCAGTTTTGACCGTCCGAATATTCGCTACATGCTGATGGAGAAGTTCAAACCGCTCGATCAGTTGATGCGCTACGTGCAGGAACAGCGCGGTAAGTCAGGCATTATCTACTGCaa                                                                                    >  1:2891657/1‑149 (MQ=255)
                                                                  ctgctgGGGCTGAACGATCCGCTGATTCAAATCAGCAGTTTTGACCGTCCGAATATTCGCTACATGCTGATGGAGAAGTTCAAACCGCTCGATCAGTTGATGCGCTACGTGCAGGAACAGCGCGGTAAGTCAGGCATTATCTACTGCaa                                                                                    <  2:487241/149‑1 (MQ=255)
                                                                   tgctgGGGCTGAACGATCCGCTGATTCAAATCAGCAGTTTTGACCGTCCGAATATTCGCTACATGCTGATGGAGAAGTTCAAACCGCTCGATCAGTTGATGCGCTACGTGCAGGAACAGCGCGGTAAGTCAGGCATTATCTACTGCAAc                                                                                   >  2:5149388/1‑149 (MQ=255)
                                                                   tgctgGGGCTGAACGATCCGCTGATTCAAATCAGCAGTTTTGACCGTCCGAATATTCGCTACATGCTGATGGAGAAGTTCAAACCGCTCGATCAGTTGATGCGCTACGTGCAGGAACAGCGCGGTAAGTCAGGCATTATCTACTGCAAc                                                                                   >  2:4587450/1‑149 (MQ=255)
                                                                    gctgGGGCTGAACGATCCGCTGATTCAAATCAGCAGTTTTGACCGTCCGAATATTCGCTACATGCTGATGGAGAAGTTCAAACCGCTCGATCAGTTGATGCGCTACGTGCAGGAACAGCGCGGTAAGTCAGGCATTATCTACTgcagca                                                                                  <  1:2992148/149‑1 (MQ=255)
                                                                     ctgGGGCTGAACGATCCGCTGATTCAAATCAGCAGTTTTGACCGTCCGAATATTCGCTCCATGCTGATGGAGAAGTTCAAACCGCTCGATCAGTTGATGCGCTACGTGCAGGAACAGCGCGGTAAGTCAGGCATTATCTACTGCAACAg                                                                                 >  2:1175726/1‑149 (MQ=255)
                                                                       ggggCTGAACGATCCGCTGATTCAAATCAGCAGTTTTGACCGTCCGAATATTCGCTACATGCTGATGGAGAAGTTCAGACCGCTCGATCAGTTGATGCGCTACGTGCAGGAACAGCGCGGTAAGTCAGGCATTATCTACTGCAACAGcc                                                                               >  1:3952888/1‑149 (MQ=255)
                                                                          gCTGAACGATCCGCTGATTCAAATCAGCAGTTTTGACCGTCCGAATATTCGCTACATGCTGATGGAGAAGTTCAAACCGCTCGATCAGTTGATGCGCTACGTGCAGGAACAGCGCGGTAAGTCAGGCATTATCTACTGCAACAGCcgcg                                                                            >  2:5106478/1‑149 (MQ=255)
                                                                          gCTGAACGATCCGCTGATTCAAATCAGCAGTTTTGACCGTCCGAATATTCGCTACATGCTGATGGAGAAGTTCAAACCGCTCGATCAGTTGATGCGCTACGTGCAGGAACAGCGCGGTAAGTCAGGCATTATCTACTGCAACAGCcgcg                                                                            >  1:4293389/1‑149 (MQ=255)
                                                                            tGAACGATCCGCTGATTCAAATCAGCAGTTTTGACCGTCCGAATATTCGCTACATGCTGATGGAGAAGTTCAAACCGCTCGATCAGTTGATGCGCTACGTGCAGGAACAGCGCGGTAAGTCAGGCATTATCTACTGCAACAGCcgcgcg                                                                          >  1:3348199/1‑149 (MQ=255)
                                                                            tGAACGATCCGCTGATTCAAATCAGCAGTTTTGACCGTCCGAATATTCGCTACATGCTGATGGAGAAGTTCAAACCGCTCGATCAGTTGATGCGCTACGTGCAGGAACAGCGCGGTAAGTCAGGCATTATCTACTGCAACAGCcgcgcg                                                                          >  2:208230/1‑149 (MQ=255)
                                                                             gAACGATCCGCTGATTCAAATCAGCAGTTTTGACCGTCCGAATATTCGCTACATGCTGATGGAGAAGTTCAAACCGCTCGATCAGTTGATGCGCTACGTGCAGGAACAGCGCGGTAAGTCAGGCATTATCTACTGCAACAGCCGCGCGa                                                                         >  2:57788/1‑149 (MQ=255)
                                                                              aaCGATCCGCTGATTCAAATCAGCAGTTTTGACCGTCCGAATATTCGCTACATGCTGATGGAGAAGTTCAAACCGCTCGATCAGTTGATGCGCTACGTGCAGGAACAGCGCGGTAAGTCAGGCATTATCTACTGCAACAGCCGCGCGaa                                                                        >  1:128843/1‑149 (MQ=255)
                                                                                cGATCCGCTGATTCAAATCAGCAGTTTTGACCGTCCGAATATTCGCTACATGCTGATGGAGAAGTTCAAACCGCTCGATCAGTTGATGCGCTACGTGCAGGAACAGCGCGGTAAGTCAGGCATTATCTACTGCAACAGCCGCGCGAAAg                                                                      <  1:4848479/149‑1 (MQ=255)
                                                                                  aTCCGCTGATTCAAATCAGCAGTTTTGACCGTCCGAATATTCGCTACATGCTGATGGAGAAGTTCAAACCGCTCGATCAGTTGATGCGCTACGTGCAGGAACAGCGCGGTAAGTCAGGCATTATCTACTGCAACAGCCGCGCGAAAGTa                                                                    <  1:954703/149‑1 (MQ=255)
                                                                                  aTCCGCTGATTCAAATCAGCAGTTTTGACCGTCCGAATATTCGCTACATGCTGATGGAGAAGTTCAAACCGCTCGATCAGTTGATGCGCTACGTGCAGGAACAGCGCGGTAAGTCAGGCATTATCTACTGCAACAGCCGCGCGAAAGTa                                                                    <  2:2500626/149‑1 (MQ=255)
                                                                                   tCCGCTGATTCAAATCAGCAGTTTTGACCGTCCGAATATTCGCTACATGCTGATGGAGAAGTTCAAACCGCTCGATCAGTTGATGCGCTACGTGCAGGAACAGCGCGGTAAGTCAGGCATTATCTACTGCAACAGCCGCGCGAAAGTag                                                                   >  1:3157172/1‑149 (MQ=255)
                                                                                      gCTGATTCAAATCAGCAGTTTTGACCGTCCGAATATTCGCTACATGCTGATGGAGAAGTTCAAACCGCTCGATCAGTTGATGCGCTACGTGCAGGAACAGCGCGGTAAGTCAGGCATTATCTACTGCAACAGCCGCGCGAAAGTagaag                                                                <  2:2053516/149‑1 (MQ=255)
                                                                                       cTGATTCAAATCAGCAGTTTTGACCGTCCGAATATTCGCTACATGCTGATGGAGAAGTTCAAACCGCTCGATCAGTTGATGCGCTACGTGCAGGAACAGCGCGGTAAGTCAGGCATTATCTACTGCAACAGCCGCGCGAAAGTagaaga                                                               <  2:10986/149‑1 (MQ=255)
                                                                                        tGATTCAAATCAGCAGTTTTGACCGTCCGAATATTCGCTACATGCTGATGGAGAAGTTCAAACCGCTCGATCAGTTGATGCGCTACGTGCAGGAACAGCGCGGTAAGTCAGGCATTATCTACTGCAACAGCCGCGCGAAAGTAGAAGac                                                              <  1:1969059/149‑1 (MQ=255)
                                                                                         gATTCAAATCAGCAGTTTTGACCGTCCGAATATTCGCTACATGCTGATGGAGAAGTTCAAACCGCTCGATCAGTTGATGCGCTACGTGCAGGAACAGCGCGGTAAGTCAGGCATTATCTACTGCAACAGCCGCGCGAAAGTAGAAGaca                                                             >  1:710640/1‑149 (MQ=255)
                                                                                          aTTCAAATCAGCAGTTTTGACCGTCCGAATATTCGCTACATGCTGATGGAGAAGTTCAAACCGCTCGATCAGTTGATGCGCTACGTGCAGGAACAGCGCGGTAAGTCAGGCATTATCTACTGCAACAGCCGCGCGAAAGTAGAAGacac                                                            >  2:2171063/1‑149 (MQ=255)
                                                                                          aTTCAAATCAGCAGTTTTGACCGTCCGAATATTCGCTACATGCTGATGGAGAAGTTCAAACCGCTCGATCAGTTGATGCGCTACGTGCAGGAACAGCGCGGTAAGTCAGGCATTATCTACTGCAACAGCCGCGCGAAAGTAGAAGacac                                                            <  1:4710075/149‑1 (MQ=255)
                                                                                               aaTCAGCAGTTTTGACCGTCCGAATATTCGCTACATGCTGATGGAGAAGTTCAAACCGCTCGATCAGTTGATGCGCTACGTGCAGGAACAGCGCGGTAAGTCAGGCATTATCTACTGCAACAGCCGCGCGAAAGTAGAAGACACCGCTg                                                       <  1:1311351/149‑1 (MQ=255)
                                                                                               aaTCAGCAGTTTTGACCGTCCGAATATTCGCTACATGCTGATGGAGAAGTTCAAACCGCTCGATCAGTTGATGCGCTACGTGCAGGAACAGCGCGGTAAGTCAGGCATTATCTACTGCAACAGCCGCGCGAAAGTAGAAGACACCGCTg                                                       >  1:3728912/1‑149 (MQ=255)
                                                                                               aaTCAGCAGTTTTGACCGTCCGAATATTCGCTACATGCTGATGGAGAAGTTCAAACCGCTCGATCAGTTGATGCGCTACGTGCAGGAACAGCGCGGTAAGTCAGGCATTATCTACTGCAACAGCCGCGCGAAAGTAGAAGACACCGCTg                                                       >  1:1133222/1‑149 (MQ=255)
                                                                                                aTCAGCAGTTTTGACCGTCCGAATATTCGCTACATGCTGATGGAGAAGTTCAAACCGCTCGATCAGTTGATGCGCTACGTGCAGGAACAGCGCGGTAAGTCAGGCATTATCTACTGCAACAGCCGCGCGAAAGTAGAAGACACCGCTgc                                                      >  2:2277895/1‑149 (MQ=255)
                                                                                                aTCAGCAGTTTTGACCGTCCGAATATTCGCTACATGCTGATGGAGAAGTTCAAACCGCTCGATCAGTTGATGCGCTACGTGCAGGAACAGCGCGGTAAGTCAGGCATTATCTACTGCAACAGCCGCGCGAAAGTAGAAGACACCGCTgc                                                      >  2:3463240/1‑149 (MQ=255)
                                                                                                 tCAGCAGTTTTGACCGTCCGAATATTCGCTACATGCTGATGGAGAAGTTCAAACCGCTCGATCAGTTGATGCGCTACGTGCAGGAACAGCGCGGTAAGTCAGGCATTATCTACTGCAACAGCCGCGCGAAAGTAGAAGACACCGCTgcg                                                     >  1:2819473/1‑149 (MQ=255)
                                                                                                    gcagTTTTGACCGTCCGAATATTCGCTACATGCTGATGGAGAAGTTCAAACCGCTCGATCAGTTGATGCGCTACGTGCAGGAACAGCGCGGTAAGTCAGGCATTATCTACTGCAACAGCCGCGCGAAAGTAGAAGACACCGCTgcgcgc                                                  >  1:2210150/1‑149 (MQ=255)
                                                                                                      agTTTTGACCGTCCGAATATTCGCTACATGCTGATGGAGAAGTTCAAACCGCTCGATCAGTTGATGCGCTACGTGCAGGAACAGCGCGGTAAGTCAGGCATTATCTACTGCAACAGCCGCGCGAAAGTAGAAGACACCGCTGCGCGCCt                                                <  1:3779644/149‑1 (MQ=255)
                                                                                                            gACCGTCCGAATATTCGCTACATGCTGATGGAGAAGTTCAAACCGCTCGATCAGTTGATGCGCTACGTGCAGGAACAGCGCGGTAAGTCAGGCATTATCTACTGCAACAGCCGCGCGAAAGTAGAAGACACCGCTGCGCGCCTGCAAAg                                          >  1:5037401/1‑149 (MQ=255)
                                                                                                                gTCCGAATATTCGCTACATGCTGATGGAGAAGTTCAAACCGCTCGATCAGTTGATGCGCTACGTGCAGGAACAGCGCGGTAAGTCAGGCATTATCTACTGCAACAGCCGCGCGAAAGTAGAAGACACCGCTGCGCGCCTGCAAAGCAAg                                      >  2:26446/1‑149 (MQ=255)
                                                                                                                gTCCGAATATTCGCTACATGCTGATGGAGAAGTTCAAACCGCTCGATCAGTTGATGCGCTACGTGCAGGAACAGCGCGGTAAGTCAGGCATTATCTACTGCAACAGCCGCGCGAAAGTAGAAGACACCGCTGCGCGCCTGCAAAGCAAg                                      >  1:2685927/1‑149 (MQ=255)
                                                                                                                 tCCGAATATTCGCTACATGCTGATGGAGAAGTTCAAACCGCTCGATCAGTTGATGCGCTACGTGCAGGAACAGCGCGGTAAGTCAGGCATTATCTACTGCAACAGCCGCGCGAAAGTAGAAGACACCGCTGCGCGCCTGCAAAGCAAgg                                     >  1:2583792/1‑149 (MQ=255)
                                                                                                                    gAATATTCGCTACATGCTGATGGAGAAGTTCAAACCGCTCGATCAGTTGATGCGCTACGTGCAGGAACAGCGCGGTAAGTCAGGCATTATCTACTGCAACAGCCGCGCGAAAGTAGAAGACACCGCTGCGCGCCTGCAAAGCAAGGGaa                                  <  2:3137566/149‑1 (MQ=255)
                                                                                                                    gAATATTCGCTACATGCTGATGGAGAAGTTCAAACCGCTCGATCAGTTGATGCGCTACGTGCAGGAACAGCGCGGTAAGTCAGGCATTATCTACTGCAACAGCCGCGCGAAAGTAGAAGACACCGCTGCGCGCCTGCAAAGCAAGGGaa                                  >  2:3903749/1‑149 (MQ=255)
                                                                                                                     aaTATTCGCTACATGCTGATGGAGAAGTTCAAACCGCTCGATCAGTTGATGCGCTACGTGCAGGAACAGCGCGGTAAGTCAGGCATTATCTACTGCAACAGCCGCGCGAAAGTAGAAGACACCGCTGCGCGCCTGCAAAGCAAGGGAAt                                 >  2:3238372/1‑149 (MQ=255)
                                                                                                                     aaTATTCGCTACATGCTGATGGAGAAGTTCAAACCGCTCGATCAGTTGATGCGCTACGTGCAGGAACAGCGCGGTAAGTCAGGCATTATCTACTGCAACAGCCGCGCGAAAGTAGAAGACACCGCTGCGCGCCTGCAAAGCAAGGGAAt                                 >  1:1698221/1‑149 (MQ=255)
                                                                                                                      atatTCGCTACATGCTGATGGAGAAGTTCAAACCGCTCGATCAGTTGATGCGCTACGTGCAGGAACAGCGCGGTAAGTCAGGCATTATCTACTGCAACAGCCGCGCGAAAGTAGAAGACACCGCTGCGCGCCTGCAAAGCAAGGGAAtt                                <  2:2977531/149‑1 (MQ=255)
                                                                                                                       tatTCGCTACATGCTGATGGAGAAGTTCAAACCGCTCGATCAGTTGATGCGCTACGTGCAGGAACAGCGCGGTAAGTCAGGCATTATCTACTGCAACAGCCGCGCGAAAGTAGAAGACACCGCTGCGCGCCTGCAAAGCAAGGGAATTa                               <  1:2707667/149‑1 (MQ=255)
                                                                                                                         ttCGCTACATGCTGATGGAGAAGTTCAAACCGCTCGATCAGTTGATGCGCTACGTGCAGGAACAGCGCGGTAAGTCAGGCATTATCTACTGCAACAGCCGCGCGAAAGTAGAAGACACCGCTGCGCGCCTGCAAAGCAAGGGAATTAgc                             <  1:4661184/149‑1 (MQ=255)
                                                                                                                         ttCGCTACATGCTGATGGAGAAGTTCAAACCGCTCGATCAGTTGATGCGCTACGTGCAGGAACAGCGCGGTAAGTCAGGCATTATCTACTGCAACAGCCGCGCGAAAGTAGAAGACACCGCTGCGCGCCTGCAAAGCAAGGGAATTAgc                             <  1:620991/149‑1 (MQ=255)
                                                                                                                           cGCTACATGCTGATGGAGAAGTTCAAACCGCTCGATCAGTTGATGCGCTACGTGCAGGAACAGCGCGGTAAGTCAGGCATTATCTACTGCAACAGCCGCGCGAAAGTAGAAGACACCGCTGCGCGCCTGCAAAGCAAGGGAATTAgcgc                           >  1:461889/1‑149 (MQ=255)
                                                                                                                              tACATGCTGATGGAGAAGTTCAAACCGCTCGATCAGTTGATGCGCTACGTGCAGGAACAGCGCGGTAAGTCAGGCATTATCTACTGCAACAGCCGCGCGAAAGTGGAAGACACCGCTGCGCGCCTGCAAAGCAAGGGAATTAGCgcggc                        <  2:3753717/149‑1 (MQ=255)
                                                                                                                                 aTGCTGATGGAGAAGTTCAAACCGCTCGATCAGTTGATGCGCTACGTGCAGGAACAGCGCGGTAAGTCAGGCATTATCTACTGCAACAGCCGCGCGAAAGTAGAAGACACCGCTGCGCGCCTGCAAAGCAAGGGAATTAGCGCggcggc                     >  2:1213739/1‑149 (MQ=255)
                                                                                                                                    cTGATGGAGAAGTTCAAACCGCTCGATCAGTTGATGCGCTACGTGCAGGAACAGCGCGGTAAGTCAGGCATTATCTACTGCAACAGCCGCGCGAAAGTAGAAGACACCGCTGCGCGCCTGCAAAGCAAGGGAATTAGCGCGGCGGCCTa                  >  1:2063177/1‑149 (MQ=255)
                                                                                                                                    cTGATGGAGAAGTTCAAACCGCTCGATCAGTTGATGCGCTACGTGCAGGAACAGCGCGGTAAGTCAGGCATTATCTACTGCAACAGCCGCGCGAAAGTAGAAGACACCGCTGCGCGCCTGCAAAGCAAGGGAATTAGCGCGGCGGCCTa                  >  1:2935516/1‑149 (MQ=255)
                                                                                                                                     tGATGGAGAAGTTCAAACCGCTCGATCAGTTGATGCGCTACGTGCAGGAACAGCGCGGTAAGTCAGGCATTATCTACTGCAACAGCCGCGCGAAAGTAGAAGACACCGCTGCGCGCCTGCAAAGCAAGGGAATTAGCGCGGCGGCCTAt                 <  1:4157885/149‑1 (MQ=255)
                                                                                                                                          gagaAGTTCAAACCGCTCGATCAGTTGATGCGCTACGTGCAGGAACAGCGCGGTAAGTCAGGCATTATCTACTGCAACAGCCGCGCGAAAGTAGAAGACACCGCTGCGCGCCTGCAAAGCAAGGGAATTAGCGCGGCGGCCTATCATGc            <  1:2072163/149‑1 (MQ=255)
                                                                                                                                           agaAGTTCAAACCGCTCGATCAGTTGATGCGCTACGTGCAGGAACAGCGCGGTAAGTCAGGCATTATCTACTGCAACAGCCGCGCGAAAGTAGAAGACACCGCTGCGCGCCTGCAAAGCAAGGGAATTAGCGCGGCGGCCTATCATGcc           >  1:651774/1‑149 (MQ=255)
                                                                                                                                           agaAGTTCAAACCGCTCGATCAGTTGATGCGCTACGTGCAGGAACAGCGCGGTAAGTCAGGCATTATCTACTGCAACAGCCGCGCGAAAGTAGAAGACACCGCTGCGCGCCTGCAAAGCAAGGGAATTAGCGCGGCGGCCTATCATGcc           >  2:1693831/1‑149 (MQ=255)
                                                                                                                                               gTTCAGACCGCTCGATCAGTTGATGCGCTACGTGCAGGAACAGCGCGGTAAGTCAGGCATTATCTACTGCAACAGCCGCGCGAAAGTAGAAGACACCGCTGCGCGCCTGCAAAGCAAGGGAATTAGCGCGGCGGCCTATCATGCCGGGc       <  2:246020/149‑1 (MQ=255)
                                                                                                                                               gTTCAAACCGCTCGATCAGTTGATGCGCTACGTGCAGGAACAGCGCGGTAAGTCAGGCATTATCTACTGCAACAGCCGCGCGAAAGTAGAAGACACCGCTGCGCGCCTGCAAAGCAAGGGAATTAGCGCGGCGGCCTATCATGCCGGGc       <  1:2402787/149‑1 (MQ=255)
                                                                                                                                               gTTCAAACCGCTCGATCAGTTGATGCGCTACGTGCAGGAACAGCGCGGTAAGTCAGGCATTATCTACTGCAACAGCCGCGCGAAAGTAGAAGACACCGCTGCGCGCCTGCAAAGCAAGGGAATTAGCGCGGCGGCCTATCATGCCGGGc       >  2:4890570/1‑149 (MQ=255)
                                                                                                                                                    aaCCGCTCGATCAGTTGATGCGCTACGTGCAGGAACAGCGCGGTAAGTCAGGCATTATCTACTGCAACAGCCGCGCGAAAGTAGAAGACACCGCTGCGCGCCTGCAAAGCAAGGGAATTAGCGCGGCGGCCTATCATGCCGGGCTGGaa  >  1:992328/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CCGACGCTGCCGTTTATGGCGCTGACCGCCACAGCCGACGACACCACGCGCCAGGATATCGTGCGCCTGCTGGGGCTGAACGATCCGCTGATTCAAATCAGCAGTTTTGACCGTCCGAATATTCGCTACATGCTGATGGAGAAGTTCAAACCGCTCGATCAGTTGATGCGCTACGTGCAGGAACAGCGCGGTAAGTCAGGCATTATCTACTGCAACAGCCGCGCGAAAGTAGAAGACACCGCTGCGCGCCTGCAAAGCAAGGGAATTAGCGCGGCGGCCTATCATGCCGGGCTGGAA  >  CP009273/3999740‑4000036

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: