Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,223,839 G→A 18.2% D52N (GAT→AAT)  pgi → glucosephosphate isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,223,8390GA18.2% 504.3 / 88.0 192D52N (GAT→AAT) pgiglucosephosphate isomerase
Reads supporting (aligned to +/- strand):  ref base G (84/73);  new base A (13/22);  total (97/95)
Fisher's exact test for biased strand distribution p-value = 9.39e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

Only 100 of 201 total aligned reads displayed.

AACATCAATCCAACGCAGACCGCTGCCTGGCAGGCACTACAGAAACACTTCGATGAAATGAAAGACGTTACGATCGCCGATCTTTTTGCTAAAGACGGCGATCGTTTTTCTAAGTTCTCCGCAACCTTCGACGATCAGATGCTGGTGGATTACTCCAAAAACCGCATCACTGAAGAGACGCTGGCGAAATTACAGGATCTGGCGAAAGAGTGCGATCTGGCGGGCGCGATTAAGTCGATGTTCTCTGGCGAGAAGATCAACCGCACTGAAAACCGCGCCGTGCTGCACGTAGCGC  >  CP009273/4223692‑4223986
                                                                                                                                                   |                                                                                                                                                   
aaCATCAATCCAACGCAGACCGCTGCCTGGCAGGCACTACAGAAACACTTCGATGAAATGAAAGACGTTACGATCGCCGATCTTTTTGCTAAAGACGGCGATCGTTTTTCTAAGTTCTCCGCAACCTTCGACGATCAGATGCtggtgg                                                                                                                                                     <  1:3151900/148‑1 (MQ=255)
 aCATCAATCCAACGCAGACCGCTGCCTGGCAGGCACTACAGAAACACTTCGATGAAATGAAAGACGTTACGATCGCCGATCTTTTTGCTAAAGACGGCGATCGTTTTTCTAAGTTCTCCGCAACCTTCGACGATCAGATGCTGGTGGa                                                                                                                                                    <  2:3262190/148‑1 (MQ=255)
        tCCAACGCAGACCGCTGCCTGGCAGGCACTACAGAAACACTTCGATGAAATGAAAGACGTTACGATCGCCGATCTTTTTGCTAAAGACGGCGATCGTTTTTCTAAGTTCTCCGCAACCTTCGACGATCAGATGCTGGTGGATTACTcc                                                                                                                                             >  1:2260763/1‑148 (MQ=255)
         ccAACGCAGACCGCTGCCTGGCAGGCACTACAGAAACACTTCGATGAAATGAAAGACGTTACGATCGCCGATCTTTTTGCTAAAGACGGCGATCGTTTTTCTAAGTTCTCCGCAACCTTCGACGATCAGATGCTGGTGGATTACTCCa                                                                                                                                            >  2:3377086/1‑148 (MQ=255)
         ccAACGCAGACCGCTGCCTGGCAGGCACTACAGAAACACTTCGATGAAATGAAAGACGTTACGATCGCCGATCTTTTTGCTAAAGACGGCGATCGTTTTTCTAAGTTCTCCGCAACCTTCGACGATCAGATGCTGGTGACTTACTCCa                                                                                                                                            >  1:2858968/1‑148 (MQ=255)
         ccAACGCAGACCGCTGCCTGGCAGGCACTACAGAAACACTTCGATGAAATGAAAGACGTTACGATCGCCGATCTTTTTGCTAAAGACGGCGATCGTTTTTCTAAGTTCTCCGCAACCTTCGACGATCAGATGCGGGTGGATTACTCCa                                                                                                                                            >  1:1258947/1‑148 (MQ=255)
                 gACCGCTGCCTGGCAGGCACTACAGAAACACTTCGATGAAATGAAAGACGTTACGATCGCCGATCTTTTTGCTAAAGACGGCGATCGTTTTTCTAAGTTCTCCGCAACCTTCGACGATCAGATGCTGGTGGATTACTCCAAAAACCGc                                                                                                                                    >  1:93725/1‑148 (MQ=255)
                  aCCGCTGCCTGGCAGGCACTACAGAAACACTTCGATGAAATGAAAGACGTTACGATCGCCGATCTTTTTGCTAAAGACGGCGATCGTTTTTCTAAGTTCTCCGCAACCTTCGACGATCAGATGCTGGTGGATTACTCCAAAAACCGCa                                                                                                                                   >  1:3222959/1‑148 (MQ=255)
                  aCCGCTGCCTGGCAGGCACTACAGAAACACTTCGATGAAATGAAAGACGTTACGATCGCCGATCTTTTTGCTAAAGACGGCGATCGTTTTTCTAAGTTCTCCGCAACCTTCGACGATCAGATGCTGGTGGATTACTCCAAAAACCGCa                                                                                                                                   <  1:1285096/148‑1 (MQ=255)
                  aCCGCTGCCTGGCAGGCACTACAGAAACACTTCGATGAAATGAAAGACGTTACGATCGCCGATCTTTTTGCTAAAGACGGCGATCGTTTTTCTAAGTTCTCCGCAACCTTCGACGATCAGATGCTGGTGGATTACTCCAAAAACCGCa                                                                                                                                   >  1:3413171/1‑148 (MQ=255)
                  aCCGCTGCCTGGCAGGCACTACAGAAACACTTCGATGAAATGAAAGACGTTACGATCGCCGATCTTTTTGCTAAAGACGGCGATCGTTTTTCTAAGTTCTCCGCAACCTTCGACGATCAGATGCTGGTGGATTACTCCAAAAACCGCa                                                                                                                                   >  1:2935252/1‑148 (MQ=255)
                         ccTGGCAGGCACTACAGAAACACTTCGATGAAATGAAAGACGTTACGATCGCCGATCTTTTTGCTAAAGACGGCGATCGTTTTTCTAAGTTCTCCGCAACCTTCGACGATCAGATGCTGGTGGATTACTCCAAAAACCGCATCACTGa                                                                                                                            >  2:1292682/1‑148 (MQ=255)
                          cTGGCAGGCACTACAGAAACACTTCGATGAAATGAAAGACGTTACGATCGCCGATCTTTTTGCTAAAGACGGCGATCGTTTTTCTAAGTTCTCCGCAACCTTCGACGATCAGATGCTGGTGGATTACTCCAAAAACCGCATCACTGaa                                                                                                                           <  1:3029107/148‑1 (MQ=255)
                           tGGCAGGCACTACAGAAACACTTCGATGAAATGAAAGACGTTACGATCGCCGATCTTTTTGCTAAAGACGGCGATCGTTTTTCTAAGTTCTCCGCAACCTTCGACGATCAGATGCTGGTGGATTACTCCAAAAACCGCATCACTGAag                                                                                                                          <  2:3382050/148‑1 (MQ=255)
                            ggcaggcaCTACAGAAACACTTCGATGAAATGAAAGACGTTACGATCGCCGATCTTTTTGCTAAAGACGGCGATCGTTTTTCTAAGTTCTCCGCAACCTTCGACGATCAGATGCTGGTGAATTACTCCAAAAACCGCATCACTGAaga                                                                                                                         <  2:2548345/148‑1 (MQ=255)
                             gcaggcaCTACAGAAACACTTCGATGAAATGAAAGACGTTACGATCGCCGATCTTTTTGCTAAAGACGGCGATCGTTTTTCTAAGTTCTCCGCAACCTTCGACGATCAGATGCTGGTGGATTACTCCAAAAACCGCATCACTGAagag                                                                                                                        >  1:219749/1‑148 (MQ=255)
                              caggcaCTACAGAAACACTTCGATGAAATGAAAGACGTTACGATCGCCGATCTTTTTGCTAAAGACGGCGATCGTTTTTCTAAGTTCTCCGCAACCTTCGACGATCAGATGCTGGTGGATTACTCCAAAAACCGCATCACTGAAgaga                                                                                                                       >  1:2546988/1‑148 (MQ=255)
                              caggcaCTACAGAAACACTTCGATGAAATGAAAGACGTTACGATCGCCGATCTTTTTGCTAAAGACGGCGATCGTTTTTCTAAGTTCTCCGCAACCTTCGACGATCAGATGCTGGTGGATTACTCCAAAAACCGCATCACTGAAgaga                                                                                                                       >  2:1921381/1‑148 (MQ=255)
                                 gcaCTACAGAAACACTTCGATGAAATGAAAGACGTTACGATCGCCGATCTTTTTGCTAAAGACGGCGATCGTTTTTCTAAGTTCTCCGCAACCTTCGACGATCAGATGCTGGTGGATTACTCCAAAAACCGCATCACTGAAGAGACGc                                                                                                                    >  2:1346191/1‑148 (MQ=255)
                                  caCTACATAAACACTTCGATGAAATGAAAGACGTTACGATCGCCGATCTTTTTGCTAAAGACGGCGATCGTTTTTCTAAGTTCTCCGCAACCTTCGACGATCAGATGCTGGTGAATTACTCCAAAAACCGCATCACTGAAGAGACGCt                                                                                                                   <  1:40879/148‑1 (MQ=255)
                                  caCTACAGAAACACTTCGATGAAATGAAAGACGTTACGATCGCCGATCTTTTTGCTAAAGACGGCGATCGTTTTTCTAAGTTCTCCGCAACCTTCGACGATCAGATGCTGGTGAATTACTCCAAAAACCGCATCACTGAAGAGACGCt                                                                                                                   <  1:1505383/148‑1 (MQ=255)
                                      aCAGAAACACTTCGATGAAATGAAAGACGTTACGATCGCCGATCTTTTTGCTAAAGACGGCGATCGTTTTTCTAAGTTCTCCGCAACCTTCGACGATCAGATGCTGGTGGATTACTCCAAAAACCGCATCACTGAAGAGACGCTGGCg                                                                                                               >  2:2642879/1‑148 (MQ=255)
                                        aGAAACACTTCGATGAAATGAAAGACGTTACGATCGCCGATCTTTTTGCTAAAGACGGCGATCGTTTTTCTAAGTTCTCCGCAACCTTCGACGATCAGATGCTGGTGGATTACTCCAAAAACCGCATCACTGAAGAGACGCTGGCGaa                                                                                                             <  2:1525743/148‑1 (MQ=255)
                                        aGAAACACTTCGATGAAATGAAAGACGTTACGATCGCCGATCTTTTTGCTAAAGACGGCGATCGTTTTTCTAAGTTCTCCGCAACCTTCGACGATCAGATGCTGGTGGATTACTCCAAAAACCGCATCACTGAAGAGACGCTGGCGaa                                                                                                             <  1:1480264/148‑1 (MQ=255)
                                         gAAACACTTCGATGAAATGAAAGACGTTACGATCGCCGATCTTTTTGCTAAAGACGGCGATCGTTTTTCTAAGTTCTCCGCAACCTTCGACGATCAGATGCTGGTGGATTACTCCAAAAACCGCATCACTGAAGAGACGCTGGCGaaa                                                                                                            <  2:2865318/148‑1 (MQ=255)
                                          aaaCACTTCGATGAAATGAAAGACGTTACGATCGCCGATCTTTTTGCTAAAGACGGCGATCGTTTTTCTAAGTTCTCCGCAACCTTCGACGATCAGATGCTGGTGGATTACTCCAAAAACCGCATCACTGAAGAGACGCTGGCGAAAt                                                                                                           >  1:595233/1‑148 (MQ=255)
                                            acacTTCGATGAAATGAAAGACGTTACGATCGCCGATCTTTTTGCTAAAGACGGCGATCGTTTTTCTAAGTTCTCCGCAACCTTCGACGATCAGATGCTGGTGGATTACTCCAAAAACCGCATCACTGAAGAGACGCTGGCGAAATTa                                                                                                         <  2:3299878/148‑1 (MQ=255)
                                            acacTTCGATGAAATGAAAGACGTTACGATCGCCGATCTTTTTGCTAAAGACGGCGATCGTTTTTCTAAGTTCTCCGCAACCTTCGACGATCAGATGCTGGTGGATTACTCCAAAAACCGCATCACTGAAGAGACGCTGGCGAAATTa                                                                                                         <  1:1025695/148‑1 (MQ=255)
                                             cacTTCGATGAAATGAAAGACGTTACGATCGCCGATCTTTTTGCTAAAGACGGCGATCGTTTTTCTAAGTTCTCCGCAACCTTCGACGATCAGATGCTGGTGAATTACTCCAAAAACCGCATCACTGAAGAGACGCTGGCGAAATTAc                                                                                                        <  1:494000/148‑1 (MQ=255)
                                               cTTCGATGAAATGAAAGACGTTACGATCGCCGATCTTTTTGCTAAAGACGGCGATCGTTTTTCTAAGTTCTCCGCAACCTTCGACGATCAGATGCTGGTGGATTACTCCAAAAACCGCATCACTGAAGAGACGCTGGCGAAATTACAg                                                                                                      >  1:1340066/1‑148 (MQ=255)
                                                ttCGATGAAATGAAAGACGTTACGATCGCCGATCTTTTTGCTAAAGACGGCGATCGTTTTTCTAAGTTCTCCGCAACCTTCGACGATCAGATGCTGGTGGATTACTCCAAAAACCGCATCACTGAAGAGACGCTGGCGAAATTACAgg                                                                                                     >  1:468588/1‑148 (MQ=255)
                                                 tCGATGAAATGAAAGACGTTACGATCGCCGATCTTTTTGCTAAAGACGGCGATCGTTTTTCTAAGTTCTCCGCAACCTTCGACGATCAGATGCTGGTGGATTACTCCAAAAACCGCATCACTGAAGAGACGCTGGCGAAATTACAGGa                                                                                                    <  1:2615633/148‑1 (MQ=255)
                                                    atgaaatgaaAGACGTTACGATCGCCGATCTTTTTGCTAAAGACGGCGATCGTTTTTCTAAGTTCTCCGCAACCTTCGACGATCAGATGCTGGTGGATTACTCCAAAAACCGCATCACTGAAGAGACGCTGGCGAAATTACAGGATCt                                                                                                 <  1:1548932/148‑1 (MQ=255)
                                                     tgaaatgaaaGACGTTACGATCGCCGATCTTTTTGCTAAAGACGGCGATCGTTTTTCTAAGTTCTCCGCAACCTTCGACGATCAGATGCTGGTGGATTACTCCAAAAACCGCATCACTGAAGAGACGCTGGCGAAATTACAGGATCTg                                                                                                >  1:3457442/1‑148 (MQ=255)
                                                     tgaaatgaaaGACGTTACGATCGCCGATCTTTTTGCTAAAGACGGCGATCGTTTTTCTAAGTTCTCCGCAACCATCGACGATCAGATGCTGGTGGATTACTCCAAAAACCGCATCACTGAAGAGACGCTGGCGAAATTACAGGATCTg                                                                                                >  1:1228705/1‑148 (MQ=255)
                                                       aaatgaaaGACGTTACGATCGCCGATCTTTTTGCTAAAGACGGCGATCGTTTTTCTAAGTTCTCCGCAACCTTCGACGATCAGATGCTGGTGGATTACTCCAAAAACCGCATCACTGAAGAGACGCTGGCGAAATTACAGGATCTGGc                                                                                              >  2:1936421/1‑148 (MQ=255)
                                                       aaatgaaaGACGTTACGATCGCCGATCTTTTTGCTAAAGACGGCGATCGTTTTTCTAAGTTCTCCGCAACCTTCGACGATCAGATGCTGGTGGATTACTCCAAAAACCGCATCACTGAAGAGACGCTGGCGAAATTACAGGATCTGGc                                                                                              >  2:679096/1‑148 (MQ=255)
                                                            aaaGACGTTACGATCGCCGATCTTTTTGCTAAATACGGCGATCGTTTTTCTAATTTCTCCGCAACCTTCGACGATCAGATGCTGGTGGATTACTCCAAAAAACGCATCACTGAAGAGACGCTGGCGAAATTACAGGATCTGGCGAAag                                                                                         >  1:466807/1‑148 (MQ=255)
                                                             aaGACGTTACGATCGCCGATCTTTTTGCTAAAGACGGCGATCGTTTTTCTAAGTTCTCCGCAACCTTCGACGATCAGATGCTGGTGGATTACTCCAAAAACCGCATCACTGAAGAGACGCTGGCGAAATTACAGGATCTGGCGAAaga                                                                                        <  1:2713918/148‑1 (MQ=255)
                                                             aaGACGTTACGATCGCCGATCTTTTTGCTAAAGACGGCGATCGTTTTTCTAAGTTCTCCGCAACCTTCGACGATCAGATGCTGGTGGATTACTCCAAAAACCGCATCACTGAAGAGACGCTGGCGAAATTACAGGATCTGGCGAAaga                                                                                        <  1:1402381/148‑1 (MQ=255)
                                                               gACGTTACGATCGCCGATCTTTTTGCTAAAGACGGCGATCGTTTTTCTAAGTTCTCCGCAACCTTCGACGATCAGATGCTGGTGGATTACTCCAAAAACCGCATCACTGAAGAGACGCTGGCGAAATTACAGGATCTGGCGAAAGAGt                                                                                      <  1:790127/148‑1 (MQ=255)
                                                               gACGTTACGATCGCCGATCTTTTTGCTAAAGACGGCGATCGTTTTTCTAAGTTATCCGCAACCATCGACGATCAGATGCTGGTGGATTACTACAAAAAACGCATCACTGAAGAGACGCTGGCGAAATTAAAGGATATGGCGAAAGAGt                                                                                      >  1:1681109/1‑148 (MQ=255)
                                                                 cGTTACGATCGCCGATCTTTTTGCTAAAGACGGCGATCGTTTTTCTAAGTTCTCCGCAACCTTCGACGATCAGATGCTGGTGGATTACTCCAAAAACCGCATCACTGAAGAGACGCTGGCGAAATTACAGGATCTGGCGAAAGAGTGc                                                                                    >  1:3096275/1‑148 (MQ=255)
                                                                  gTTACGATCGCCGATCTTTTTGCTAAAGACGGCGATCGTTTTTCTAAGTTCTCCGCAACCTTCGACGATCAGATGCTGGTGGATTACTCCAAAAACCGCATCACTGAAGAGACGCTGGCGAAATTACAGGATCTGGCGAAAGAGTGCg                                                                                   <  2:2743086/148‑1 (MQ=255)
                                                                   ttACGATCGCCGATCTTTTTGCTAAAGACGGCGATCGTTTTTCTAAGTTCTCCGCAACCTTCGACGATCAGATGCTGGTGGATTACTCCAAAAACCGCATCACTGAAGAGACGCTGGCGAAATTACAGGATCTGGCGAAAGAGTGCGa                                                                                  >  1:3331522/1‑148 (MQ=255)
                                                                   ttACGATCGCCGATCTTTTTGCTAAAGACGGCGATCGTTTTTCTAAGTTCTCCGCAACCTTCGACGATCAGATGCTGGTGGATTACTCCAAAAACCGCATCACTGAAGAGACGCTGGCGAAATTACAGGATCTGGCGAAAGAGTGCGa                                                                                  <  2:1567458/148‑1 (MQ=255)
                                                                   ttACGATCGCCGATCTTTTTGCTAAAGACGGCGATCGTTTTTCTAAGTTCTCCGCAACCTTCGACGATCAGATGCTGGTGGATTACTCCAAAAACCGCATCACTGAAGAGACGCTGGCGAAATTACAGGATCTGGCGAAAGAGTGCGa                                                                                  <  2:3336408/148‑1 (MQ=255)
                                                                   ttACGATCGCCGATCTTTTTGCTAAAGACGGCGATCGTTTTTCTAAGTTCTCCGCAACCTTCGACGATCAGATGCTGGTGAATTACTCCAAAAACCGCATCACTGAAGAGACGCTGGCGAAATTACAGGATCTGGCGAAAGAGTGCGa                                                                                  <  2:1811816/148‑1 (MQ=255)
                                                                             cGATCTTTTTGCTAAAGACGGCGATCGTTTTTCTAAGTTCTCCGCAACCTTCGACGATCAGATGCTGGTGAATTACTCCAAAAACCGCATCACTGAAGAGACGCTGGCGAAATTACAGGATCTGGCGAAAGAGTGCGATCTggcgggc                                                                        <  1:2308031/148‑1 (MQ=255)
                                                                               aTCTTTTTGCTAAAGACGGCGATCGTTTTTCTAAGTTCTCCGCAACCTTCGACGATCAGATGCTGGTGAATTACTCCAAAAACCGCATCACTGAAGAGACGCTGGCGAAATTACAGGATCTGGCGAAAGAGTGCGATCTGGCGGgcgc                                                                      >  1:1851727/1‑148 (MQ=255)
                                                                                tCTTTTTGCTAAAGACGGCGATCGTTTTTCTAAGTTCTCCGCAACCTTCGACGATCAGATGCTGGTGGATTACTCCAAAAACCGCATCACTGAAGAGACGCTGGCGAAATTACAGGATCTGGCGAAAGAGTGCGATCTGGCGGGcgcg                                                                     >  1:1012227/1‑148 (MQ=255)
                                                                                tCTTTTTGCTAAAGACGGCGATCGTTTTTCTAAGTTCTCCGCAACCTTCGACGATCAGATGCTGGTGGATTACTCCAAAAACCGCATCACTGAAGAGACGCTGGCGAAATTACAGGATCTGGCGAAAGAGTGCGATCTGGCGGGcgcg                                                                     <  2:1706188/148‑1 (MQ=255)
                                                                                 cTTTTTGCTAAAGACGGCGATCGTTTTTCTAAGTTCTCCGCAACCTTCGACGATCAGATGCTGGTGGATTACTCCAAAAACCGCATCACTGAAGAGACGCTGGCGAAATTACAGGATCTGGCGAAAGAGTGCGATCTGGCGGGCGCGa                                                                    >  2:1187233/1‑148 (MQ=255)
                                                                                        cTAAAGACGGCGATCGTTTTTCTAAGTTCTCCGCAACCTTCGACGATCAGATGCTGGTGAATTACTCCAAAAACCGCATCACTGAAGAGACGCTGGCGAAATTACAGGATCTGGCGAAAGAGTGCGATCTGGCGGGCGCGATTAAGTc                                                             >  1:2659327/1‑148 (MQ=255)
                                                                                         tAAAGACGGCGATCGTTTTTCTAAGTTCTCCGCAACCTTCGACGATCAGATGCTGGTGGATTACTCCAAAAACCGCATCACTGAAGAGACGCTGGCGAAATTACAGGATCTGGCGAAAGAGTGCGATCTGGCGGGCGCGATTAAGTCg                                                            <  1:759563/148‑1 (MQ=255)
                                                                                          aaaGACGGCGATCGTTTTTCTAAGTTCTCCGCAACCTTCGACGATCAGATGCTGGTGGATTACTCCAAAAACCGCATCACTGAAGAGACGCTGGCGAACTTACAGGATCTGGCGAAAGAGTGCGATCTGGCGGGCGCGATTAAGTCGa                                                           >  2:282866/1‑148 (MQ=255)
                                                                                          aaaGACGGCGATCGTTTTTCTAAGTTCTCCGCAACCTTCGACGATCAGATGCTGGTGAATTACTCCAAAAACCGCATCACTGAAGAGACGCTGGCGAAATTACAGGATCTGGCGAAAGAGTGCGATCTGGCGGGCGCGATTAAATCGa                                                           >  2:300128/1‑148 (MQ=255)
                                                                                             gACGGCGATCGTTTTTCTAAGTTCTCCGCAACCTTCGACGATAAGATGATGGTGGATTACTCCAAAAACCGCATCACTGAAGAGACGCTGGCGAAATTACAGGATCTGGCGAAAGAGTTCGATCTGTCGGGCGCGATTAAGTCGATGt                                                        <  1:507410/148‑1 (MQ=255)
                                                                                             gACGGCGATCGTTTTTCTAAGTTCTCCGCAACCTTCGACAATCAGATGCTGGTGGATTACTCCAAAAACCGCATCACTGAAGAGACGCTGGCGAAATTACAGGATCTGGCGAAAGAGTGCGATCTGGCGGGCGCGATTAAGTCGATGt                                                        >  2:470287/1‑148 (MQ=255)
                                                                                              aCGGCGATCGTTTTTCTAAGTTCTCCGCAACCTTCGACGATCAGATGCTGGTGGATTACTCCAAAAACCGCATCACTGAAGAGACGCTGGCGAAATTACAGGATCTGGCGAAAGAGTGCGATCTGGCGGGCGCGATTAAGTCGATGtt                                                       >  1:1975636/1‑148 (MQ=255)
                                                                                                ggCGATCGTTTTTCTAAGTTCTCCGCAACCTTCGACGATCAGATGCTGGTGGATTACTCCAAAAACCGCATCACTGAAGAGACGCTGGCGAAATTACAGGATCTGGCGAAAGAGTGCGATCTGGCGGGCGCGATTAAGTCGATGTtct                                                     <  2:3066479/148‑1 (MQ=255)
                                                                                                      cGTTTTTCTAAGTTCTCCGCAACCTTCGACGATCAGATGCTGGTGGATTACTCCAAAAACCGCATCACTGAAGAGACGCTGGCGAAATTACAGGATCTGGCGAAAGAGTGCGATCTGGCGGGCGCGATTAAGTCGATGTTCTCTGGCg                                               >  1:817957/1‑148 (MQ=255)
                                                                                                      cGTTTTTCTAAGTTCTCCGCAACCTTCGACGATCAGATGCTGGTGGATTACTCCAAAAACCGCATCACTGAAGAGACGCTGGCGAAATTACAGGATCTGGCGAAAGAGTGCGATCTGGCGGGCGCGATTAAGTCGATGTTCTCTGGCg                                               >  2:65896/1‑148 (MQ=255)
                                                                                                        tttttCTAAGTTCTCCGCAACCTTCGACGATCAGATGCTGGTGGATTACTCCAAAAACCGCATCACTGAAGAGACGCTGGCGAAATTACAGGATCTGGCGAAAGAGTGCGATCTGGCGGGCGCGATTAAGTCGATGTTCTCTGGCgag                                             <  2:2765566/148‑1 (MQ=255)
                                                                                                        tttttCTAAGTTCTCCGCAACCTTCGACGATCAGATGCTGGTGGATTACTCCAAAAACCGCATCACTGAAGAGACGCTGGCGAAATTACAGGATCTGGCGAAAGAGTGCGATCTGGCGGGCGCGATTAAGTCGATGTTCTCTGGCgag                                             >  1:924708/1‑148 (MQ=255)
                                                                                                        tttttCTAAGTTCTCCGCAACCTTCGACGATCAGATGCTGGTGGATTACTCCAAAAACCGCATCACTGAAGAGACGCTGGCGAAATTACAGGATCTGGCGAAAGAGTGCGATCTGGCGGGCGCGATTAAGTCGATGTTCTCTGGCgag                                             <  2:2907140/148‑1 (MQ=255)
                                                                                                        tttttCTAAGTTCTCCGCAACCTTCGACGATCAGATGCTGGTGAATTACTCCAAAAACCGCATCACTGAAGAGACGCTGGCGAAATTACAGGATCTGGCGAAAGAGTGCGATCTGGCGGGCGCGATTAAGTCGATGTTCTCTGGCgag                                             >  2:2351996/1‑148 (MQ=255)
                                                                                                         ttttCTAAGTTCTCCGCAACCTTCGACGATCAGATGCTGGTGGATTACTCCAAAAACCGCATCACTGAAGAGACGCTGGCGAAATTACAGGATCTGGCGAAAGAGTGCGATCTGGCGGGCGCGATTAAGTCGATGTTCTCTGGCgaga                                            >  2:2443861/1‑148 (MQ=255)
                                                                                                            tCTAAGTTCTCCGCAACCTTCGACGATCAGATGCTGGTGGATTACTCCAAAAACCGCATCACTGAAGAGACGCTGGCGAAATTACAGGATCTGGCGAAAGAGTGCGATCTGGCGGGCGCGATTAAGTCGATGTTCTCTGGCGagaaga                                         <  1:657527/148‑1 (MQ=255)
                                                                                                            tCTAAGTTCTCCGCAACCTTCGACGATCAGATGCTGGTGGATTACTCCAAAAACCGCATCACTGAAGAGACGCTGGCGAAATTACAGGATCTGGCGAAAGAGTGCGATCTGGCGGGCGCGATTAAGTCGATGTTCTCTGGCGagaaga                                         <  2:194619/148‑1 (MQ=255)
                                                                                                                aggtctcCGCAACCTTCGACGATCAGATGCTGGTGGATTACTCCAAAAACCGCATCACTGAAGAGACGCTGGCGAAATTACAGGATCTGGCGAAAGAGTGCGATCTGGCGGGCGCGATTAAGTCGATGTTCTCTGGCGAGAAGATCaa                                     <  2:2690148/145‑1 (MQ=255)
                                                                                                                aGTTCTCCGCAACCTTCGACGCTCAGATGCTGGTGAATTACTCCAAAAACCGCATCACTGAAGAGACGCTGGCGAAATTACAGGATCTGGCGAAAGAGTGCGATCTGGCGGGCGCGATTAAGTCGATGTTCTCTGGCGAGAAGATCaa                                     <  2:2310211/148‑1 (MQ=255)
                                                                                                                aGTTCTCCGCAACCTTCGACGATCAGATGCTGGTGGATTACTCCAAAAACCGCATCACTTAAGAGACGCTGGCGAAATTACAGGATCTGGCGAAAGAGTGCGATCTGGCGGGCGCGATTAAGTCGATGTTCTCTGGCGAGAAGATCaa                                     <  2:212331/148‑1 (MQ=255)
                                                                                                                aGTTCTCCGCAACCTTCGACGATCAGATGCTGGTGGATTACTCCAAAAACCGCATCACTGAAGAGACGCTGGCGAAATTACAGGATCTGGCGAAAGAGTGCGATCTGGCGGGCGCGATTAAGTCGATGTTCTCTGGCGAGAAGATCaa                                     <  2:2195824/148‑1 (MQ=255)
                                                                                                                aGTTCTCCGCAACCTTCGACGATCAGATGCTGGTGGATTACTCCAAAAACCGCATCACTGAAGAGACGCTGGCGAAATTACAGGATCTGGCGAAAGAGTGCGATCTGGCGGGCGCGATTAAGTCGATGTTCTCTGGCGAGAAGATCaa                                     <  2:2205097/148‑1 (MQ=255)
                                                                                                                aGTTCTCCGCAACCTTCGACGATCAGATGCTGGTGGATTACTCCAAAAACCGCATCACTGAAGAGACGCTGGCGAAATTACAGGATCTGGCGAAAGAGTGCGATCTGGCGGGCGCGATTAAGTCGATGTTCTCTGGCGAGAAGATCaa                                     >  2:2224274/1‑148 (MQ=255)
                                                                                                                  ttCTCCGCAACCTTCGACGATCAGATGCTGGTGGATTACTCCAAAAACCGCATCACTGAAGAGACGCTGGCGAAATTACAGGATCTGGCGAAAGAGTGCGATCTGGCGGGCGCGATTAAGTCGATGTTCTCTGGCGAGAAGATCAAcc                                   >  2:1785956/1‑148 (MQ=255)
                                                                                                                  ttCTCCGCAACCTTCGACGATCAGATGCTGGTGAATTACTCCAAAAACCGCATCACTGAAGAGACGCTGGCGAAATTACAGGATCTGGCGAAAGAGTGCGATCTGGCGGGCGCGATTAAGTCGATGTTCTCTGGCGAGAAGATCAAcc                                   >  1:47764/1‑148 (MQ=255)
                                                                                                                   tctcCGCAACCTTCGACGATCAGATGCTGGTGGATTACTCCAAAAACCGCATCACTGAAGAGACGCTGGCGAAATTACAGGATCTGGCGAAAGAGTGCGATCTGGCGGGCGCGATTAAGTCGATGTTCTCTGGCGAGAAGATCAACCg                                  >  1:2490821/1‑148 (MQ=255)
                                                                                                                    ctcCGCAACCTTCGACGATCAGATGCTGGTGGATTACTCCAAAAACCGCATCACTGAAGAGACGCTGGCGAAATTACAGGATCTGGCGAAAGAGTGCGATCTGGCGGGCGCGATTAAGTCGATGTTCTCTGGCGAGAAGATCAACCGc                                 >  1:2323152/1‑148 (MQ=255)
                                                                                                                    ctcCGCAACCTTCGACGATCAGATGCTGGTGGATTACTCCAAAAACCGCATCACTGAAGAGACGCTGGCGAAATTACAGGATCTGGCGAAAGAGTGCGATCTGGCGGGCGCGATTAAGTCGATGTTCTCTGGCGAGAAGATCAACCGc                                 <  2:1241698/148‑1 (MQ=255)
                                                                                                                    ctcCGCAACCTTCGACGATCAGATGCTGGTGAATTACTCCAAAAACCGCATCACTGAAGAGACGCTGGCGAAATTACAGGATCTGGCGAAAGAGTGCGATCTGGCGGGCGCGATTAAGTCGATGTTCTCTGGCGAGAAGATCAACCGc                                 >  2:321394/1‑148 (MQ=255)
                                                                                                                      ccGCAACCTTCGACGATCAGATGCTGGTGAATTACTCCAAAAACCGCATCACTGAAGAGACGCTGGCGAAATTACAGGATCTGGCGAAAGAGTGCGATCTGGCGGGCGCGATTAAGTCGATGTTCTCTGGCGAGAAGATCAACCGCAc                               >  1:1314129/1‑148 (MQ=255)
                                                                                                                           aCCTTCGACGATCAGATGCTGGTGGATTACTCCAAAAACCGCATCACTGAAGAGACGCTGGCGAAATTACAGGATCTGGCGAAAGAGTGCGATCTGGCGGGCGCGATTAAGTCGATGTTCTCTGGCGAGAAGATCAACCGCACTGaaa                          >  1:2441792/1‑148 (MQ=255)
                                                                                                                               tCGACGATCAGATGCTGGTGGATTACTCCAAAAACCGCATCACTGAAGAGACGCTGGCGAAATTACAGGATCTGGCGAAAGAGTGCGATCTGGCGGGCGCGATTAAGTCGATGTTCTCTGGCGAGAAGATCAACCGCACTGAAAACcg                      >  1:1912315/1‑148 (MQ=255)
                                                                                                                               tCGACGATCAGATGCTGGTGAATTACTCCAAAAACCGCATCACTGAAGAGACGCTGGCGAAATTACAGGATCTGGCGAAAGAGTGCGATCTGGCGGGCGCGATTAAGTCGATGTTCTCTGGCGAGAAGATCAACCGCACTGAAAACcg                      >  1:533380/1‑148 (MQ=255)
                                                                                                                                cgacgaTCAGATGCTGGTGGATTACTCCAAAAACCGCATCACTGAAGAGACGCTGGCGAAATTACAGGATCTGGCGAAAGAGTGCGATCTGGCGGGCGCGATTAAGTCGATGTTCTCTGGCGAGAAGATCAACCGCACTGAAAACcgc                     >  2:584558/1‑148 (MQ=255)
                                                                                                                                 gacgtTCAGATGCTGGTGAATTACTCCAAAAACCGCATCACTGAAGAGACGCTGGCGAATTTACTGGATCTGGCGAAATAGTGCGATCTGGCGGGCGCGATTAAGTCGATGTTCTCTGGGGAGATGATCAACCGCACTGAAAACcgcg                    <  2:74833/148‑1 (MQ=255)
                                                                                                                                  acgaTCAGATGCTGGTGGATTACTCCAAAAACCGCATCACTGAAGAGACGCTGGCGAAATTACAGGATCTGGCGAAAGAGTGCGATCTGGCGGGCGCGATTAAGTCGATGTTCTCTGGCGAGAAGATCAACCGCACTGAAAACCgcgc                   >  1:2661450/1‑148 (MQ=255)
                                                                                                                                    gaTCAGATGCTGGTGGATTACTCCAAAAACCGCATCACTGAAGAGACGCTGGCGAAATTACAGGATCTGGCGAANGAGTGCGATCTGGCGGGCGCGATTAAGTCGATGTTCTCTGGCGAGAAGATCAACCGCACTGAAAACCGCGCCg                 <  1:1808505/148‑1 (MQ=255)
                                                                                                                                    gaTCAGATGCTGGTGGATTACTCCAAAAACCGCATCACTGAAGAGACGCTGGCGAAATTACAGGATCTGGCGAAAGAGTGCGATCTGGCGGGCGCGATTAAGTCGATGTTCTCTGGCGAGAAGATCAACCGCACTGAAAACCGCGCCg                 <  2:1228705/148‑1 (MQ=255)
                                                                                                                                    gaTCAGATGCTGGTGGATTACTCCAAAAACCGCATCACTGAAGAGACGCTGGCGAAATTACAGGATCTGGCGAAAGAGTGCGATCTGGCGGGCGCGATTAAGTCGATGTTCTCTGGCGAGAAGATCAACCGCACTGAAAACCGCGCCg                 <  1:3060209/148‑1 (MQ=255)
                                                                                                                                      tCAGATGCTGGTGGATTACTCCAAAAACCGCATCACTGAAGAGACGCTGGCGAAATTACAGGATCTGGCGAAAGAGTGCGATCTGGCGGGCGCGATTAAGTCGATGTTCTCTGGCGAGAAGATCAACCGCACTGAAAACCGCGCCGtg               <  1:2914054/148‑1 (MQ=255)
                                                                                                                                       cAGATGCTGGTGGATTACTCCAAAAACCGCATCACTGAAGAGACGCTGGCGAAATTACAGGATCTGGCGAAAGAGTGCGATCTGGCGGGCGCGATTAAGTCGATGTTCTCTGGCGAGAAGATCAACCGCACTGAAAACCGCGCCGtgc              <  2:924708/148‑1 (MQ=255)
                                                                                                                                         gATGCTGGTGGATTACTCCAAAAACCGCATCACTGAAGAGACGCTGGCGAAATTATAGGATCTGGCGAAAGAGTGCGATCTGGCGGGCGCGATTAAGTCGATGTTCTCTGGCGAGAAGATCAACCGCACTGAAAACCGCGCCGtgctg            >  2:876457/1‑148 (MQ=255)
                                                                                                                                           tGCTGGTGGATTACTCCAAAAACCGCATCACTGAAGAGACGCTGGCGAAATTACAGGATCTGGCGAAAGAGTGCGATCTGGCGGGCGCGATTAAGTCGATGTTCTCTGGCGAGAAGATCTACCGCACTGAAAACCGCGCCGTGCTGCa          <  2:1209337/148‑1 (MQ=255)
                                                                                                                                              tggtggATTACTCCAAAAACCGCATCACTGAAGAGACGCTGTCGAAAGTACAGGATCTGGCGAAAGAGTGCGATCAGGCGCGCGCGATTAAGTCGATGTTCTCTCGCTAGAAGATCAACAGCACTGAAAACAGCGCCGTGATGCACGt       >  1:515516/1‑148 (MQ=255)
                                                                                                                                               ggtggATTACTCCAAAAACCGCATCACTGAAGAGACGCTGGCGAAATTACAGGATCTGGCGAAAGAGTGCGATCTGGCGGGCGCGATTAAGTCGATGTTCTCTGGCGAGAAGATCAACCGCACTGAAAACCGCGCCGTGCTGCACGTa      >  2:2079225/1‑148 (MQ=255)
                                                                                                                                               ggtgaATTACTCCAAAAACCGCATCACTGAAGAGACGCTGGCGAAATTACAGGATCTGGCGAAAGAGTGCGATCTGGCGGGCGCGATTAAGTCGATGTTCTCTGGCGAGAAGATCAACCGCACTGAAAACCGCGCCGTGCTGCACGTa      >  1:2471301/1‑148 (MQ=255)
                                                                                                                                                   gATTACTCCAAAAACCGCATCACTGAAGAGACGCTGGCGAAATTACAGGATCTGGCGAAAGAGTGCGATCTGGCGGGCGCGATTAAGTCGATGTTCTCTGGCGAGAAGATCAACCGCACTGAAAACCGCGCCGTGCTGCACGTAgcgc  >  2:347843/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
AACATCAATCCAACGCAGACCGCTGCCTGGCAGGCACTACAGAAACACTTCGATGAAATGAAAGACGTTACGATCGCCGATCTTTTTGCTAAAGACGGCGATCGTTTTTCTAAGTTCTCCGCAACCTTCGACGATCAGATGCTGGTGGATTACTCCAAAAACCGCATCACTGAAGAGACGCTGGCGAAATTACAGGATCTGGCGAAAGAGTGCGATCTGGCGGGCGCGATTAAGTCGATGTTCTCTGGCGAGAAGATCAACCGCACTGAAAACCGCGCCGTGCTGCACGTAGCGC  >  CP009273/4223692‑4223986

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: