Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,318,228 C→A 44.2% S251* (TCA→TAA)  yciV → PHP domain protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,318,2280CA44.2% 203.8 / 267.1 230S251* (TCA→TAA) yciVPHP domain protein
Reads supporting (aligned to +/- strand):  ref base C (59/69);  new base A (57/45);  total (116/114)
Fisher's exact test for biased strand distribution p-value = 1.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 232 total aligned reads displayed.

ACAATCTTAGTGCTAAATGGCTGAAAAGATTGGTTGCGCATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATCACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACAGAGAGGGAATTATGAGCCAGTTTTTTTATAT  >  CP009273/1318081‑1318374
                                                                                                                                                   |                                                                                                                                                  
 cAATCTTAGTGCTAAATGGCTGAAAAGATTGGTTGCGCATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATca                                                                                                                                                   >  1:2754846/1‑148 (MQ=255)
     cTTAGTGCTAAATGGCTGAAAAGATTGGTTGCGCATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAg                                                                                                                                               >  1:1077116/1‑148 (MQ=255)
      ttAGTGCTAAATGGCTGAAAAGATTGGTTGCGCATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATCACAAgg                                                                                                                                              <  2:619017/148‑1 (MQ=255)
       tAGTGCTAAATGGCTGAAAAGATTGGTTGCGCATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGa                                                                                                                                             <  2:26869/148‑1 (MQ=255)
       tAGTGCTAAATGGCTGAAAAGATTGGTTGCGCATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGa                                                                                                                                             >  2:2227328/1‑148 (MQ=255)
       tAGTGCTAAATGGCTGAAAAGATTGGTTGCGCATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGa                                                                                                                                             >  1:1344121/1‑148 (MQ=255)
          tGCTAAATGGCTGAAAAGATTGGTTGCGCATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTCCCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCt                                                                                                                                          >  2:933396/1‑148 (MQ=255)
          tGCTAAATGGCTGAAAAGATTGGTTGCGCATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATCACAAGGATCt                                                                                                                                          >  2:1549786/1‑148 (MQ=255)
              aaaTGGCTGAAAAGATTGGTTGCGCATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATCACAAGGATCTGAtt                                                                                                                                      <  1:3247362/148‑1 (MQ=255)
                 tGGCTGAAAAGATTGGTTGCGCATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTtc                                                                                                                                   >  2:3769628/1‑148 (MQ=255)
                        aaaGATTGGTTGCGCATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATCACAAGGATCTGATTTTCATCAGcc                                                                                                                            >  2:584589/1‑148 (MQ=255)
                         aaGATTGGTTGCGCATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCa                                                                                                                           >  2:900670/1‑148 (MQ=255)
                         aaGATTGGTTGCGCATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCa                                                                                                                           >  1:1776501/1‑148 (MQ=255)
                           gATTGGTTGCGCATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATCACAAGGATCTGATTTTCATCAGCCATg                                                                                                                         >  2:3462919/1‑148 (MQ=255)
                              tGGTTGCGCATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGccc                                                                                                                      >  1:622104/1‑148 (MQ=255)
                              tGGTTGCGCATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGATGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGccc                                                                                                                      >  2:815880/1‑148 (MQ=255)
                                 ttGCGCATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATCACAAGGATCTGATTTTCATCAGCCATGCCCTTg                                                                                                                   <  2:2507938/148‑1 (MQ=255)
                                    cgcATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATCACAAGGATCTGATTTTCATCAGCCATGCCCGTGGAt                                                                                                                >  1:858485/1‑148 (MQ=255)
                                       aTTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATCACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGa                                                                                                             <  1:3169786/148‑1 (MQ=255)
                                         tttGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAAc                                                                                                           >  2:2118535/1‑148 (MQ=255)
                                         tttGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAAc                                                                                                           >  1:2304537/1‑148 (MQ=255)
                                               gAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGTCCGTGGATCGAACTGGGTc                                                                                                     >  1:3531497/1‑148 (MQ=255)
                                                 acacCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATCACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGt                                                                                                   >  1:896817/1‑148 (MQ=255)
                                                 acacCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATCACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGt                                                                                                   >  1:254585/1‑148 (MQ=255)
                                                    ccacGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATCACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTaaa                                                                                                <  1:2439771/148‑1 (MQ=255)
                                                    ccacGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATCACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTaaa                                                                                                <  1:2049785/148‑1 (MQ=255)
                                                     cacGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAAc                                                                                               >  1:257241/1‑148 (MQ=255)
                                                      acGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAAct                                                                                              >  2:492754/1‑148 (MQ=255)
                                                      acGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAAct                                                                                              >  2:3389843/1‑148 (MQ=255)
                                                      acGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAAct                                                                                              <  2:1235339/148‑1 (MQ=255)
                                                       cGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATCACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAActc                                                                                             <  2:2754846/148‑1 (MQ=255)
                                                             cgcgATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATCACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTg                                                                                       >  1:1873614/1‑148 (MQ=255)
                                                             cgcgATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATCACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTg                                                                                       >  1:2998987/1‑148 (MQ=255)
                                                               cgATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGcc                                                                                     <  2:3310339/148‑1 (MQ=255)
                                                                  tGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATCACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGc                                                                                  >  2:3381529/1‑148 (MQ=255)
                                                                  tGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATCACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGc                                                                                  >  2:683740/1‑148 (MQ=255)
                                                                     aaGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATCACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAgg                                                                               >  1:335424/1‑148 (MQ=255)
                                                                       gTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCg                                                                             >  2:1108255/1‑148 (MQ=255)
                                                                       gTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCg                                                                             <  1:3769628/148‑1 (MQ=255)
                                                                       gTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGATGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCg                                                                             <  1:426618/148‑1 (MQ=255)
                                                                       gTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGATGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCg                                                                             <  1:815880/148‑1 (MQ=255)
                                                                           cgcAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGa                                                                         >  2:1342169/1‑148 (MQ=255)
                                                                            gcAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATCACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGaa                                                                        <  2:3049677/148‑1 (MQ=255)
                                                                             cAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATCACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAg                                                                       <  2:3290922/148‑1 (MQ=255)
                                                                              aGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATCACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAgg                                                                      <  2:2854926/148‑1 (MQ=255)
                                                                                 gCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATCACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGt                                                                   <  2:2816898/148‑1 (MQ=255)
                                                                                 gCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATCACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGt                                                                   <  2:3358157/148‑1 (MQ=255)
                                                                                 gCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGt                                                                   <  1:3633968/148‑1 (MQ=255)
                                                                                 gCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGt                                                                   <  2:1130194/148‑1 (MQ=255)
                                                                                 gCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGt                                                                   <  2:824457/148‑1 (MQ=255)
                                                                                   cAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATCACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCt                                                                 <  1:1178832/148‑1 (MQ=255)
                                                                                    aGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCAGTTATGGGCATCACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTg                                                                >  2:2162761/1‑148 (MQ=255)
                                                                                      cAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATCACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGc                                                              <  2:2614140/148‑1 (MQ=255)
                                                                                      cAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATCACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGc                                                              <  1:1299159/148‑1 (MQ=255)
                                                                                         tCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATCACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGc                                                           >  2:838218/1‑148 (MQ=255)
                                                                                           gCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATCACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTa                                                         <  2:2232557/148‑1 (MQ=255)
                                                                                                aTGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATCACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGa                                                    >  2:632321/1‑148 (MQ=255)
                                                                                                aTGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATCACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGa                                                    >  2:1285446/1‑148 (MQ=255)
                                                                                                aTGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATCACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGa                                                    >  2:1566821/1‑148 (MQ=255)
                                                                                                 tGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATCACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGaa                                                   >  1:1001561/1‑148 (MQ=255)
                                                                                                  gAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAAc                                                  <  1:1367371/148‑1 (MQ=255)
                                                                                                     cGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATCACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGc                                               >  2:2431389/1‑148 (MQ=255)
                                                                                                     cGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGc                                               >  2:2004893/1‑148 (MQ=255)
                                                                                                       tACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCg                                             <  1:2225036/148‑1 (MQ=255)
                                                                                                               tGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATCACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAacac                                     <  2:1654533/148‑1 (MQ=255)
                                                                                                                 gcggcCCTTGCGCGTCAGCATCATTTATGGGCATCACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAAcacca                                   <  2:254585/148‑1 (MQ=255)
                                                                                                                 gcggcCCTTGCGCGTCAGCATCATTTATGGGCATCACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAAcacca                                   <  2:896797/148‑1 (MQ=255)
                                                                                                                  cggcCCTTGCGCGTCAGCATCATTTATGGGCATCACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAAcaccac                                  <  1:2641465/148‑1 (MQ=255)
                                                                                                                  cggcCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAAcaccac                                  <  2:1168021/148‑1 (MQ=255)
                                                                                                                  cggcCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAAcaccac                                  >  1:2649321/1‑148 (MQ=255)
                                                                                                                      ccTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACagag                              >  1:451468/1‑148 (MQ=255)
                                                                                                                      ccTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACagag                              <  1:3134208/148‑1 (MQ=255)
                                                                                                                       cTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACagaga                             >  1:2833661/1‑148 (MQ=255)
                                                                                                                       cTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACagaga                             >  1:487507/1‑148 (MQ=255)
                                                                                                                       cTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCAGCAGAACACCACagaga                             <  1:2343652/148‑1 (MQ=255)
                                                                                                                        ttGCGCGTCAGCATCATTTATGGGCATCACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACagagag                            >  2:1319620/1‑148 (MQ=255)
                                                                                                                          gcgcGTCAGCATCATTTATGGGCATCACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACAGAGAggg                          <  2:2189795/148‑1 (MQ=255)
                                                                                                                            gcgTCAGCATCATTTATGGGCATCACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACAGAGAGGGaa                        <  2:3276808/148‑1 (MQ=255)
                                                                                                                            gcgTCAGCATCATTTATGGGCATCACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACAGAGAGGGaa                        <  1:3398846/148‑1 (MQ=255)
                                                                                                                              gTCAGCATCATTTATGGGCATCACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACAGAGAGGGAAtt                      <  1:921109/148‑1 (MQ=255)
                                                                                                                               tCAGCAGCATTTATGGGCATCACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACAGAGAGGGAATTa                     >  2:3567585/1‑148 (MQ=255)
                                                                                                                                cAGCATCATTTATGGGCATCACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACAGAGAGGGAATTAt                    >  1:1125523/1‑148 (MQ=255)
                                                                                                                                cAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACAGAGAGGGAATTAt                    >  2:2455677/1‑148 (MQ=255)
                                                                                                                                  gCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACAGAGAGGGAATTATGa                  <  2:1776501/148‑1 (MQ=255)
                                                                                                                                   catcatTTATGGGCATCACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACAGAGAGGGAATTATGAg                 >  1:1237936/1‑148 (MQ=255)
                                                                                                                                    atcatTTATTGGCATCACAAGGATCTGATTTTGATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGTAACAGCCGCAGAACACCACAGAGAGGGAATTATGAGc                <  2:2945932/148‑1 (MQ=255)
                                                                                                                                    atcatTTATGGGCATCACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACAGAGAGGGAATTATGAGc                >  1:181309/1‑148 (MQ=255)
                                                                                                                                       atTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACAGAGAGGGAATTATGAGCCAg             <  1:1915835/148‑1 (MQ=255)
                                                                                                                                        tttATGGGCATCACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACAGAGAGGGAATTATGAGCCAGt            <  2:688621/148‑1 (MQ=255)
                                                                                                                                        tttATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACAGAGAGGGAATTATGAGCCAGt            <  1:1992512/148‑1 (MQ=255)
                                                                                                                                         ttATGGGCATCACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACAGAGAGGGAATTATGAGCCAGtt           >  2:2819967/1‑148 (MQ=255)
                                                                                                                                          tATGGGCATCACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACAGAGAGGGAATTATGAGCCAGttt          <  1:3437181/148‑1 (MQ=255)
                                                                                                                                          tATGGGCATCACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACAGAGAGGGAATTATGAGCCAGttt          <  1:1643342/148‑1 (MQ=255)
                                                                                                                                             gggCATCACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACAGAGAGGGAATTATGAGCCAGtttttt       <  1:637067/148‑1 (MQ=255)
                                                                                                                                              ggCATCACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACAGAGAGGGAATTATGAGCCAGttttttt      >  1:2924829/1‑148 (MQ=255)
                                                                                                                                                cATCACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACAGAGAGGGAATTATGAGCCAGTTTTTTtat    <  1:298680/148‑1 (MQ=255)
                                                                                                                                                cATCACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACAGAGAGGGAATTATGAGCCAGTTTTTTtat    <  2:2998987/148‑1 (MQ=255)
                                                                                                                                                cATCACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACAGAGAGGGAATTATGAGCCAGTTTTTTtat    <  2:1874617/148‑1 (MQ=255)
                                                                                                                                                cATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACAGAGAGGGAATTATGAGCCAGTTTTTTtat    <  2:2893392/148‑1 (MQ=255)
                                                                                                                                                  tCACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACAGAGAGGGAATTATGAGCCAGTTTTTTTatat  <  1:1549786/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
ACAATCTTAGTGCTAAATGGCTGAAAAGATTGGTTGCGCATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATCACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACAGAGAGGGAATTATGAGCCAGTTTTTTTATAT  >  CP009273/1318081‑1318374

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: