Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 961,060 G→T 9.9% D429Y (GAC→TAC)  ycaI → ComEC family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273961,0600GT 9.9% 589.8 / 28.5 172D429Y (GAC→TAC) ycaIComEC family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (75/79);  new base T (6/11);  total (82/90)
Fisher's exact test for biased strand distribution p-value = 3.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

Only 100 of 173 total aligned reads displayed.

TTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCTTTGGCTTTACTTTTCTGGGGGTTAAAGAGTTTGCCGGAAGGGTGGATCAACATTGCTGAATGTTGGCAATGGCTATCATTTTCCCCATGGTTCTTACTGGTGGTATGGCGATTAAACGCCTGGCGAACGTTGCCAGCA  >  CP009273/960913‑961206
                                                                                                                                                   |                                                                                                                                                  
     tGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCgg                                                                                                                                               >  1:2983003/1‑148 (MQ=255)
     tGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCgg                                                                                                                                               >  1:2982924/1‑148 (MQ=255)
         aTTAGTCTGACCTCGTCTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCtt                                                                                                                                           <  2:2395698/148‑1 (MQ=255)
               cTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTGATGGTTTCTTGCCGACCGGTCTTTGGCtt                                                                                                                                     >  2:1142128/1‑148 (MQ=255)
                tGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCTTTGGCttt                                                                                                                                    >  2:1785542/1‑148 (MQ=255)
                 gACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCTACCGGTCTTTGGCTTTa                                                                                                                                   >  2:2769607/1‑148 (MQ=255)
                  aCCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCGTGCCGACCGGTCTTTGGCTTTAc                                                                                                                                  <  2:2339542/148‑1 (MQ=255)
                   ccTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCTTTGGCTTTACt                                                                                                                                 >  1:492366/1‑148 (MQ=255)
                     tCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCTTTGGCTTTACttt                                                                                                                               <  1:519580/148‑1 (MQ=255)
                     tCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCTTTGGCTTTACttt                                                                                                                               >  2:930928/1‑148 (MQ=255)
                     tCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCTTTGGCTTTACttt                                                                                                                               <  2:1216179/148‑1 (MQ=255)
                      cGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCTACCGGTCTTTGGCTTTACtttt                                                                                                                              >  2:1182090/1‑148 (MQ=255)
                       gTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCTTTGGCTTTACTTTTc                                                                                                                             <  1:805583/148‑1 (MQ=255)
                        tttATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCTACCGGTCTTTGGCTTTACTTTTCt                                                                                                                            <  1:2908584/148‑1 (MQ=255)
                        tttATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCTTTGGCTTTACTTTTCt                                                                                                                            >  2:2894786/1‑148 (MQ=255)
                        tttATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCTTTGGCTTTACTTTTCt                                                                                                                            >  2:255793/1‑148 (MQ=255)
                        tttATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCTTTGGCTTTACTTTTCt                                                                                                                            >  2:2377904/1‑148 (MQ=255)
                         ttATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCTTTGGCTTTACTTTTCTg                                                                                                                           >  2:2763989/1‑148 (MQ=255)
                           aTTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCTTTGGCTTTACTTTTCTggg                                                                                                                         <  1:825356/148‑1 (MQ=255)
                            ttGCAAATCTTTTTGCTTTTCCCTTTGTGTAATTTATCTCTGTTTCGTTGATGCTGGAGGCGATGGTTGTGTATTTATGCGTGGCGTTATTTCTGGATCATGGGTTATGGTTTCTTGCCGACAGGTCTTTGGCTTTACTTTTCTgggt                                                                                                                        <  1:2084509/148‑2 (MQ=255)
                             tGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCTTTGGCTTTACTTTTCTggggg                                                                                                                       <  2:1759604/148‑1 (MQ=255)
                             tGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTATTGCCGACCGGTCTTTGGCTTTACTTTTCTggggg                                                                                                                       <  1:1887661/148‑1 (MQ=255)
                                aaaTCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCTTTGGCTTTACTTTTCTGGGGGTTa                                                                                                                    <  1:1156110/148‑1 (MQ=255)
                                 aaTCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCTTTGGCTTTACTTTTCTGGGGGTTaa                                                                                                                   >  2:61708/1‑148 (MQ=255)
                                 aaTCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCTTTGGCTTTACTTTTCTGGGGGTTaa                                                                                                                   >  2:61974/1‑148 (MQ=255)
                                   tCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCTTTGGCTTTACTTTTCTGGGGGTTAAag                                                                                                                 <  1:2624891/148‑1 (MQ=255)
                                    cTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCTTTGGCTTTACTTTTCTGGGGGTTAAaga                                                                                                                >  1:1380952/1‑148 (MQ=255)
                                    cTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCTTTGGCTTTACTTTTCTGGGGGTTAAaga                                                                                                                <  1:3212686/148‑1 (MQ=255)
                                    cTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCTTTGGCTTTAATTTTCTGGGGGTTAAaga                                                                                                                >  1:2654951/1‑148 (MQ=255)
                                    cTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATATAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGGTTCTTGCCGACAGGTCTTTTGCTTTACGTTTCTGGGGGTTAAaga                                                                                                                >  1:2001598/1‑148 (MQ=255)
                                      aTTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCTTTGGCTTTACTTTTCTGGGGGTTAAAGAGt                                                                                                              >  2:2104303/1‑148 (MQ=255)
                                       ttAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCTTTGGCTTTACTTTTCTGGGGGTTAAAGAGtt                                                                                                             <  1:730516/148‑1 (MQ=255)
                                       ttAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCTTTGGCTTTACTTTTCTGGGGGTTAAAGAGtt                                                                                                             <  1:494002/148‑1 (MQ=255)
                                         aGCAAGTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCTACCGGTCTTTGGCTTTACTTTTCTGGGGGTTAAAGAGTTTg                                                                                                           <  2:2185060/148‑1 (MQ=255)
                                             aTTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCTTTGGCTTTACTTTTCTGGGGGTTAAAGAGTTTGCCgg                                                                                                       >  2:1486317/1‑148 (MQ=255)
                                             aTTCCCTTGGTGACATATATCACGGTTCCGTTGATCCTCGCCGAGATGCTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCCACCGGTCATAGGCTTTACTTTTCTGGGGGTTAAATAGTTTGCCgg                                                                                                       >  2:29534/1‑148 (MQ=255)
                                               tCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCTACCGGTCTTTGGCTTTACTTTTCTGGGGGTTAAAGAGTTTGCCGGaa                                                                                                     >  1:1021528/1‑148 (MQ=255)
                                               tCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATACTGGAGCAAGGGTTATGGTTTCTTTCCTACCGGTCTTTTGCTTTACTTTTCTGGGGGTTTAATAGTTTTCCGGaa                                                                                                     >  1:1021515/1‑148 (MQ=255)
                                                cccTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCTTTGGCTTTACTTTTCTGGGGGTTAAAGAGTTTGCCGGAAg                                                                                                    >  1:1407385/1‑148 (MQ=255)
                                                cccTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATACTGGAGCAAGGGTTATGGTTTCTTGACGAACGGTCTTTGGCTTTAATTGTCTGGGGGTTAAAGAGTGTTCCGGGAg                                                                                                    >  1:869988/1‑148 (MQ=255)
                                                 ccTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCTACCGGTCTTTGGCTTTACTTTTCTGGGGGTTAAAGAGTTTGCCGGAAgg                                                                                                   <  1:1949083/148‑1 (MQ=255)
                                                 ccTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCTTTGGCTTTACTTTTCTGGGGGTTAAAGAGTTTGCCGGAAgg                                                                                                   >  1:2681990/1‑148 (MQ=255)
                                                  cTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCTTTGGCTTTACTTTTCTGGGGGTTAAAGAGTTTGACGGAAggg                                                                                                  >  1:1895017/1‑148 (MQ=255)
                                                   ttGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCTTTGGCTTTACTTTTCTGGGGGTTAAAGAGTTTGCCGGAAGGGt                                                                                                 <  2:505941/148‑1 (MQ=255)
                                                         aCATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCTTTGGCTTTACTTTTCTGGGGGTTAAAGAGTTTGCCGGAAGGGTGGATCa                                                                                           >  1:545359/1‑148 (MQ=255)
                                                             ttATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCTTTGGCTTTACTTTTCTGGGGGTTAAAGAGTTTGCCGGAAGGGTGGATCAACAt                                                                                       >  1:1717281/1‑148 (MQ=255)
                                                             ttATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCTTTGGCTTTACTTTTCTGGGGGTTAAAGAGTTTGCCGGAAGGGTGGATCAACAt                                                                                       >  1:2264380/1‑148 (MQ=255)
                                                             ttATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCTTTGGCTTTACTTTTCTGGGGGTTAAAGAGTTTGCCGGAAGGGTGGATCAACAt                                                                                       >  2:811825/1‑148 (MQ=255)
                                                                 cACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCTTTGGCTTTACTTTTCTGGGGGTTAAAGAGTTTGCCGGAAGGGTGGATCAACATTGCt                                                                                   >  2:573502/1‑148 (MQ=255)
                                                                      ttCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCTTTGGCTTTACTTTTCTGGGGGTTAAAGAGTTTGTCGGAAGGGTGGATCAACATTGCTGAATg                                                                              >  1:259118/1‑148 (MQ=255)
                                                                        ccGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCTTTGGCTTTACTTTTCTGGGGGTTAAAGAGTTTGCCGGAAGGGTGGATCAACATTGCTGAATGtt                                                                            <  1:2763989/148‑1 (MQ=255)
                                                                          gTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCTTTGGCTTTACTTTTCTGGGGGTTAAAGAGTTTGCCGGAAGGGTGGATCAACATTGCTGAATGTTgg                                                                          >  1:1833725/1‑148 (MQ=255)
                                                                          gTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCTTTGGCTTTACTTTTCTGGGGGTTAAAGAGTTTGCCGGAAGGGTGGATCAACATTGCTGAATGTTgg                                                                          >  2:3000728/1‑148 (MQ=255)
                                                                            tGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCTTTGGCTTTACTTTTCTGGGGGTTAAAGAGTTTGCCGGAAGGGTGGATCAACATTGCTGAATGTTGGCa                                                                        <  1:1061419/148‑1 (MQ=255)
                                                                              aTCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCTTTGGCTTTACTTTTCTGGGGGTTAAAGAGTTTGCCGGAAGGGTGGATCAACATTGCTGAATGTTGGCAAt                                                                      >  2:2618577/1‑148 (MQ=255)
                                                                               tCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCTTTGGCTTTACTTTTCTGGGGGTTAAAGAGTTTGCCGGAAGGGTGGATCAACATTGCTGAATGTTGGCAATg                                                                     <  1:2270853/148‑1 (MQ=255)
                                                                               tCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCTTTGGCTTTACTTTTCTGGGGGTTAAAGAGTTTGCCGGAAGGGTGGATCAACATTGCTGAATGTTGGCAATg                                                                     <  2:1975514/148‑1 (MQ=255)
                                                                                 cTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCTTTGGCTTTACTTTTCTGGGGGTTAAAGAGTTTGCCGGAAGGGTGGATCAACATTGCTGAATGTTGGCAATGGc                                                                   >  2:454638/1‑148 (MQ=255)
                                                                                  tCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCTTTGGCTTTACTTTTCTGGGGGTTAAAGAGTTTGCCGGAAGGGTGGATCAACATTGCTGAATGTTGGCAATGGCt                                                                  >  2:2937812/1‑148 (MQ=255)
                                                                                  tCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGGTATGGTTTCTTGCCGACCGGTCTTTGGCTTTACTTTTCTGGGGGGTAAAGAGTTTGCCGGAAGGGTGGATCAACATTGCTGAATGTTGGCAATGGCt                                                                  >  2:1288165/1‑148 (MQ=255)
                                                                                   cgccgcGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCTTTGGCTTTACTTTTCTGGGGGTTAAAGAGTTTGCCGGAAGGGTGGATCAACATTGCTGAATGTTGGCAATGGCTa                                                                 <  1:302244/148‑1 (MQ=255)
                                                                                    gccgcGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCTACCGGTCTTTGGCTTTACTTTTCTGGGGGTTAAAGAGTTTGCCGGAAGGGTGGATCAACATTGCTGAATGTTGGCAATGGCTAt                                                                <  2:2786206/148‑1 (MQ=255)
                                                                                    gccgcGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCTACCGGTCTTTGGCTTTACTTTTCTGGGGGTTAAAGAGTTTGCCGGAAGGGTGGATCAACATTGCTGAATGTTGGCAATGGCTAt                                                                <  2:1134213/148‑1 (MQ=255)
                                                                                    gccgcGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCTTTGGCTTTACTTTTCTGGGGGTTAAAGAGTTTGCCGGAAGGGTGGATCAACATTGCTGAATGTTGGCAATGGCTAt                                                                <  2:1869737/148‑1 (MQ=255)
                                                                                        cgATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCTTTGGCTTTACTTTTCTGGGGGTTAAAGAGTTTGCCGGAAGGGTGGATCAACATTGCTGAATGTTGGCAATGGCTATCAtt                                                            >  1:2272749/1‑148 (MQ=255)
                                                                                          aTGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCTTTGGCTTTACTTTTCTGGGGGTTAAAGAGTTTGCCGGAAGGGTGGATCAACATTGCTGAATGTTGGCAATGGCTATCAtttt                                                          <  1:3279556/148‑1 (MQ=255)
                                                                                                     tttAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCTTTGGCTTTACTTTTCTGGGGGTTAAAGAGTTTGCCGGAAGGGTGGATCAACATTGCTGAATGTTGGCAATGGCTATCATTTTCCCCATGGTTc                                               <  1:2326565/148‑1 (MQ=255)
                                                                                                     tttAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCTTTGGCTTTACTTTTCTGGGGGTTAAAGAGTTTGCCGGAAGGGTGGATCAACATTGCTGAATGTTGGCAATGGCTATCATTTTCCCCATGGTTc                                               <  1:1947711/148‑1 (MQ=255)
                                                                                                     tttAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCTTTGGCTTTACTTTTCTGGGGGTTAAAGAGTTTGCCGGAAGGGTGGATCAACATTGCTGAATGTTGGCAATGGCTATCATTTTCCCCATGGTTc                                               <  2:2477117/148‑1 (MQ=255)
                                                                                                      ttAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCTTTGGCTTTACTTTTCTGGGGGTTAAAGAGTTTGCCGGAAGGGTGGATCAACATTGCTGAATGTTGGCAATGGCTATCATTTTCCACATGGTTCt                                              >  2:3002714/1‑148 (MQ=255)
                                                                                                       tAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCTTTGGCTTTACTTTTCTGGGGGTTAAAGAGTTTGCCGGAAGGGTGGATCAACATTGCTGAATGTTGGCAATGGCTATAATTTTCCCCATGGTTCtt                                             <  1:3129592/148‑1 (MQ=255)
                                                                                                        aaGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCTTTGGCTTTACTTTTCTGGGGGTTAAAGAGTTTGCCGGAAGGGTGGATCAACATTGCTGAATGTTGGCAATGGCTATCATTTTCCCCATGGTTCTTa                                            <  1:3056171/148‑1 (MQ=255)
                                                                                                        aaGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCTTTGGCTTTACTTTTCTGGGGGTTAAAGAGTTTGCCGGAAGGGTGGATCAACATTGCTGAATGTTGGCAATGGCTATCATTTTCCCCATGGTTCTTa                                            <  1:606698/148‑1 (MQ=255)
                                                                                                         aGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCTTTGGCTTTACTTTTCTGGGGGTTAAAGAGTTTGCCGGAAGGGTGGATCAACATTGCTGAATGTTGGCAATGGCTATCATTTTCCCCATGGTTCTTAc                                           <  2:2121475/148‑1 (MQ=255)
                                                                                                              gCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCTTTGGCTTTACTTTTCTGGGGGTTAAAGAGTTTGCCGGAAGGGTGGATCAACATTGCTGAATGTTGGCAATGGCTATCATTTTCCCCATGGTTCTTACtggtg                                      <  1:1940146/148‑1 (MQ=255)
                                                                                                               ccGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCTACCGGTCTTTGGCTTTACTTTTCTGGGGGTTAAAGAGTTTGCCGGAAGGGTGGATCAACATTGCTGAATGTTGGCAATGGCTATCATTTTCCCCATGGTTCTTACtggtgg                                     >  2:1660252/1‑148 (MQ=255)
                                                                                                               ccGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCTTTGGCTTTACTTTTCTGGGGGTTAAAGAGTTTGCCGGAAGGGTGGATCAACATTGCTGAATGTTGGCAATGGCTATCATTTTCCCCATGGTTCTTACtggtgg                                     <  1:3200721/148‑1 (MQ=255)
                                                                                                                     aTCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCTTTGGCTTTACTTTTCTGGGGGTTAAAGAGTTTGCCGGAAGGGTGGATCAACATTGCTGAATGTTGGCAATGGCTATCATTTTCCCCATGGTTCTTACTGGTGGTATGGc                               >  1:34903/1‑148 (MQ=255)
                                                                                                                       ccTGGAGCAAGGGTTATGGTTTCTTGCCTACCGGTCTTTGGCTTTACTTTTCTGGGGGTTAAAGAGTTTGCCGGAAGGGTGGATCAACATTGCTGAATGTTGGCAATGGCTATCATTTTCCCCATGGTTCTTACTGGTGGTATGGCGa                             <  1:1670448/148‑1 (MQ=255)
                                                                                                                          ggAGCAAGGGTTATGGTTTCTTGCCGACCGGTCTTTGGCTTTACTTTTCTGGGGGTTAAAGAGTTTGCCGGAAGGGTGGATCAACATTGCTGAATGTTGGCAATGGCTATCATTTTCCCCATGGTTCTTACTGGTGGTATGGCGATTa                          >  1:1893123/1‑148 (MQ=255)
                                                                                                                            aGCAAGGGTTATGGTTTCTTGCCGACCGGTCTTTGGCTTTACTTTTCTGGGGGTTAAAGAGTTTGCCGGAAGGGTGGATCAACATTGCTGAATGTTGGCAATGGCTATCATTTTCCCCATGGTTCTTACTGGTGGTATGGCGATTaaa                        >  2:2251036/1‑148 (MQ=255)
                                                                                                                            aGCAAGGGTTATGGTTTCTTGCCGACCGGTCTTTGGCTTTACTTTTCTGGGGGTTAAAGAGTTTGCCGGAAGGGTGGATCAACATTGCTGAATGTTGGCAATGGCTATCATTTTCCCCATGGTTCTTACTGGTGGTATGGCGATTaaa                        <  1:547362/148‑1 (MQ=255)
                                                                                                                                aGGGTTATGGTTTCTTGCCTACCGGTCTTTGGCTTTACTTTTCTGGGGGTTAAAGAGTTTGCCGGAAGGGTGGATCAACATTGCTGAATGTTGGCAATGGCTATCATTTTCCCCATGGTTCTTACTGGTGGTATGGCGATTAAACGcc                    <  1:1776774/148‑1 (MQ=255)
                                                                                                                                aGGGTTATGGTTTCTTGCCGACCGGTCTTTGGCTTTACTTTTCTGGGGGTTAAAGAGTTTGCCGGAAGGGTGGATCAACATTGCTGAATGTTGGCAATGGCTATCATTTTCCCCATGGTTCTTACTGGTGGTATGGCGATTAAACGcc                    <  1:3108456/148‑1 (MQ=255)
                                                                                                                                  ggTTATGGTTTCTTGCCGACCGGTCTTTGGCTTTACTTTTCTGGGGGTTAAAGAGTTTGCCGGAAGGGTGGATCAACATTGCTGAATGTTGGCAATGGCTATCATTTTCCCCATGGTTCTTACTGGTGGTATGGCGATTAAACGCCTg                  >  1:1564540/1‑148 (MQ=255)
                                                                                                                                  ggTTATGGTTTCTTGCCGACCGGTCTTTGGCTTTACTTTTCTGGGGGTTAAAGAGTTTGCCGGAAGGGTGGATCAACATTGCTGAATGTTGGCAATGGCTATCATTTTCCCCATGGTTCTTACTGGTGGTATGGCGATTAAACGCCTg                  >  2:424954/1‑148 (MQ=255)
                                                                                                                                  ggTTATGGTTTCTTGCCGACCGGTCTTTGGCTTTACTTTTCTGGGGGTTAAAGAGTTTGCCGGAAGGGTGGATCAACATTGCTGAATGTTGGCAATGGCTATCATTTTCCCCATGGTTCTTACTGGTGGTATGGCGATTAAACGCCTg                  <  2:2982924/148‑1 (MQ=255)
                                                                                                                                  ggTTATGGGTTCGTGCCGACCGGTCTTTGGCTTTACTTTTCTGGGGGTTAAAGAGTGTGCCGGAAGGGTGGATCAACATTGCTGAATGTTGGCAATGGCGATCATTGTCCCCATGGTTCTTACTGGTGGTATGGCGATTAAACGCCTg                  <  2:3166312/148‑1 (MQ=255)
                                                                                                                                   gTTATGGTTTCTTGCCGACCGGTCTTTGGCTTTACTTTTCTGGGGGTTAAAGAGTTTGCCGGAAGGGTGGATCAACATTGCTGAATGTTGGCAATGGCTATCATTTTCCCCATGGTTCTTACTGGTGGTATGGCGATTAAACGCCTgg                 >  1:261702/1‑148 (MQ=255)
                                                                                                                                   gTTATGGTTTCTTGCCGACCGGTCTTTGGCTTTACTTTTCTGGGGGTTAAAGAGTTTGCCGGAAGGGTGGATCAACATTGCTGAATGTTGGCAATGGCTATCATTTTCCCCATGGTTCTTACTGGTGGTATGGCGATTAAACGCCTgg                 >  1:261624/1‑148 (MQ=255)
                                                                                                                                     tATGGTTTCTTGCCGACCGGTCTTTGGCTTTACTTTTCTGGGGGTTAAAGAGTTTGCCGGAAGGGTGGATCAACATTGCTGAATGTTGGCAATGGCTATCATTTTCCCCATGGTTCTTACTGGTGGTATGGCGATTAAACGCCTGGCg               >  1:2821479/1‑148 (MQ=255)
                                                                                                                                          tttCTTGCCGACCGGTCTTTGTCTTTACTTTTCTGGGGGTTAAAGAGTTTGCCGGAAGGGTGGATCAACATTGCTGAATGTTGGCAATGGCTATCATTTTCCCCATGGTTCTTACTGGTGGTATGGCGATTAAACGCCTGGCGAACGt          <  1:1142128/148‑1 (MQ=255)
                                                                                                                                          tttCTTGCCGACCGGTCTTTGGCTTTACTTTTCTGGGGGTTAAAGAGTTTGCCGGAAGGGTGGATCAACATTGCTGAATGTTGGCAATGGCTATCATTTTCCCCATGGTTCTTACTGGTGGTATGGCGATTAAACGCCTGGCGAACGt          <  1:1141899/148‑1 (MQ=255)
                                                                                                                                            tCTTGCCGACCGGTCTTTGGCTTTACTTTTCTGGGGGTTAAAGAGTTTGCCGGAAGGGTGGATCAACATTGCTGAATGTTGGCAATGGCTATCATTTTCCCCATGGTTCTTACTGGTGGTATGGCGATTAAACGCCTGGCGAACGttt        <  2:326537/148‑2 (MQ=255)
                                                                                                                                            tCTTGCCGACCGGTCTTTGGCTTTACTTTTCTGGGGGTTAAAGAGTTTGCCGGAAGGGTGGATCAACATTGCTGAATGTTGGCAATGGCTATCATTTTCCCCATGGTTCTTACTGGTGGTATGGCGATTAAACGCCTGGCGAACGTTg        <  1:255927/148‑1 (MQ=255)
                                                                                                                                            tCTTGCCGACCGGTCTTTGGCTTTACTTTTCTGGGGGTTAAAGAGTTTGCCGGAAGGGTGGATCAACATTGCTGAATGTTGGCAATGGCTATCATTTTCCCCATGGTTCTTACTGGTGGTATGGCGATTAAACGCCTGGCGAACGTTg        >  2:1261804/1‑148 (MQ=255)
                                                                                                                                            tCTTGCCGACCGGTCTTTGGCTTTACTTTTCTGGGGGTTAAAGAGTTTGCCGGAAGGGTGGATCAACATTGCTGAATGTTGGCAATGGCTATCATTTTCCCCATGGTTCTTACTGGTGGTATGGCGATTAAACGCCTGGCGAACGTTg        <  1:2194874/148‑1 (MQ=255)
                                                                                                                                            tCTTGCCGACCGGTCTTTGGCTTTACTTTTCTGGGGGTTAAAGAGTTTGCCGGAAGGGTGGATCAACATTGCTGAATGTTGGCAATGGCTATCATTTTCCCCATGGTTCTTACTGGTGGTATGGCGATTAAACGCCTGGCGAACGTTg        >  1:3073058/1‑148 (MQ=255)
                                                                                                                                               tGCCGACCGGTCTTTGGCTTTACTTTTCTGGGGGTTAAAGAGTTTGCCGGAAGGGTGGATCAACATTGCTGAATGTTGGCAATGGCTATCATTTTCCCCATGGTTCTTACTGGTGGTATGGCGATTAAACGCCTGGCGAACGTTGCCa     >  1:2434995/1‑148 (MQ=255)
                                                                                                                                                  cGACCGGTCTTTGGCTTTACTTTTCTGGGGGTTAAAGAGTTTGCCGGAAGGGTGGATCAACATTGCTGAATGTTGGCAATGGCTATCATTTTCCCCATGGTTCTTACTGGTGGTATGGCGATTAAACGCCTGGCGAACGTTGCCAGCa  <  2:263706/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCTTTGGCTTTACTTTTCTGGGGGTTAAAGAGTTTGCCGGAAGGGTGGATCAACATTGCTGAATGTTGGCAATGGCTATCATTTTCCCCATGGTTCTTACTGGTGGTATGGCGATTAAACGCCTGGCGAACGTTGCCAGCA  >  CP009273/960913‑961206

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: