Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,219,527 C→A 7.1% A182S (GCG→TCG)  pepE ← (alpha)‑aspartyl dipeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,219,5270CA 7.1% 1097.7 / 43.9 286A182S (GCG→TCG) pepE(alpha)‑aspartyl dipeptidase
Reads supporting (aligned to +/- strand):  ref base C (139/126);  new base A (17/4);  total (156/130)
Fisher's exact test for biased strand distribution p-value = 1.22e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 310 total aligned reads displayed.

ACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCATATCGTTGGTG  >  CP009273/4219379‑4219675
                                                                                                                                                    |                                                                                                                                                    
aCTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGgcgc                                                                                                                                                      >  1:10411/1‑149 (MQ=255)
 cTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGcgcg                                                                                                                                                     >  1:807235/1‑149 (MQ=255)
  ttAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGcga                                                                                                                                                    >  1:1342988/1‑149 (MQ=255)
   tAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGcgac                                                                                                                                                   <  2:1963093/149‑1 (MQ=255)
     aaaaCGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGcgacga                                                                                                                                                 <  2:2342904/149‑1 (MQ=255)
           gTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGAGACGACcagca                                                                                                                                           <  2:3414510/149‑1 (MQ=255)
           gTGACCAGCTTCCAGCGGAAACGCTTCTTCACCACCCATAAACACATAAGTGGTGTTAGGTCCACCCAGCGTAGCGTGACAATTACTCACTGTTATCCACTAACCATCCCGTAGACAAAAAATCGTCAATTCTGTCGCGACCACcagca                                                                                                                                           >  1:3534374/1‑149 (MQ=255)
             gACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGt                                                                                                                                         <  2:437874/149‑1 (MQ=255)
             gACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGAGACGACCAGCAGt                                                                                                                                         >  1:4218136/1‑149 (MQ=255)
              aCCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGtt                                                                                                                                        >  1:4298251/1‑149 (MQ=255)
                  gCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTcgcg                                                                                                                                    >  2:557225/1‑149 (MQ=255)
                      ccAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATa                                                                                                                                <  2:2400182/149‑1 (MQ=255)
                       cAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATAc                                                                                                                               <  2:887738/149‑1 (MQ=255)
                       cAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATAc                                                                                                                               >  1:4288155/1‑149 (MQ=255)
                         gCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACgc                                                                                                                             >  1:2950472/1‑149 (MQ=255)
                         gCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGAGACGACCAGCAGTTCGCGAATACgc                                                                                                                             >  1:4326531/1‑149 (MQ=255)
                           ggAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACgctg                                                                                                                           >  2:1685292/1‑149 (MQ=255)
                            gAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACgctgc                                                                                                                          <  1:3142644/149‑1 (MQ=255)
                             aaCCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACgctgct                                                                                                                         <  1:1182858/149‑1 (MQ=255)
                             aaCCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACgctgct                                                                                                                         >  2:156044/1‑149 (MQ=255)
                             aaCCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACgctgct                                                                                                                         >  2:2759076/1‑149 (MQ=255)
                               ccGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACATTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCa                                                                                                                       >  2:4013309/1‑149 (MQ=255)
                                cGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCAc                                                                                                                      >  2:4957966/1‑149 (MQ=255)
                                    tcttcACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGt                                                                                                                  <  2:1328121/149‑1 (MQ=255)
                                      ttcACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGttt                                                                                                                >  2:4305480/1‑149 (MQ=255)
                                      ttcACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGttt                                                                                                                <  2:4350121/149‑1 (MQ=255)
                                      ttcACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGttt                                                                                                                >  2:4836463/1‑149 (MQ=255)
                                              ccTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTa                                                                                                        >  2:187582/1‑149 (MQ=255)
                                               cTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAAAAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTAt                                                                                                       >  2:1068869/1‑149 (MQ=255)
                                                   aaCACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGcc                                                                                                   <  1:4569976/149‑1 (MQ=255)
                                                    acacATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCt                                                                                                  <  2:451333/149‑1 (MQ=255)
                                                    acacATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCt                                                                                                  <  2:1769138/149‑1 (MQ=255)
                                                     cacaTAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCtt                                                                                                 >  1:2561326/1‑149 (MQ=255)
                                                       caTAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTcc                                                                                               <  1:3828735/149‑1 (MQ=255)
                                                         tAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCgg                                                                                             >  1:1507351/1‑149 (MQ=255)
                                                          aaGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGc                                                                                            >  1:1670200/1‑149 (MQ=255)
                                                           aGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCa                                                                                           >  1:2388137/1‑149 (MQ=255)
                                                           aGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGAGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCa                                                                                           >  1:1012287/1‑149 (MQ=255)
                                                              ggtgTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGAGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAgcg                                                                                        <  2:5201617/149‑1 (MQ=255)
                                                               gtgtTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAgcgc                                                                                       <  1:668181/149‑1 (MQ=255)
                                                                tgtTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGcgcg                                                                                      <  2:710814/149‑1 (MQ=255)
                                                                 gtTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGt                                                                                     <  1:247891/149‑1 (MQ=255)
                                                                  ttCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGtt                                                                                    <  1:38379/149‑1 (MQ=255)
                                                                  ttCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGtt                                                                                    <  1:2642001/149‑1 (MQ=255)
                                                                  ttCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGtt                                                                                    <  1:3275536/149‑1 (MQ=255)
                                                                    cGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTgg                                                                                  >  1:2020804/1‑149 (MQ=255)
                                                                     gggCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGACTACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGt                                                                                 >  2:3244910/1‑149 (MQ=255)
                                                                     gggCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGt                                                                                 >  2:1518479/1‑149 (MQ=255)
                                                                       gCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGa                                                                               <  1:2732220/149‑1 (MQ=255)
                                                                        ccACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGaa                                                                              <  1:1835589/149‑1 (MQ=255)
                                                                          aCCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGt                                                                            >  2:1308784/1‑149 (MQ=255)
                                                                          aCCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGt                                                                            >  2:3656050/1‑149 (MQ=255)
                                                                               gcgtagcgtGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCggg                                                                       >  2:1490249/1‑149 (MQ=255)
                                                                                 gtagcgtGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGtt                                                                     >  2:3060210/1‑149 (MQ=255)
                                                                                   agcgtGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGAGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGa                                                                   >  1:337515/1‑149 (MQ=255)
                                                                                    gcgtGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGAt                                                                  <  1:4073174/149‑1 (MQ=255)
                                                                                    gcgtGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGAt                                                                  >  2:4759721/1‑149 (MQ=255)
                                                                                      gtGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGAttt                                                                >  1:2187719/1‑149 (MQ=255)
                                                                                      gtGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGAttt                                                                <  2:3873154/149‑1 (MQ=255)
                                                                                      gtGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGAttt                                                                >  2:2610146/1‑149 (MQ=255)
                                                                                      gtGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGAttt                                                                >  1:4760316/1‑149 (MQ=255)
                                                                                        gACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGc                                                              >  1:3057941/1‑149 (MQ=255)
                                                                                        gACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGc                                                              <  2:3932775/149‑1 (MQ=255)
                                                                                          ccTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAg                                                            >  1:402672/1‑149 (MQ=255)
                                                                                          ccTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAg                                                            <  2:4672473/149‑1 (MQ=255)
                                                                                            tttACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCg                                                          >  1:3546765/1‑149 (MQ=255)
                                                                                             ttACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCgg                                                         <  1:3206524/149‑1 (MQ=255)
                                                                                                cTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGaa                                                      <  1:1448786/149‑1 (MQ=255)
                                                                                                 tCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAAc                                                     >  1:4372991/1‑149 (MQ=255)
                                                                                                      gtgATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGAtt                                                >  2:4070295/1‑149 (MQ=255)
                                                                                                      gtgATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGAtt                                                >  1:4716690/1‑149 (MQ=255)
                                                                                                       tgATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGAttt                                               <  2:2886608/149‑1 (MQ=255)
                                                                                                        gATCCAGTTACCTTCCGGTAGACCAATACTCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTa                                              <  1:4022208/149‑1 (MQ=255)
                                                                                                         aTCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAg                                             >  2:4816223/1‑149 (MQ=255)
                                                                                                          tCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAgc                                            <  1:4957966/149‑1 (MQ=255)
                                                                                                           ccAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAgcg                                           >  2:335465/1‑149 (MQ=255)
                                                                                                             aGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGAGTTAGCGCa                                         >  2:3503520/1‑149 (MQ=255)
                                                                                                                  ccTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGaa                                    >  1:4259411/1‑149 (MQ=255)
                                                                                                                     tCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAAcc                                 <  1:2845512/149‑1 (MQ=255)
                                                                                                                       cGGTAGACCAATAATCGTCAGTTCTGGCGAGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCtt                               <  2:2704257/149‑1 (MQ=255)
                                                                                                                        ggTAGACCAAGAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTg                              <  1:93388/149‑1 (MQ=255)
                                                                                                                          tAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCg                            <  2:513318/149‑1 (MQ=255)
                                                                                                                          tAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCg                            >  1:3408311/1‑149 (MQ=255)
                                                                                                                          tAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCg                            >  1:808472/1‑149 (MQ=255)
                                                                                                                             aCCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGAt                         <  2:1397976/149‑1 (MQ=255)
                                                                                                                              ccAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATc                        <  1:5177305/149‑1 (MQ=255)
                                                                                                                               cAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCg                       >  1:3645094/1‑149 (MQ=255)
                                                                                                                                 ataatCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATCCGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGAc                     <  2:1834925/149‑1 (MQ=255)
                                                                                                                                  taatCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACa                    >  1:482104/1‑149 (MQ=255)
                                                                                                                                    atCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAAt                  <  1:2816516/149‑1 (MQ=255)
                                                                                                                                    atCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAAt                  <  1:1888857/149‑1 (MQ=255)
                                                                                                                                       gTCAGTTCTGGCGAGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCgg               >  2:1245888/1‑149 (MQ=255)
                                                                                                                                         cAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCa             >  1:4660462/1‑149 (MQ=255)
                                                                                                                                           gTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCata           <  2:2035887/149‑1 (MQ=255)
                                                                                                                                              cTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTAACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACACTCGGCATATCg        >  1:3698571/1‑149 (MQ=255)
                                                                                                                                               tGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCATATCGt       <  2:2947637/149‑1 (MQ=255)
                                                                                                                                                ggCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCATATCGtt      <  1:1753452/149‑1 (MQ=255)
                                                                                                                                                 gcgcGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCATATCGTtg     <  1:4466939/149‑1 (MQ=255)
                                                                                                                                                  cgcgACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCATATCGTtgg    <  2:3791026/149‑1 (MQ=255)
                                                                                                                                                   tcgacgaACAGCAGTTCGCGAATACGCTTCTCACGGGTTTAACCTTTATGTCCTTCCGGCAGAGAGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCATATCGTtggt   <  1:2071806/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
ACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCATATCGTTGGTG  >  CP009273/4219379‑4219675

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: