Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 456,298 T→A 100% L652Q (CTG→CAG)  lon → DNA‑binding ATP‑dependent protease La

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273456,2980TA97.5% 1111.8 / 11.7 320L652Q (CTG→CAG) lonDNA‑binding ATP‑dependent protease La
Reads supporting (aligned to +/- strand):  ref base T (2/6);  new base A (131/181);  total (133/187)
Fisher's exact test for biased strand distribution p-value = 4.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 320 total aligned reads displayed.

GGACGGAAGTGGGCGGTGACTTGCTGACCATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCGATG  >  CP009273/456151‑456444
                                                                                                                                                   |                                                                                                                                                  
ggACGGAAGTGGGCGGTGACTTGCTGACCATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAgg                                                                                                                                                  <  2:2763474/150‑1 (MQ=255)
  aCGGAAGTGGGCGGTGACTTGCTGACCATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAgggg                                                                                                                                                <  1:839664/150‑1 (MQ=255)
   cGGAAGTGGGCGGTGACTTGCTGACCATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGa                                                                                                                                               <  1:1532047/150‑1 (MQ=255)
     gAAGTGGGCGGTGACTTGCTGACCATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATc                                                                                                                                             <  2:2696925/150‑1 (MQ=255)
     gAAGTGGGCGGTGACTTGCTGACCATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATc                                                                                                                                             <  2:2623745/150‑1 (MQ=255)
      taGTGGGCGGTGACTTGCTGACCTTTGAAACCGCATGTGGTCCGGGTAAAGGCTATCTGACCTATTCAGGTTCGCTCGGCGAAGTGATGCAGGTGTCCATTCAGGCGGCGTTAACGGGGGTTCGTGCGCGTGCGGAAAAACAGGGGATCa                                                                                                                                            <  2:682218/149‑1 (MQ=255)
        gTGGGCGGTGACTTGCTGACCATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCTGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAAc                                                                                                                                          <  1:3255379/150‑1 (MQ=255)
        gTGGGCGGTGACTTGCTGACCATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAAc                                                                                                                                          <  2:772276/150‑1 (MQ=255)
         tGGGCGGTGACTTGCTGACCATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAAcc                                                                                                                                         <  1:2213443/150‑1 (MQ=255)
          gggCGGTGACTTGCTGACCATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAAccc                                                                                                                                        <  1:2049316/150‑1 (MQ=255)
           ggCGGTGACGTGCGGACCATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGACAAACAGGGGATCAACCCt                                                                                                                                       <  2:4517556/150‑1 (MQ=255)
            gCGGTGACTTGCTGACCATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTg                                                                                                                                      <  2:1180076/150‑1 (MQ=255)
             cGGTGACTTGCTGACCATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGa                                                                                                                                     <  1:4544095/150‑1 (MQ=255)
                tGACTTGCTGACCATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGAttt                                                                                                                                  >  2:841171/1‑150 (MQ=255)
                tGACTTGCTGACCATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGAttt                                                                                                                                  >  2:4321356/1‑150 (MQ=255)
                   cTTGCTGACCATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTa                                                                                                                               <  1:4395971/150‑1 (MQ=255)
                   cTTGCTGACCATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTa                                                                                                                               <  2:358825/150‑1 (MQ=255)
                    ttGCTGACCATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTAc                                                                                                                              <  2:407444/150‑1 (MQ=255)
                     tGCTGACCATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACg                                                                                                                             >  1:467664/1‑150 (MQ=255)
                     tGCTGACCATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACg                                                                                                                             >  1:1871007/1‑150 (MQ=255)
                          aCCATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGaaaaa                                                                                                                        >  2:4619496/1‑150 (MQ=255)
                           ccATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAAc                                                                                                                       <  1:400214/150‑1 (MQ=255)
                              ttGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGTAGTTATGCAGGAGTCCATTCAGACGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACATGGGATCTACCCTGATTTTTACGAAAAACGTg                                                                                                                    <  2:3031501/150‑1 (MQ=255)
                               tGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGa                                                                                                                   >  1:4241458/1‑150 (MQ=255)
                               tGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGa                                                                                                                   >  1:4151681/1‑150 (MQ=255)
                                 aaaCCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACa                                                                                                                 >  2:692298/1‑150 (MQ=255)
                                 aaaCCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACa                                                                                                                 >  2:1900086/1‑150 (MQ=255)
                                   aCCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACAtc                                                                                                               >  2:3466856/1‑150 (MQ=255)
                                   aCCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACAtc                                                                                                               >  2:2417516/1‑150 (MQ=255)
                                   aCCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACAtc                                                                                                               <  2:3823382/150‑1 (MQ=255)
                                      gCATGTGTTCCGGGTAAAGGCAAACTGTCCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACAtccac                                                                                                            <  1:2030250/150‑1 (MQ=255)
                                      gCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACAtccac                                                                                                            >  1:3968040/1‑150 (MQ=255)
                                      gCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACAtccac                                                                                                            >  2:4394546/1‑150 (MQ=255)
                                       cATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACAtccacg                                                                                                           >  2:3667233/1‑150 (MQ=255)
                                          gtgtTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACAtccacgtcc                                                                                                        <  1:2846264/150‑1 (MQ=255)
                                           tgtTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACAtccacgtcca                                                                                                       <  1:128806/150‑1 (MQ=255)
                                            gtTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACAtccacgtccac                                                                                                      >  2:2166884/1‑150 (MQ=255)
                                            gtTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACAtccacgtccac                                                                                                      >  1:1574804/1‑150 (MQ=255)
                                              tCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATccacgtccacgt                                                                                                    >  1:97455/1‑150 (MQ=255)
                                                    tAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGa                                                                                              <  2:256438/150‑1 (MQ=255)
                                                    tAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGa                                                                                              >  2:2367761/1‑150 (MQ=255)
                                                    tAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGa                                                                                              >  2:1273928/1‑150 (MQ=255)
                                                       aGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAgg                                                                                           >  1:1269708/1‑150 (MQ=255)
                                                          cAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGc                                                                                        <  2:3515463/150‑1 (MQ=255)
                                                           aaaCTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCg                                                                                       >  2:148080/1‑150 (MQ=255)
                                                             aCTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGAc                                                                                     <  1:4206948/150‑1 (MQ=255)
                                                                 aCCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCg                                                                                 <  1:312128/150‑1 (MQ=255)
                                                                   cTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGaa                                                                               >  1:2626431/1‑150 (MQ=255)
                                                                    tataCCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGaaa                                                                              >  1:3900460/1‑150 (MQ=255)
                                                                    tataCCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGaaa                                                                              <  1:3801519/150‑1 (MQ=255)
                                                                     ataCCGGTGCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCTCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGTAAAACGTGACATCCACGTCCTCGTACCGGAAGGTGCGACGCCGAAAg                                                                             <  1:2498372/150‑1 (MQ=255)
                                                                       aCCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGAt                                                                           <  1:4574692/150‑1 (MQ=255)
                                                                       aCCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGAt                                                                           <  1:2696117/150‑1 (MQ=255)
                                                                          ggTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGt                                                                        <  2:3444891/150‑1 (MQ=255)
                                                                          ggTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGt                                                                        >  1:592314/1‑150 (MQ=255)
                                                                          ggTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGt                                                                        >  1:207877/1‑150 (MQ=255)
                                                                          ggTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGt                                                                        <  2:2332962/150‑1 (MQ=255)
                                                                           gTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTc                                                                       <  1:1065501/150‑1 (MQ=255)
                                                                              cGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGa                                                                    <  1:1686429/150‑1 (MQ=255)
                                                                                cTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGt                                                                  <  1:2164227/150‑1 (MQ=255)
                                                                                    gCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCg                                                              <  1:3343051/150‑1 (MQ=255)
                                                                                       aaGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTa                                                           <  1:3648581/150‑1 (MQ=255)
                                                                                         gTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTAtt                                                         <  1:1749358/150‑1 (MQ=255)
                                                                                              gCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTAt                                                    <  2:3729523/150‑1 (MQ=255)
                                                                                               cAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTAtg                                                   <  1:3416668/150‑1 (MQ=255)
                                                                                               cAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGACGAAAGATGGTCCGAGTGCCGGTATTGCTAtg                                                   <  1:4233399/150‑1 (MQ=255)
                                                                                                 ggAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTAtgtg                                                 <  2:3572037/150‑1 (MQ=255)
                                                                                                  gAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGc                                                >  2:2996133/1‑150 (MQ=255)
                                                                                                     tCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCAcc                                             <  1:3548879/150‑1 (MQ=255)
                                                                                                       cATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGAAATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACcgc                                           >  2:292965/1‑150 (MQ=255)
                                                                                                             ggcggcGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGt                                     <  1:3981618/150‑1 (MQ=255)
                                                                                                               cggcgTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGttt                                   <  1:2457834/150‑1 (MQ=255)
                                                                                                                  cgTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTGTCtt                                >  1:1930242/1‑150 (MQ=255)
                                                                                                                   gTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTg                               <  2:2029283/150‑1 (MQ=255)
                                                                                                                       actgtggtTCGTGCGAGTGCGGAAAAACAGGGGATCAACACTGATTTTAACTAAAAACGAGACATCCACGTCCACTTACCGGAAGGTGAGACGACGAAAGATGGTCCTAGTGACGGTATTGATATGTGCACCGCGCTGGTTTCTTGCCTg                           >  2:4454122/4‑150 (MQ=255)
                                                                                                                       aCGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTg                           <  2:855389/150‑1 (MQ=255)
                                                                                                                        cGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGa                          >  1:735939/1‑150 (MQ=255)
                                                                                                                         ggtggtTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGAc                         >  2:2549699/1‑150 (MQ=255)
                                                                                                                          gtggtTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGATGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGAcc                        <  2:2459872/150‑1 (MQ=255)
                                                                                                                          gtggtTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGAcc                        <  2:4108393/150‑1 (MQ=255)
                                                                                                                           tggtTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCg                       >  2:2251023/1‑150 (MQ=255)
                                                                                                                            ggtTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCgg                      >  1:1353320/1‑150 (MQ=255)
                                                                                                                             gtTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGt                     <  2:294623/150‑1 (MQ=255)
                                                                                                                                cgtgcgcgtgcgGAAAAACAGGGGAACAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAAc                  <  1:2424525/150‑1 (MQ=255)
                                                                                                                                 gtgcgcgtgcgGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAAcc                 >  2:2084624/1‑150 (MQ=255)
                                                                                                                                    cgcgtgcgGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCgg              >  2:3151070/1‑150 (MQ=255)
                                                                                                                                    cgcgtgcgGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCgg              >  1:3748154/1‑150 (MQ=255)
                                                                                                                                       gtgcgGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTc           >  1:2675741/1‑150 (MQ=255)
                                                                                                                                       gtgcgGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTc           <  2:1344007/150‑1 (MQ=255)
                                                                                                                                       gtgcgGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTc           <  2:3330624/150‑1 (MQ=255)
                                                                                                                                        tgcgGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCg          <  2:2861360/150‑1 (MQ=255)
                                                                                                                                          cgGAAAAACAGGTGATCAACCCTGATTTTTACGAAAAACGTTACATCCACGTCCACGTACCGGAAGGTGCGACGCCGTAAGATGGTCAGAGTGCCGGTATTGATATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTg        <  2:504360/150‑1 (MQ=255)
                                                                                                                                          cgGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTg        <  2:1933010/150‑1 (MQ=255)
                                                                                                                                          cgGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTg        >  1:1673006/1‑150 (MQ=255)
                                                                                                                                           ggAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGc       <  1:3216086/150‑1 (MQ=255)
                                                                                                                                             aaaaaCAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCg     >  1:3976170/1‑150 (MQ=255)
                                                                                                                                             aaaaaCAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCg     <  2:3456609/150‑1 (MQ=255)
                                                                                                                                              aaaaCAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCGa    <  1:4017294/150‑1 (MQ=255)
                                                                                                                                               aaaCAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCGAt   >  2:399849/1‑150 (MQ=255)
                                                                                                                                               aaaCAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCGAt   <  1:1057022/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GGACGGAAGTGGGCGGTGACTTGCTGACCATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCGATG  >  CP009273/456151‑456444

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