Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,219,550 C→A 16.2% R174L (CGT→CTT)  pepE ← (alpha)‑aspartyl dipeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,219,5500CA16.2% 1100.4 / 125.1 337R174L (CGT→CTT) pepE(alpha)‑aspartyl dipeptidase
Reads supporting (aligned to +/- strand):  ref base C (127/155);  new base A (36/19);  total (163/174)
Fisher's exact test for biased strand distribution p-value = 7.57e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 341 total aligned reads displayed.

CCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCATATCGTTGGTGGTACGAATAGTTGGGCAAGCAAGG  >  CP009273/4219401‑4219699
                                                                                                                                                     |                                                                                                                                                     
ccAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATAc                                                                                                                                                       >  1:4032699/1‑150 (MQ=255)
  aGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACgc                                                                                                                                                     >  1:2161899/1‑150 (MQ=255)
   gCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATAAgct                                                                                                                                                    <  1:339131/150‑1 (MQ=255)
     ggAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACgctgc                                                                                                                                                  >  1:3735704/1‑150 (MQ=255)
      gAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACgctgct                                                                                                                                                 >  1:1698094/1‑150 (MQ=255)
         ccGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATAAGCTGCTCAc                                                                                                                                              >  1:823993/1‑150 (MQ=255)
          cGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACg                                                                                                                                             <  2:1577027/150‑1 (MQ=255)
             ttcttcACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGt                                                                                                                                          <  2:1219127/150‑1 (MQ=255)
               cttcACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGttt                                                                                                                                        >  1:1459343/1‑150 (MQ=255)
                ttcACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTc                                                                                                                                       >  2:439394/1‑150 (MQ=255)
                  cACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCAc                                                                                                                                     >  1:2956410/1‑150 (MQ=255)
                  cACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCAc                                                                                                                                     <  1:4229488/150‑1 (MQ=255)
                  cACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCAc                                                                                                                                     <  1:4229656/150‑1 (MQ=255)
                   aCCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCAcc                                                                                                                                    >  1:3237694/1‑150 (MQ=255)
                     cAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATAAGCTGCTCACGGGTTTCACCtt                                                                                                                                  >  2:436815/1‑150 (MQ=255)
                     cAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATAAGCTGCTCACGGGTTTCACCtt                                                                                                                                  >  2:2754037/1‑150 (MQ=255)
                      aGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCttt                                                                                                                                 >  1:716056/1‑150 (MQ=255)
                      aGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATAAGCTGCTCACGGGTTTCACCttt                                                                                                                                 >  1:1322176/1‑150 (MQ=255)
                       gCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTa                                                                                                                                <  1:1735763/150‑1 (MQ=255)
                        ccTGAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTAt                                                                                                                               <  2:170646/150‑1 (MQ=255)
                              acacATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCtt                                                                                                                         <  2:3805624/150‑1 (MQ=255)
                                   tAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGc                                                                                                                    <  1:4310495/150‑1 (MQ=255)
                                    aaGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCa                                                                                                                   <  2:3518031/150‑1 (MQ=255)
                                     aGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATAAGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAg                                                                                                                  >  1:1769888/1‑150 (MQ=255)
                                        ggtgTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTGTATGGCCTTCCGGCAgcgc                                                                                                               >  1:3496196/1‑150 (MQ=255)
                                           gtTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGtt                                                                                                            <  2:537480/150‑1 (MQ=255)
                                             tCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTgg                                                                                                          <  1:3851072/150‑1 (MQ=255)
                                              cGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGt                                                                                                         <  1:1069931/150‑1 (MQ=255)
                                               gggCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTg                                                                                                        <  1:2102612/150‑1 (MQ=255)
                                               gggCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATAAGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCCGCGCGTTGGTg                                                                                                        >  2:2094740/1‑150 (MQ=255)
                                               gggCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATAAGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTg                                                                                                        >  2:1386843/1‑150 (MQ=255)
                                                 gCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGaa                                                                                                      <  2:664292/150‑1 (MQ=255)
                                                  ccACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAg                                                                                                     >  2:2944854/1‑150 (MQ=255)
                                                   cACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGt                                                                                                    <  1:4588313/150‑1 (MQ=255)
                                                    aCCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTg                                                                                                   <  1:3097384/150‑1 (MQ=255)
                                                       cAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCgg                                                                                                >  2:1314891/1‑150 (MQ=255)
                                                       cAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATAAGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCgg                                                                                                <  2:4109105/150‑1 (MQ=255)
                                                         gcgtagcgtGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGt                                                                                              <  2:122234/150‑1 (MQ=255)
                                                          cgtagcgtGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGtt                                                                                             <  1:3526836/150‑1 (MQ=255)
                                                          cgtagcgtGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGtt                                                                                             <  1:1769273/150‑1 (MQ=255)
                                                             agcgtGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGAt                                                                                          >  1:656763/1‑150 (MQ=255)
                                                              gcgtGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATAAGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGAtt                                                                                         <  2:1950454/150‑1 (MQ=255)
                                                                      tttACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATAAGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCgg                                                                                 <  2:2501434/150‑1 (MQ=255)
                                                                       ttACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCggg                                                                                <  1:3002873/150‑1 (MQ=255)
                                                                        tACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATAAGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGa                                                                               >  1:2005401/1‑150 (MQ=255)
                                                                                   aTCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAgc                                                                    <  2:4032699/150‑1 (MQ=255)
                                                                                    tCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAgcg                                                                   <  2:2559112/150‑1 (MQ=255)
                                                                                       aGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCAt                                                                <  2:4358835/150‑1 (MQ=255)
                                                                                        gTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACTGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATc                                                               <  1:2773254/150‑1 (MQ=255)
                                                                                         ttACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCg                                                              <  1:3061161/150‑1 (MQ=255)
                                                                                         gtACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCg                                                              <  2:2070894/149‑1 (MQ=255)
                                                                                          tACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGa                                                             >  2:1191723/1‑150 (MQ=255)
                                                                                           aCCTTCCGGTAGACCAATACTCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGaa                                                            <  1:1005811/150‑1 (MQ=255)
                                                                                           aCCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGaa                                                            >  1:895995/1‑150 (MQ=255)
                                                                                           aCCTTCCGGGAGACCACTAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGaa                                                            <  2:3132775/150‑1 (MQ=255)
                                                                                             cTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAAc                                                          >  2:3060507/1‑150 (MQ=255)
                                                                                               tCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATAAGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCt                                                        >  2:4642650/1‑150 (MQ=255)
                                                                                                ccGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCtt                                                       <  2:3865327/150‑1 (MQ=255)
                                                                                                     aGACCAATACTCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGa                                                  <  2:3823806/150‑1 (MQ=255)
                                                                                                      gACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGAt                                                 <  2:1065861/150‑1 (MQ=255)
                                                                                                       aCCAATAATCGTCAGTTCTGGCGCGTCGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATc                                                <  2:1058799/150‑1 (MQ=255)
                                                                                                       aCCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATc                                                <  2:3343499/150‑1 (MQ=255)
                                                                                                       aCCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATc                                                <  1:3830076/150‑1 (MQ=255)
                                                                                                        ccAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCg                                               >  1:2911989/1‑150 (MQ=255)
                                                                                                        ccAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCg                                               <  1:1507951/150‑1 (MQ=255)
                                                                                                        ccAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCg                                               >  1:4238652/1‑150 (MQ=255)
                                                                                                         cAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGa                                              >  2:1784406/1‑150 (MQ=255)
                                                                                                            taatCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACaa                                           <  2:3774885/150‑1 (MQ=255)
                                                                                                             actCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAAt                                          <  2:1280362/148‑1 (MQ=255)
                                                                                                             aatCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATAAGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAAt                                          <  1:2596136/150‑1 (MQ=255)
                                                                                                              atCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATc                                         >  2:899780/1‑150 (MQ=255)
                                                                                                              atCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATc                                         >  2:263789/1‑150 (MQ=255)
                                                                                                                 gTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGc                                      >  2:73884/1‑150 (MQ=255)
                                                                                                                   cAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCat                                    >  2:2368374/1‑150 (MQ=255)
                                                                                                                    aGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCata                                   >  1:1052396/1‑150 (MQ=255)
                                                                                                                     gTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCatat                                  >  1:3679527/1‑150 (MQ=255)
                                                                                                                       tCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTGATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCATATCg                                >  2:397656/1‑150 (MQ=255)
                                                                                                                        cTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCATATCGt                               >  1:570653/1‑150 (MQ=255)
                                                                                                                         tGGCGCGACGACCAGCAGTTCGCGAATAAGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCATATCGtt                              >  2:3095125/1‑150 (MQ=255)
                                                                                                                          ggCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCATATCGTtg                             <  2:3253051/150‑1 (MQ=255)
                                                                                                                               gacgacCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCATATCGTTggtggt                        <  1:3291752/150‑1 (MQ=255)
                                                                                                                                acgacCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCATATCGTTGGTGGTa                       <  1:2546090/150‑1 (MQ=255)
                                                                                                                                 cgacCAGCAGTTCGCGAATAAGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGGTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCATATCGTTGGTGGTAc                      >  1:1717262/1‑150 (MQ=255)
                                                                                                                                   acCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCATATCGTTGGTGGTACGa                    >  1:3921114/1‑150 (MQ=255)
                                                                                                                                     cagcagTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCATATCGTTGGTGGTACGAAt                  >  1:2203992/1‑150 (MQ=255)
                                                                                                                                     cagcagTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCATATCGTTGGTGGTACGAAt                  <  2:1733563/150‑1 (MQ=255)
                                                                                                                                     cagcagTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCATATCGTTGGTGGTACGAAt                  <  2:453242/150‑1 (MQ=255)
                                                                                                                                      agcagTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCATATCGTTGGTGGTACGAATa                 <  1:2652912/150‑1 (MQ=255)
                                                                                                                                      agcagTTCGCGAATAAGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCATATCGTTGGTGGTACGAATa                 >  1:4132249/1‑150 (MQ=255)
                                                                                                                                       gcagTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCATATCGTTGGTGGTACGAATAg                <  2:592922/150‑1 (MQ=255)
                                                                                                                                         agTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCATATCGTTGGTGGTACGAATAGtt              >  2:1702858/1‑150 (MQ=255)
                                                                                                                                             cgcgAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCATATCGTTGGTGGTACGAATAGTTGGgc          >  1:1119494/1‑150 (MQ=255)
                                                                                                                                              gcgAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCATATCGTTGGTGGTACGAATAGTTGGgca         >  1:1269990/1‑150 (MQ=255)
                                                                                                                                               cgAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCATATCGTTGGTGGTACGAATAGTTGGgcaa        >  1:1896610/1‑150 (MQ=255)
                                                                                                                                                gAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGAAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCATATCGTTGGTGGTACGAATAGTTGGgcaag       >  2:3187525/1‑150 (MQ=255)
                                                                                                                                                 aaTACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCATATCGTTGGTGGTACGAATAGTTGGgcaagc      <  1:4532836/150‑1 (MQ=255)
                                                                                                                                                  aTAAGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCATATCGTTGGTGGTACGAATAGTTGGgcaagca     <  1:436815/150‑1 (MQ=255)
                                                                                                                                                   tACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCATATCGTTGGTGGTACGAATAGTTGGgcaagcaa    >  2:3454431/1‑150 (MQ=255)
                                                                                                                                                    aCGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCATATCGTTGGTGGTACGAATAGTTGGGcaagcaag   >  1:4076242/1‑150 (MQ=255)
                                                                                                                                                     cGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCATATCGTTGGTGGTACGAATAGTTGGGCAAGCAAgg  >  2:626341/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCATATCGTTGGTGGTACGAATAGTTGGGCAAGCAAGG  >  CP009273/4219401‑4219699

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: