Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,882,392 C→A 18.7% W356L (TGG→TTG)  cysI ← sulfite reductase, beta subunit, NAD(P)‑binding, heme‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,882,3920CA18.7% 519.7 / 72.2 189W356L (TGG→TTG) cysIsulfite reductase, beta subunit, NAD(P)‑binding, heme‑binding
Reads supporting (aligned to +/- strand):  ref base C (71/82);  new base A (23/13);  total (94/95)
Fisher's exact test for biased strand distribution p-value = 6.60e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 189 total aligned reads displayed.

CGGTACACCGGCGATGATCAGATTCTGGTTCGCCGTAATGCGGAAATCGCCTTTGTGGATCTTCGCGATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATT  >  CP009273/2882247‑2882536
                                                                                                                                                 |                                                                                                                                                
   tACACCGGCGATGATCAGATTCTGGTTCGCCGTAATGCGGAAATCGCCTTTGTGGATCTTCGCGATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTa                                                                                                                                             >  2:324263/1‑148 (MQ=255)
      acCGGCGATGATCAGATTCTGGTTCGCCGTAATGCGGAAATCGCCTTTGTGGATCTTCGCGATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTatca                                                                                                                                          >  1:3062521/1‑148 (MQ=255)
      acCGGCGATGATCAGATTCTGGTTCGCCGTAATGCGGAAATCGCCTTTGTGGATCTTCGCGATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTatca                                                                                                                                          >  1:1287216/1‑148 (MQ=255)
       ccGGCGATGATCAGATTCTGGTTCGCCGTAATGCGGAAATCGCCTTTGTGGATCTTCGCGATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTatcat                                                                                                                                         <  2:1798483/148‑1 (MQ=255)
          gCGATGATCAGATTCTGGTTCGCCGTAATGCGGAAATCGCCTTTGTGGATCTTCGCGATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCaa                                                                                                                                      <  2:2325388/148‑1 (MQ=255)
           cGATGATCAGATTCTGGTTCGCCGTAATGCGGAAATCGCCTTTGTGGATCTTCGCGATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAAt                                                                                                                                     <  2:3165361/148‑1 (MQ=255)
               gatCAGATTCTGGTTCGCCGTAATGCGGAAATCGCCTTTGTGGATCTTCGCGATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGccc                                                                                                                                 <  2:1656947/148‑1 (MQ=255)
                 tCAGATTCTGGTTCGCCGTAATGCGGAAATCGCCTTTGTGGATCTTCGCGATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCAAGTTATCATCAATGCCCtt                                                                                                                               >  2:1934433/1‑148 (MQ=255)
                   aGATTCTGGTTCGCCGTAATGCGGAAATCGCCTTTGTGGATCTTCGCGATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTaa                                                                                                                             <  1:3103270/148‑1 (MQ=255)
                   aGATTCTGGTTCGCCGTAATGCGGAAATCGCCTTTGTGGATCTTCGCGATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTaa                                                                                                                             <  1:1581731/148‑1 (MQ=255)
                     aTTCTGGTTCGCCGTAATGCGGAAATCGCCTTTGTGGATCTTCGCGATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAAcc                                                                                                                           <  1:1910537/148‑1 (MQ=255)
                         tGGTTCGCCGTAATGCGGAAATCGCCTTTGTGGATCTTCGCGATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGc                                                                                                                       <  2:2818372/148‑1 (MQ=255)
                           gTTCGCCGTAATGCGGAAATCGCCTTTGTGGATCTTCGCGATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCa                                                                                                                     <  2:38667/148‑1 (MQ=255)
                           gTTCGCCGTAATGCGGAAATCGCCTTTGTGGATCTTCGCGATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCa                                                                                                                     <  2:38653/148‑1 (MQ=255)
                              cGCCGTAATGCGGAAATCGCCTTTGTGGATCTTCGCGATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATa                                                                                                                  <  1:1578373/148‑1 (MQ=255)
                              cGCCGTAATGCGGAAATCGCCTTTGTGGATCTTCGCGATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATa                                                                                                                  <  1:3347220/148‑1 (MQ=255)
                                ccGTAATGCGGAAATCGCCTTTGTGGATCTTCGCGATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACg                                                                                                                >  1:1624355/1‑148 (MQ=255)
                                ccGTAATGCGGAAATCGCCTTTGTGGATCTTCGCGATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACg                                                                                                                <  1:2360435/148‑1 (MQ=255)
                                 cGTAATGCGGAAATCGCCTTTGTGGATCTTCGCGATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGa                                                                                                               <  1:526103/148‑1 (MQ=255)
                                 cGTAATGCGGAAATCGCCTTTGTGGATCTTCGCGATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCAAGTTATCATCAATGCCCTTAACCCAGCCAATACGa                                                                                                               <  2:1152889/148‑1 (MQ=255)
                                  gTAATGCGGAAATCGCCTTTGTGGATCTTCGCGATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGAt                                                                                                              >  2:3023404/1‑148 (MQ=255)
                                       gCGGAAATCGCCTTTGTGGATCTTCGCGATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCt                                                                                                         >  2:2457647/1‑148 (MQ=255)
                                          gAAATCGCCTTTGTGGATCTTCGCGATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGt                                                                                                      <  1:2099077/148‑1 (MQ=255)
                                           aaaTCGCCTTTGTGGATCTTCGCGATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTc                                                                                                     >  1:3200932/1‑148 (MQ=255)
                                           aaaTCGCCTTTGTGGATCTTCGCGATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTc                                                                                                     >  1:2376070/1‑148 (MQ=255)
                                           aaaTCGCCTTTGTGGATCTTCGCGATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTc                                                                                                     >  1:2375978/1‑148 (MQ=255)
                                             aTCGCCTTTGTGGATCTTCGCGATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCg                                                                                                   >  2:859360/1‑148 (MQ=255)
                                              tCGCCTTTGTGGATCTTCGCGATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCgg                                                                                                  <  2:1755537/148‑1 (MQ=255)
                                               cGCCTTTGTGGATCTTCGCGATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCAAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGt                                                                                                 >  1:1904786/1‑148 (MQ=255)
                                                   tttGTGGATCTTCGCGATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAAc                                                                                             >  1:2178204/1‑148 (MQ=255)
                                                          aTCTTCGCGATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATg                                                                                      <  1:1712330/148‑1 (MQ=255)
                                                            cTTCGCGATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGa                                                                                    >  2:3516684/1‑148 (MQ=255)
                                                            cTTCGCGATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGa                                                                                    >  2:1204072/1‑148 (MQ=255)
                                                            cTTCGCGATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGa                                                                                    <  2:3145591/148‑1 (MQ=255)
                                                             ttCGCGATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCAAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGAc                                                                                   <  1:2925757/148‑1 (MQ=255)
                                                              tCGCGATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCAAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACg                                                                                  <  2:3424353/148‑1 (MQ=255)
                                                                 cgATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGAt                                                                               <  2:2135621/148‑1 (MQ=255)
                                                                  gATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCAAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATc                                                                              >  1:1054806/1‑148 (MQ=255)
                                                                      tcCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGtt                                                                          >  1:1973994/1‑148 (MQ=255)
                                                                      tcCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGtt                                                                          <  1:324263/148‑1 (MQ=255)
                                                                           cagGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAAt                                                                     >  1:1622224/1‑148 (MQ=255)
                                                                              gCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTg                                                                  <  2:1659644/148‑1 (MQ=255)
                                                                              gCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTg                                                                  <  2:1277047/148‑1 (MQ=255)
                                                                              gCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCAAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTg                                                                  <  1:1466288/148‑1 (MQ=255)
                                                                               ccGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGa                                                                 <  1:2593712/148‑1 (MQ=255)
                                                                               ccGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGa                                                                 <  1:383499/148‑1 (MQ=255)
                                                                               ccGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGa                                                                 <  1:2381999/148‑1 (MQ=255)
                                                                               ccGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGa                                                                 <  1:662421/148‑1 (MQ=255)
                                                                                cGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGAt                                                                <  1:1017914/148‑1 (MQ=255)
                                                                                 ggTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATc                                                               <  1:1611815/148‑1 (MQ=255)
                                                                                  gTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATcc                                                              <  1:1956060/148‑1 (MQ=255)
                                                                                   ttttCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATccc                                                             <  1:1682282/148‑1 (MQ=255)
                                                                                   ttttCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATccc                                                             <  1:1682301/148‑1 (MQ=255)
                                                                                   ttttCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCACATGGACGGATCGGTTCAAATTTGATccc                                                             >  2:2275913/1‑148 (MQ=255)
                                                                                   ttttCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCAAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATccc                                                             <  2:1619027/148‑1 (MQ=255)
                                                                                   ttttCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCAAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATccc                                                             <  1:1314046/148‑1 (MQ=255)
                                                                                    tttCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATcccc                                                            >  1:500295/1‑148 (MQ=255)
                                                                                          cGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGa                                                      >  2:1467671/1‑148 (MQ=255)
                                                                                              cgcgCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCt                                                  >  1:1498646/1‑148 (MQ=255)
                                                                                                 gcCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCa                                               <  1:1408969/148‑1 (MQ=255)
                                                                                                  ccGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCAc                                              >  2:1497109/1‑148 (MQ=255)
                                                                                                    ggATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCAAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACtt                                            <  2:1083407/148‑1 (MQ=255)
                                                                                                     gATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTc                                           >  2:2061102/1‑148 (MQ=255)
                                                                                                      aTAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCAAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTcc                                          >  2:173361/1‑148 (MQ=255)
                                                                                                       tAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCg                                         <  2:2735036/148‑1 (MQ=255)
                                                                                                         aTCAAGGATGCGACCATTTTCGATAAACAGAGTCAGGTGCAAGTTATCATCAATGCCATTAACCCAGCAAATACGATCGCCTCGTCCGTTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCt                                       >  1:2998135/1‑148 (MQ=255)
                                                                                                          tCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCAAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCtt                                      >  1:3023343/1‑148 (MQ=255)
                                                                                                           cAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCAAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCttt                                     >  2:266393/1‑148 (MQ=255)
                                                                                                            aaGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCAGCGCGACGCTCCACTTCCGCTTTa                                    >  2:438349/1‑148 (MQ=255)
                                                                                                             aGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTTTaa                                   >  2:337140/1‑148 (MQ=255)
                                                                                                             aGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCAAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTTTaa                                   <  2:2529312/148‑1 (MQ=255)
                                                                                                              ggATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTTTaaa                                  >  1:1619175/1‑148 (MQ=255)
                                                                                                                aTGCGACCATTTTCGATAAACAGCGTCAGGTGCAAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACg                                >  1:563427/1‑148 (MQ=255)
                                                                                                                 tGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGt                               <  2:3388986/148‑1 (MQ=255)
                                                                                                                  gCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGtc                              <  2:2753823/148‑1 (MQ=255)
                                                                                                                   cGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGtct                             <  2:2375978/148‑1 (MQ=255)
                                                                                                                   cGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGtct                             <  2:2682566/148‑1 (MQ=255)
                                                                                                                   cGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGtct                             <  2:2682596/148‑1 (MQ=255)
                                                                                                                   cGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGtct                             <  2:2376070/148‑1 (MQ=255)
                                                                                                                   cGACAATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGtct                             <  2:1973994/148‑1 (MQ=255)
                                                                                                                      ccATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCaa                          >  1:2540055/1‑148 (MQ=255)
                                                                                                                      ccATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCaa                          >  1:2547652/1‑148 (MQ=255)
                                                                                                                       cATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAAc                         >  2:1320009/1‑148 (MQ=255)
                                                                                                                         ttttCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAAccc                       >  2:1445556/1‑148 (MQ=255)
                                                                                                                          tttCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAAcccc                      <  2:475991/148‑1 (MQ=255)
                                                                                                                              gATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCAcg                  >  2:2869762/1‑148 (MQ=255)
                                                                                                                              gATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCAcg                  >  1:1170053/1‑148 (MQ=255)
                                                                                                                                 aaaCAGCGTCAGGTGCAAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACgcgc               >  2:1726877/1‑148 (MQ=255)
                                                                                                                                 aaaCAGCGTCAGGTGCAAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACgcgc               >  2:1726865/1‑148 (MQ=255)
                                                                                                                                  aaCAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCt              >  1:2513574/1‑148 (MQ=255)
                                                                                                                                    cAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTcc            >  2:1674698/1‑148 (MQ=255)
                                                                                                                                      gCGTCAGGTGCAAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAg          >  1:985288/1‑148 (MQ=255)
                                                                                                                                         tCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCgt       >  2:2550495/1‑148 (MQ=255)
                                                                                                                                         tCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCgt       >  2:746518/1‑148 (MQ=255)
                                                                                                                                          cAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCgtg      >  1:2785026/1‑148 (MQ=255)
                                                                                                                                            ggTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTa    <  1:2551299/148‑1 (MQ=255)
                                                                                                                                            ggTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTa    <  1:2743890/148‑1 (MQ=255)
                                                                                                                                            ggTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTa    >  1:2268977/1‑148 (MQ=255)
                                                                                                                                              tGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTAtt  >  1:166302/1‑148 (MQ=255)
                                                                                                                                              tGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTAtt  <  1:859360/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
CGGTACACCGGCGATGATCAGATTCTGGTTCGCCGTAATGCGGAAATCGCCTTTGTGGATCTTCGCGATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATT  >  CP009273/2882247‑2882536

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: