Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,947,002 T→A 5.6% L15Q (CTG→CAG)  cmoA → carboxy‑SAM synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,947,0020TA 5.6% 616.9 / 20.8 197L15Q (CTG→CAG) cmoAcarboxy‑SAM synthase
Reads supporting (aligned to +/- strand):  ref base T (108/78);  new base A (4/7);  total (112/85)
Fisher's exact test for biased strand distribution p-value = 2.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

Only 100 of 200 total aligned reads displayed.

GTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTCCCGGATATGATCCAGCGTTCCGTACCCGGCTATTCCAATATTATTTCCATGATTGGTATGTTAGCCGAGCGCTTCGTTCAACCTGGTACGCAGGTTTACGATCTGGG  >  CP009273/1946855‑1947149
                                                                                                                                                   |                                                                                                                                                   
gTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACt                                                                                                                                                     <  2:1556065/148‑1 (MQ=255)
  tGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTgg                                                                                                                                                   >  1:1200432/1‑148 (MQ=255)
   gCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTggg                                                                                                                                                  <  2:727354/148‑1 (MQ=255)
              tGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCt                                                                                                                                       <  1:3477110/148‑1 (MQ=255)
               ggCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCCCTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCtt                                                                                                                                      >  1:643322/1‑148 (MQ=255)
               ggCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCtt                                                                                                                                      >  1:1836719/1‑148 (MQ=255)
                   aaTCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTtgat                                                                                                                                  <  1:1241427/148‑1 (MQ=255)
                    aTCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTtgatg                                                                                                                                 <  2:1478784/148‑1 (MQ=255)
                    aTCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTtgatg                                                                                                                                 >  1:512679/1‑148 (MQ=255)
                    aTCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTtgatg                                                                                                                                 <  1:1564070/148‑1 (MQ=255)
                     tCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTtgatga                                                                                                                                >  1:1095250/1‑148 (MQ=255)
                     tCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTtgatga                                                                                                                                >  1:2555348/1‑148 (MQ=255)
                     tCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTtgatga                                                                                                                                <  1:2477849/148‑1 (MQ=255)
                        ttGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACg                                                                                                                             <  1:1286680/148‑1 (MQ=255)
                         tGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACAGGGCGACTGGACCTTTGATGAACgg                                                                                                                            <  1:637985/148‑1 (MQ=255)
                           gTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGt                                                                                                                          <  1:1698434/148‑1 (MQ=255)
                            tataTGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTa                                                                                                                         <  2:2060184/148‑1 (MQ=255)
                              tatGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGc                                                                                                                       <  2:3207223/148‑1 (MQ=255)
                              tatGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGc                                                                                                                       >  2:323620/1‑148 (MQ=255)
                                tGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTg                                                                                                                     >  1:722803/1‑148 (MQ=255)
                                  cccTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACAGGGCGACTGGACCTTTGATGAACGGGTAGCTGaa                                                                                                                   <  2:329365/148‑1 (MQ=255)
                                    cTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGt                                                                                                                 >  2:1839173/1‑148 (MQ=255)
                                    cTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGt                                                                                                                 <  2:278375/148‑1 (MQ=255)
                                     tGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTCTCGCCAGACAGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTc                                                                                                                <  1:1217876/148‑1 (MQ=255)
                                     tGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTc                                                                                                                >  2:2627869/1‑148 (MQ=255)
                                     tGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTc                                                                                                                <  2:1890937/148‑1 (MQ=255)
                                     tGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTc                                                                                                                <  1:1934671/148‑1 (MQ=255)
                                     tGGGAGTTGGTTTTCTCCCTGCGTGAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTc                                                                                                                <  2:1890983/148‑1 (MQ=255)
                                      gggAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCt                                                                                                               <  2:1726089/148‑1 (MQ=255)
                                       ggAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCtt                                                                                                              >  1:2837829/1‑148 (MQ=255)
                                        gAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTc                                                                                                             >  2:2333957/1‑148 (MQ=255)
                                         aGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTcc                                                                                                            >  2:1866963/1‑148 (MQ=255)
                                          gTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTccc                                                                                                           >  2:2829054/1‑148 (MQ=255)
                                             ggTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTCCCGGa                                                                                                        <  1:2259502/148‑1 (MQ=255)
                                              gTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTCCCGGat                                                                                                       >  1:1995993/1‑148 (MQ=255)
                                               ttttCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTCCCGGata                                                                                                      >  2:3555279/1‑148 (MQ=255)
                                                      ccTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTCCCGGATATGATCCa                                                                                               >  1:2560401/1‑148 (MQ=255)
                                                        tGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTCCCGGATATGATCCAGc                                                                                             >  1:167151/1‑148 (MQ=255)
                                                        tGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTCCCGGATATGATCCAGc                                                                                             >  1:2884472/1‑148 (MQ=255)
                                                        tGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTCCCGGATATGATCCAGc                                                                                             >  1:1111265/1‑148 (MQ=255)
                                                         gCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTCCCGGATATGATCCAGCg                                                                                            <  1:674405/148‑1 (MQ=255)
                                                           gTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTCCCGGATATGATCCAGCGtt                                                                                          <  2:2599486/148‑1 (MQ=255)
                                                             tAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTCCCGGATATGATCCAGCGTTcc                                                                                        >  1:2642165/1‑148 (MQ=255)
                                                             tAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTCCCGGATATGATCCAGCGTTcc                                                                                        >  1:2295477/1‑148 (MQ=255)
                                                                cgcACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTCCCGGATATGATCCAGCGTTCCGTa                                                                                     >  1:3034951/1‑148 (MQ=255)
                                                                  cacaATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTCCCGGATATGATCCAGCGTTCCGTAcc                                                                                   >  2:1243774/1‑148 (MQ=255)
                                                                      aTACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTCCCGGATATGATCCAGCGTTCCGTACCCGGc                                                                               <  2:39307/148‑1 (MQ=255)
                                                                      aTACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTCCCGGATATGATCCAGCGTTCCGTACCCGGc                                                                               <  2:2184797/148‑1 (MQ=255)
                                                                        aCGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTCCCGGATATGATCCAGCGTTCCGTACCCGGCtc                                                                             >  2:3391427/1‑147 (MQ=255)
                                                                         cGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACAGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTCCCGGATATGATCCAGCGTTCCGTACCCGGCTAt                                                                            <  2:2866882/148‑1 (MQ=255)
                                                                          gCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTCCCGGATATGATCCAGCGTTCCGTACCCGGCTAtt                                                                           >  2:3518496/1‑148 (MQ=255)
                                                                          gCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTCCCGGATATGATCCAGCGTTCCGTACCCGGCTAtt                                                                           >  2:1101693/1‑148 (MQ=255)
                                                                            cACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTCCCGGATATGATCCAGCGTTCCGTACCCGGCTATTcc                                                                         >  1:3192214/1‑148 (MQ=255)
                                                                             aCTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTCCCGGATATGATCCAGCGTTCCGTACCCGGCTATTCCa                                                                        >  2:1199269/1‑148 (MQ=255)
                                                                             aCTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTCCCGGATATGATCCAGCGTTCCGTACCCGGCTATTCCa                                                                        >  2:2601066/1‑148 (MQ=255)
                                                                              ctttctttTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTCCCGGATATGATCCAGCGTTCCGTACCCGGCTATTCCaa                                                                       >  1:792054/1‑148 (MQ=255)
                                                                               tttctttTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTCCCGGATATGATCCAGCGTTCCGTACCCGGCTATTCCAat                                                                      >  2:850109/1‑148 (MQ=255)
                                                                                ttctttTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTCCCGGATATGATCCAGCGTTCCGTACCCGGCTATTCCAata                                                                     >  1:8977/1‑148 (MQ=255)
                                                                                 tctttTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTCCCGGATATGATCCAGCGTTCCGTACCCGGCTATTCCAatat                                                                    >  2:200953/1‑148 (MQ=255)
                                                                                 tctttTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTCCCGGATATGATCCAGCGTTCCGTACCCGGCTATTCCAatat                                                                    <  2:643322/148‑1 (MQ=255)
                                                                                 tctttTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTCCCGGATATGATCCAGCGTTCCGTACCCGGCTATTCCAatat                                                                    >  2:370978/1‑148 (MQ=255)
                                                                                   tttttCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTCCCGGATATGATCCAGCGTTCCGTACCCGGCTATTCCAAtatta                                                                  >  2:1794346/1‑148 (MQ=255)
                                                                                    ttttCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTCCCGGATATGATCCAGCGTTCCGTACCCGGCTATTCCAAtattat                                                                 <  1:3376881/148‑1 (MQ=255)
                                                                                    ttttCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTCCCGGATATGATCCAGCGTTCCGTACCCGGCTATTCCAAtattat                                                                 >  1:1755042/1‑148 (MQ=255)
                                                                                    ttttCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTCCCGGATATGATCCAGCGTTCCGTACCCGGCTATTCCAAtattat                                                                 <  1:963000/148‑1 (MQ=255)
                                                                                        cccGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTCCCGGATATGATCCAGCGTTCCGTACCCGGCTATTCCAATATTATTTcc                                                             <  2:1438557/148‑1 (MQ=255)
                                                                                           ggATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTCCCGGATATGATCCAGCGTTCCGTACCCGGCTATTCCAATATTATTTCCATg                                                          <  1:297225/148‑1 (MQ=255)
                                                                                           ggATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTCCCGGATATGATCCAGCGTTCCGTACCCGGCTATTCCAATATTATTTCCATg                                                          <  1:1182876/148‑1 (MQ=255)
                                                                                             aTTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTCCCGGATATGATCCAGCGTTCCGTACCCGGCTATTCCAATATTATTTCCATGAt                                                        <  2:569037/148‑1 (MQ=255)
                                                                                             aTTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTCCCGGATATGATCCAGCGTTCCGTACCCGGCTATTCCAATATTATTTCCATGAt                                                        <  1:390876/148‑1 (MQ=255)
                                                                                               ttttACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTCCCGGATATGATCCAGCGTTCCGTACCCGGCTATTCCAATATTATTTCCATGATTg                                                      >  2:831938/1‑148 (MQ=255)
                                                                                                 ttACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTCCAGGATATGATCCAGCGTTCCGTACCCGGCTATTCCAATATTATTTCCATGATTGGt                                                    <  2:941747/148‑1 (MQ=255)
                                                                                                   aCGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTCCCGGATATGATCCAGCGTTCCGTACCCGGCTATTCCAATATTATTTCCATGATTGGTAt                                                  <  2:1005524/148‑1 (MQ=255)
                                                                                                   aCGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTCCCGGATATGATCCAGCGTTCCGTACCCGGCTATTCCAATATTATTTCCATGATTGGTAt                                                  <  2:2833875/148‑1 (MQ=255)
                                                                                                   aCGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTCCCGGATATGATCCAGCGTTCCGTACCCGGCTATTCCAATATTATTTCCATGATTGGTAt                                                  >  1:2706013/1‑148 (MQ=255)
                                                                                                     gTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTCCCGGATATGATCCAGCGTTCCGTACCCGGCTATTCCAATATTATTTCCATGATTGGTATGt                                                <  1:2829054/148‑1 (MQ=255)
                                                                                                     gTTATGTCTCACCGCGACAAGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAAAGGGTAGCTGAAGTCTTCCCGGATATGATCCAGCGTTCCGTACCCGGCTATTCCAATATTATTTCCATGATTGGTATGt                                                <  2:2766087/148‑1 (MQ=255)
                                                                                                       tATGTCTCACCGCGACAAGCTATTATCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGTTGAACGGGTATCTGAAGTCTTCCCTGATATGATCCAGCGTTACTTAACCGGCTCTTCCAATATTATTTCCATGATTGGTATGTTa                                              >  2:1627158/1‑148 (MQ=255)
                                                                                                        aTGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTCCCGGATATGATCCAGCGTTCCGTACCCGGCTATTCCAATATTATTTCCATGATTGGTATGTTAg                                             <  2:1753655/148‑1 (MQ=255)
                                                                                                            ctcACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTCCCGGATATGATCCAGCGTTCCGTACCCGGCTATTCCAATATTATTTCCATGATTGGTATGTTAGCCGa                                         >  1:161282/1‑148 (MQ=255)
                                                                                                                ccGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTCCCGGATATGATCCAGCGTTCCGTACCCGGCTATTCCAATATTATTTCCATGATTGGTATGTTAGCCGAgcgc                                     <  2:2255402/148‑1 (MQ=255)
                                                                                                                ccGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTCCCGGATATGATCCAGCGTTCCGTACCCGGCTATTCCAATATTATTTCCATGATTGGTATGTTAGCCGAgcgc                                     <  2:466822/148‑1 (MQ=255)
                                                                                                                 cgcgACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTCCCGGATATGATCCAGCGTTCCGTACCCGGCTATTCCAATATTATTTCCATGATTGGTATGTTAGCCGAGCGCt                                    >  1:2784557/1‑148 (MQ=255)
                                                                                                                     acacGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTCCCGGATATGATCCAGCGTTCCGTACCCGGCTATTCCAATATTATTTCCATGATTGGTATGTTAGCCGAGCGCTTCGt                                >  2:1809958/1‑148 (MQ=255)
                                                                                                                            aTTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTCCCGGATATGATCCAGCGTTCCGTACCCGGCTATTCCAATATTATTTCCATGATTGGTATGTTAGCCGAGCGCTTCGTTCAACCt                         >  2:1511132/1‑148 (MQ=255)
                                                                                                                             ttttCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTCCCGGATATGATCCAGCGTTCCGTACCCGGCTATTCCAATATTATTTCCATGATTGGTATGTTAGCAGAGCGCTTCGTTCAACCTg                        >  1:1059695/1‑148 (MQ=255)
                                                                                                                                 cTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTCCCGGATATGATCCAGCGTTCCGTACCCGGCTATTCCAATATTATTTCCATGATTGGTATGTTAGCCGAGCGCTTCGTTCAACCTGGTAc                    >  1:3197277/1‑148 (MQ=255)
                                                                                                                                 cTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTCCCGGATATGATCCAGCGTTCCGTACCCGGCTATTCCAATATTATTTCCATGATTGGTATGTTAGCCGAGCGCTTCGTTCAACCTGGTAc                    >  1:3197529/1‑148 (MQ=255)
                                                                                                                                    ccccTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTCCCGGATATGATCCAGCGTTCCGTACCCGGCTATTCCAATATTATTTCCATGATTGGTATGTTAGCCGAGCGCTTCGTTCAACCTGGTACGCa                 >  1:3142364/1‑148 (MQ=255)
                                                                                                                                     cccTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTCCCGGATATGATCCAGCGTTCCGTACCCGGCTATTCCAATATTATTTCCATGATTGGTATGTTAGCCGAGCGCTTCGTTCAACCTGGTACGCAg                >  1:2337987/1‑148 (MQ=255)
                                                                                                                                       cTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTCCCGGATATGATCCAGCGTTCCGTACCCGGCTATTCCAATATTATTTCCATGATTGGTATGTTAGCCGAGCGCTTCGTTCAACCTGGTACGCAGGt              >  1:2783212/1‑148 (MQ=255)
                                                                                                                                       cTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTCCCGGATATGATCCAGCGTTCCGTACCCGGCTATTCCAATATTATTTCCATGATTGGTATGTTAGCCGAGCGCTTCGTTCAACCTGGTACGCAGGt              <  1:1138283/148‑1 (MQ=255)
                                                                                                                                         aTCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTCCCGGATATGATCCAGCGTTCCGTACCCGGCTATTCCAATATTATTTCCATGATTGGTATGTTAGCCGAGCGCTTCGTTCAACCTGGTACGCAGGttt            >  1:462242/1‑148 (MQ=255)
                                                                                                                                             ccAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTCCCGGATATGATCCAGCGTTCCGTACCCGGCTATTCCAATATTATTTCCATGATTGGTATGTTAGCCGAGCGCTTCGTTCAACCTGGTACGCAGGTTTACGa        <  2:510390/148‑1 (MQ=255)
                                                                                                                                              cAGACAGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTCCCGGATATGATCCAGCGTTCCGTACCCGGCTATTCCAATATTATTTCCATGATTGGTATGTTAGCCGAGCGCTTCGTTCAACCTGGTACGCAGGTTTACGAt       >  1:1776275/1‑148 (MQ=255)
                                                                                                                                               aGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTCCCGGATATGATCCAGCGTTCCGTACCCGGCTATTCCAATATTATTTCCATGATTGGTATGTTAGCCGAGCGCTTCGTTCAACCTGGTACGCAGGTTTACGATc      >  1:1006329/1‑148 (MQ=255)
                                                                                                                                               aGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTCCCGGATATGATCCAGCGTTCCGTACCCGGCTATTCCAATATTATTTCCATGATTGGTATGTTAGCCGAGCGCTTCGTTCAACCTGGTACGCAGGTTTACGATc      >  1:1005743/1‑148 (MQ=255)
                                                                                                                                                gACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTCCCGGATATGATCCAGCGTTCCGTACCCGGCTATTCCAATATTATTTCCATGATTGGTATGTTAGCCGAGCGCTTCGTTCAACCTGGTACGCAGGTTTACGATCt     >  2:3095498/1‑148 (MQ=255)
                                                                                                                                                gACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTCCCGGATATGATCCAGCGTTCCGTACCCGGCTATTCCAATATTATTTCCATGATTGGTATGTTAGCCGAGCGCTTCGTTCAACCTGGTACGCAGGTTTACGATCt     <  2:3034951/148‑1 (MQ=255)
                                                                                                                                                   tGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTCCCGGATATGATCCAGCGTTCCGTACCCGGCTATTCCAATATTATTTCCATGATTGGTATGTTAGCCGAGCGCTTCGTTCAACCTGGTACGCAGGTTTACGATCTggg  <  2:512982/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTCCCGGATATGATCCAGCGTTCCGTACCCGGCTATTCCAATATTATTTCCATGATTGGTATGTTAGCCGAGCGCTTCGTTCAACCTGGTACGCAGGTTTACGATCTGGG  >  CP009273/1946855‑1947149

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: