Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,181,504 G→T 100% D69Y (GAT→TAT)  pepT → peptidase T

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,181,5040GT98.6% 684.8 / 0.5 213D69Y (GAT→TAT) pepTpeptidase T
Reads supporting (aligned to +/- strand):  ref base G (2/1);  new base T (86/124);  total (88/125)
Fisher's exact test for biased strand distribution p-value = 5.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 213 total aligned reads displayed.

GTGAGACAGGTTCCCAGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTGGGTATCGGCGATGAAGTTTTATCACCGGTTATG  >  CP009273/1181363‑1181647
                                                                                                                                             |                                                                                                                                               
gTGAGACAGGTTCCCAGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATcc                                                                                                                                           <  1:2083240/148‑1 (MQ=255)
  gagaCAGGTTCCCAGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATACCg                                                                                                                                         >  2:1368361/1‑148 (MQ=255)
     aCAGGTTCCCAGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCg                                                                                                                                      <  1:768623/148‑1 (MQ=255)
      cAGGTTCCCAGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGa                                                                                                                                     <  1:2257515/148‑1 (MQ=255)
      cAGGTTCCCAGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGa                                                                                                                                     <  2:3280517/148‑1 (MQ=255)
        ggTTCCCAGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGAtt                                                                                                                                   >  2:3546723/1‑148 (MQ=255)
         gTTCCCAGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTg                                                                                                                                  <  1:268474/148‑1 (MQ=255)
          ttCCCAGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTgg                                                                                                                                 <  1:2601810/148‑1 (MQ=255)
              cAGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCttt                                                                                                                             >  1:1400615/1‑148 (MQ=255)
               aGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTa                                                                                                                            <  1:2498440/148‑1 (MQ=255)
                 cACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTAtt                                                                                                                          <  1:759213/148‑1 (MQ=255)
                  aCGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTAttt                                                                                                                         <  2:2712072/148‑1 (MQ=255)
                   cGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTtc                                                                                                                        <  1:518547/148‑1 (MQ=255)
                   cGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTtc                                                                                                                        <  1:2084816/148‑1 (MQ=255)
                     gAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTtctc                                                                                                                      <  1:3307874/148‑1 (MQ=255)
                      aaGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCa                                                                                                                     >  2:3245022/1‑148 (MQ=255)
                      aaGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCGTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCGTTATTTCTCa                                                                                                                     >  2:2312728/1‑148 (MQ=255)
                         gCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCAtgt                                                                                                                  >  2:460778/1‑148 (MQ=255)
                         gCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCAtgt                                                                                                                  <  2:1212740/148‑1 (MQ=255)
                            aaTGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGa                                                                                                               <  1:2722785/148‑1 (MQ=255)
                             aTGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGAt                                                                                                              >  1:2182259/1‑148 (MQ=255)
                              tGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTCTATCCCGGCGATTGGCTTTAGTTGTCGTGTgggtg                                                                                                             >  1:3098241/1‑145 (MQ=255)
                               ggAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATAc                                                                                                            <  2:1060051/148‑1 (MQ=255)
                                gAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATAcc                                                                                                           >  1:911100/1‑148 (MQ=255)
                                 aaGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCt                                                                                                          >  1:967892/1‑148 (MQ=255)
                                    ttATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCAc                                                                                                       <  2:3351437/148‑1 (MQ=255)
                                    ttATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCAc                                                                                                       <  2:1589712/148‑1 (MQ=255)
                                    ttATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCAc                                                                                                       <  2:851213/148‑1 (MQ=255)
                                    ttATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATCGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCAc                                                                                                       >  2:2587195/1‑148 (MQ=255)
                                     tATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCAcc                                                                                                      >  2:1195500/1‑148 (MQ=255)
                                        tGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGa                                                                                                   >  1:755828/1‑148 (MQ=255)
                                        tGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATAACTCACCGGa                                                                                                   >  2:2652197/1‑148 (MQ=255)
                                        tGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGAGTGGCTTTATTTCTCATGTGGATACCTCACCGGa                                                                                                   >  1:1136829/1‑148 (MQ=255)
                                         gCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGAt                                                                                                  <  2:766870/148‑1 (MQ=255)
                                         gCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGAt                                                                                                  <  2:1468687/148‑1 (MQ=255)
                                         gCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGAt                                                                                                  <  2:2232916/148‑1 (MQ=255)
                                             ctgctgAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTAt                                                                                                                           >  2:1505314/1‑119 (MQ=255)
                                             ctgctgAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTAt                                                                                                                           <  1:1505314/119‑1 (MQ=255)
                                             ctgctgAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCAGCGGATTGCa                                                                                              >  2:1151681/1‑148 (MQ=255)
                                             ctgctgAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCa                                                                                              >  1:80811/1‑148 (MQ=255)
                                             ctgctgAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCa                                                                                              <  1:689540/148‑1 (MQ=255)
                                              tgctgAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAg                                                                                             <  1:1622023/148‑1 (MQ=255)
                                              tgctgAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAg                                                                                             <  2:527639/148‑1 (MQ=255)
                                              tgctgAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAg                                                                                             <  1:50684/148‑1 (MQ=255)
                                              tgctgAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAg                                                                                             <  1:3110611/148‑1 (MQ=255)
                                                    aaGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCaa                                                                                       >  1:2127718/1‑148 (MQ=255)
                                                         cagcTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAAtg                                                                                  >  1:2609224/1‑148 (MQ=255)
                                                            cTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCCAAAATGTGa                                                                               >  1:1310894/1‑148 (MQ=255)
                                                             tCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGaa                                                                              <  1:1458584/148‑1 (MQ=255)
                                                              cGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAAt                                                                             >  2:2201601/1‑148 (MQ=255)
                                                               gAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATc                                                                            <  2:903985/148‑1 (MQ=255)
                                                               gAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATc                                                                            <  2:161354/148‑1 (MQ=255)
                                                                aaGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGGTTTATTTCTCCTGTGGATACCTCACCCGATTGCAGCGGCAAAAATGTGAATcc                                                                           >  1:1270516/1‑148 (MQ=255)
                                                                aaGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATcc                                                                           >  1:2942084/1‑148 (MQ=255)
                                                                aaGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATcc                                                                           >  1:1881696/1‑148 (MQ=255)
                                                                 agagATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCg                                                                          <  2:1605393/148‑1 (MQ=255)
                                                                    gaTGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCaa                                                                       <  2:3437924/148‑1 (MQ=255)
                                                                    gaTGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCaa                                                                       <  2:1678910/148‑1 (MQ=255)
                                                                       ggggCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCAACGGATTTCAGCGGCAAAAATGTGAATCCGCAAAtt                                                                    >  1:415018/1‑148 (MQ=255)
                                                                        gggCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTAATTCTCATGTGGATACCACACCGGATTTCAGCGGCAAAAATGTGAATCCGCaaaaag                                                                   >  1:1907500/1‑145 (MQ=255)
                                                                         ggCTTCTCAATGGGAGCGTAAGTGAGCTGGGCACTTTGATGGCGACGTTACCGGCTCACGTACCTTGCTATATCCCTGCTATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCACAttgt                                                                  <  2:2027982/148‑1 (MQ=255)
                                                                           cTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAAttgttg                                                                <  2:3127480/148‑1 (MQ=255)
                                                                              aTCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGaaa                                                             >  2:30126/1‑148 (MQ=255)
                                                                              aTCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGaaa                                                             >  1:3551430/1‑148 (MQ=255)
                                                                              aTCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGaaa                                                             <  1:3057644/148‑1 (MQ=255)
                                                                                 aaTGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACt                                                          <  1:3245022/148‑1 (MQ=255)
                                                                                   tgtgACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTAt                                                        <  2:2392624/148‑1 (MQ=255)
                                                                                     tgACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATcg                                                      <  2:1875824/148‑1 (MQ=255)
                                                                                       aCCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGTAAACTATcgcg                                                    >  1:1341841/1‑148 (MQ=255)
                                                                                       aCCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATcgcg                                                    >  2:393317/1‑148 (MQ=255)
                                                                                       aCCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATcgcg                                                    >  2:1943510/1‑148 (MQ=255)
                                                                                        ccTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCgg                                                   >  2:1071898/1‑148 (MQ=255)
                                                                                        ccTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCgg                                                   >  2:103551/1‑148 (MQ=255)
                                                                                          ttAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTg                                                 <  1:3054777/148‑1 (MQ=255)
                                                                                           tAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTgg                                                >  2:663676/1‑148 (MQ=255)
                                                                                             aGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATACCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCg                                              <  1:1368361/148‑1 (MQ=255)
                                                                                              gTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGa                                             <  1:444708/148‑1 (MQ=255)
                                                                                              gTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGa                                             <  1:2354405/148‑1 (MQ=255)
                                                                                               tGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGat                                            <  1:2613143/148‑1 (MQ=255)
                                                                                                 agaAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGatat                                          <  1:48000/148‑1 (MQ=255)
                                                                                                  gaAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATAtt                                         >  2:3469400/1‑148 (MQ=255)
                                                                                                     gggCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTgcg                                      <  1:1750785/148‑1 (MQ=255)
                                                                                                      ggCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTgcgc                                     >  2:2708196/1‑148 (MQ=255)
                                                                                                        cACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTGTCATGTGGATACCTCACCCGATTTCAGCGGCAAAAATCTGAATACGCAAAATGTTGAAAACAATCGCGGTGGCGCTATTGCGCTg                                   >  1:3404809/1‑148 (MQ=255)
                                                                                                         aCTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAACACTGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTgg                                  <  1:2167950/148‑1 (MQ=255)
                                                                                                           tttGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTGGGt                                <  1:2345539/148‑1 (MQ=255)
                                                                                                             tGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTGGGTAt                              >  1:697234/1‑148 (MQ=255)
                                                                                                             tGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTGGGTAt                              <  1:2404057/148‑1 (MQ=255)
                                                                                                               aTGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCAAATGGAAAACTCACCCGAATGCAGCGGCAAAAATGTTAAACCGCAAATTGTTTAAAACTATAGCGGTTGCGATATTGCGCTGGGTATCg                            >  1:3029943/1‑148 (MQ=255)
                                                                                                                  gCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTGGGTATCGGCg                         <  2:2829132/148‑1 (MQ=255)
                                                                                                                  gCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTGGGTATCGGCg                         <  1:1711808/148‑1 (MQ=255)
                                                                                                                       gTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTGGGTATCGGCGATGaa                    <  2:428322/148‑1 (MQ=255)
                                                                                                                         tACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTGGGTATCGGCGATGAAGt                  <  2:94482/148‑1 (MQ=255)
                                                                                                                          aCCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTGGGTATCGGCGATGAAGtt                 <  2:1593474/148‑1 (MQ=255)
                                                                                                                              gCTAACGTCCCTGGCTATATCCCGGCGAGTGGCTTTATTTCTCATGTGGATACCTCACCGGATTTCAGCGGCAAAAATGTGAATACGCAAATTGTTGAAAACTATAGCGGTGGCGCGATTGCGCTGGGTATCGGCGATGAAGTTTTAt             >  1:2125550/1‑148 (MQ=255)
                                                                                                                                tAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTGGGTATCGGCGATGAAGTTTTATCa           <  1:461124/148‑1 (MQ=255)
                                                                                                                                  aCGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTGGGTATCGGCGATGAAGTTTTATCAcc         >  2:765363/1‑148 (MQ=255)
                                                                                                                                    gTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTGGGTATCGGCGATGAAGTTTTATCACCgg       <  1:2252949/148‑1 (MQ=255)
                                                                                                                                        cTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTGGGTATCGGCGATGAAGTTTTATCACCGGTTAt   <  1:2001662/148‑1 (MQ=255)
                                                                                                                                         tGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTGGGTATCGGCGATGAAGTTTTATCACCGGTTATg  <  1:3449943/148‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                               
GTGAGACAGGTTCCCAGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTGGGTATCGGCGATGAAGTTTTATCACCGGTTATG  >  CP009273/1181363‑1181647

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: