Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,219,550 C→A 7.8% 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,5500CA 7.8% 994.3 / 36.3 262R174L (CGT→CTT) pepE(alpha)‑aspartyl dipeptidase
Reads supporting (aligned to +/- strand):  ref base C (113/128);  new base A (10/11);  total (123/139)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 262 total aligned reads displayed.

CAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCATATCGTTGGTGGTACGAATAGTTGGGCAAGCAAG  >  CP009273/4219402‑4219698
                                                                                                                                                    |                                                                                                                                                    
cAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCCGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACg                                                                                                                                                     <  2:1235087/150‑1 (MQ=255)
 aGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACgc                                                                                                                                                    >  2:3528190/1‑150 (MQ=255)
  gCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACgct                                                                                                                                                   <  2:132435/150‑1 (MQ=255)
    ggAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACgctgc                                                                                                                                                 >  2:2032244/1‑150 (MQ=255)
     gAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACgctgct                                                                                                                                                >  1:328754/1‑150 (MQ=255)
            ttcttcACCATCCTTATACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGTCCTTTACTATGTGTGATCCTGTTACCTTCCGGTATACCAATAATAGTCAGTTCTGGCGCGACTACCAGCTGTTCGCGAATAAGATGCTCACTGGt                                                                                                                                         <  2:744712/150‑1 (MQ=255)
                tcACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCa                                                                                                                                     >  2:2991386/1‑150 (MQ=255)
                 cACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCTTCAGTTCAGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCAc                                                                                                                                    >  1:272273/1‑150 (MQ=255)
                     aGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCttt                                                                                                                                >  1:1828116/1‑150 (MQ=255)
                      gCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTa                                                                                                                               <  1:3828597/150‑1 (MQ=255)
                         ttAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATgg                                                                                                                            >  1:1377627/1‑150 (MQ=255)
                         ttAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATgg                                                                                                                            >  1:3485384/1‑150 (MQ=255)
                         ttAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTGTATgg                                                                                                                            >  2:2136586/1‑150 (MQ=255)
                          tAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGc                                                                                                                           >  1:1227222/1‑150 (MQ=255)
                            aaCACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATACTCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCt                                                                                                                         <  2:1886772/150‑1 (MQ=255)
                            aaCACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCt                                                                                                                         <  2:1674805/150‑1 (MQ=255)
                            aaCACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCt                                                                                                                         <  2:4145843/150‑1 (MQ=255)
                            aaCACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATAAGCTGCTCACGGGTTTCACCTTTATGGCCt                                                                                                                         <  2:2260548/150‑1 (MQ=255)
                             acacATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCtt                                                                                                                        >  2:208689/1‑150 (MQ=255)
                              ctcaTAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATACTCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTc                                                                                                                       <  1:4302211/148‑1 (MQ=255)
                              cacaTAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTc                                                                                                                       <  1:540399/150‑1 (MQ=255)
                                caTAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCg                                                                                                                     >  1:2694097/1‑150 (MQ=255)
                                 aTAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCgg                                                                                                                    <  1:57746/150‑1 (MQ=255)
                                   aaGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATAAGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCa                                                                                                                  <  1:2944350/150‑1 (MQ=255)
                                    aGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAg                                                                                                                 <  1:2594846/150‑1 (MQ=255)
                                      tggtgTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAgcg                                                                                                               <  2:1883473/150‑1 (MQ=255)
                                      tggtgTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAgcg                                                                                                               >  2:3344726/1‑150 (MQ=255)
                                        gtgtTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATAAGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGcgcg                                                                                                             <  1:4435077/150‑1 (MQ=255)
                                          gtTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGtt                                                                                                           <  2:4133713/150‑1 (MQ=255)
                                          gtTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGtt                                                                                                           <  1:383336/150‑1 (MQ=255)
                                             cGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGt                                                                                                        <  1:2576811/150‑1 (MQ=255)
                                             cGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGt                                                                                                        >  1:117602/1‑150 (MQ=255)
                                               ggCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGa                                                                                                      <  1:2111955/150‑1 (MQ=255)
                                                gCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGaa                                                                                                     <  2:4606158/150‑1 (MQ=255)
                                                  cACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTGTATGGCCTTCCGGCAGCGCGTTGGTGAAGt                                                                                                   >  2:4534757/1‑150 (MQ=255)
                                                   aCCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTg                                                                                                  >  1:414200/1‑150 (MQ=255)
                                                   aCCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATAAGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTg                                                                                                  >  2:4535028/1‑150 (MQ=255)
                                                       agcgtagcgtGACCTTTACTCACTGTGATCCTGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCggg                                                                                              <  1:1484155/150‑1 (MQ=255)
                                                          gtagcgtGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCGCGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTg                                                                                           <  1:2021090/150‑1 (MQ=255)
                                                           tagcgtGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGa                                                                                          <  1:3429692/150‑1 (MQ=255)
                                                           tagcgtGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGa                                                                                          >  2:2663232/1‑150 (MQ=255)
                                                            agcgtGACCTTTACTCCCTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGTCGACCAGCAGTTCGCTAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGTCAGCGCGTTGTTGAAGTGCGGGTTGAt                                                                                         <  2:4499670/150‑1 (MQ=255)
                                                            agcgtGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGAt                                                                                         >  2:3118730/1‑150 (MQ=255)
                                                            agcgtGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGAt                                                                                         >  2:1086494/1‑150 (MQ=255)
                                                               gtGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTg                                                                                      <  1:663515/150‑1 (MQ=255)
                                                                 gACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCa                                                                                    <  1:4017900/150‑1 (MQ=255)
                                                                 gACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCa                                                                                    <  2:1002568/150‑1 (MQ=255)
                                                                  aCCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAg                                                                                   >  1:1146069/1‑150 (MQ=255)
                                                                         cTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAAc                                                                            >  1:498313/1‑150 (MQ=255)
                                                                         cTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATAAGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAAc                                                                            <  1:3926217/150‑1 (MQ=255)
                                                                           cACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAg                                                                          <  1:2403352/150‑1 (MQ=255)
                                                                           cACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAg                                                                          >  1:1650924/1‑150 (MQ=255)
                                                                           cACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATAAGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAg                                                                          >  2:1265196/1‑150 (MQ=255)
                                                                              tgtgATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGACGTGCGGGTTGATTTGCAGCGGGAACAGAtt                                                                       >  2:479897/1‑150 (MQ=255)
                                                                                tgATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTa                                                                     <  1:4522882/150‑1 (MQ=255)
                                                                                 gATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAg                                                                    >  2:2858706/1‑150 (MQ=255)
                                                                                  aTCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAgc                                                                   <  1:2979051/150‑1 (MQ=255)
                                                                                  aTCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATAAGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAgc                                                                   >  1:738302/1‑150 (MQ=255)
                                                                                   tCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAgcg                                                                  >  1:3694846/1‑150 (MQ=255)
                                                                                     cAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCa                                                                >  1:2842470/1‑150 (MQ=255)
                                                                                        ttACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATAAGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCg                                                             <  2:813540/150‑1 (MQ=255)
                                                                                          aCCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGaa                                                           <  1:4186450/150‑1 (MQ=255)
                                                                                            cTTCCGGTAGACCAATCATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAAc                                                         <  1:801974/150‑1 (MQ=255)
                                                                                             ttCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAAcc                                                        <  2:3851966/150‑1 (MQ=255)
                                                                                                 ggTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGc                                                    >  1:4260328/1‑150 (MQ=255)
                                                                                                  gTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCg                                                   >  1:3182445/1‑150 (MQ=255)
                                                                                                  gTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCg                                                   >  1:2838038/1‑150 (MQ=255)
                                                                                                    aGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGa                                                 >  1:1343504/1‑150 (MQ=255)
                                                                                                      aCCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATc                                               <  1:1631703/150‑1 (MQ=255)
                                                                                                      aCCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATc                                               >  2:2280892/1‑150 (MQ=255)
                                                                                                       ccAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCg                                              <  2:40146/150‑1 (MQ=255)
                                                                                                        cAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGa                                             >  2:4425000/1‑150 (MQ=255)
                                                                                                          ataatCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACa                                           <  1:1503350/150‑1 (MQ=255)
                                                                                                           taatCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACaa                                          <  2:4413711/150‑1 (MQ=255)
                                                                                                             atCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATc                                        <  1:2564058/150‑1 (MQ=255)
                                                                                                              tCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCg                                       <  1:2448141/150‑1 (MQ=255)
                                                                                                              tCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCg                                       <  2:279993/150‑1 (MQ=255)
                                                                                                               cGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCgg                                      <  1:3743184/150‑1 (MQ=255)
                                                                                                                  cAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCat                                   <  1:2474328/150‑1 (MQ=255)
                                                                                                                     ttCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCATATc                                <  2:4239469/150‑1 (MQ=255)
                                                                                                                     ttCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCATATc                                <  2:545891/150‑1 (MQ=255)
                                                                                                                     ttCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCATATc                                <  1:218418/150‑1 (MQ=255)
                                                                                                                     ttCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCATATc                                <  2:1989207/150‑1 (MQ=255)
                                                                                                                         ggCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCATATCGTtg                            <  2:1458661/150‑1 (MQ=255)
                                                                                                                              gacgacCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCATATCGTTggtggt                       >  1:2856581/1‑150 (MQ=255)
                                                                                                                               acgacCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCATATCGTTGGTGGTa                      <  2:3468335/150‑1 (MQ=255)
                                                                                                                               acgacCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCATATCGTTGGTGGTa                      <  2:3744683/150‑1 (MQ=255)
                                                                                                                                 gacCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCATATCGTTGGTGGTACg                    <  2:3162654/150‑1 (MQ=255)
                                                                                                                                 gacCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCATATCGTTGGTGGTACg                    >  1:4320340/1‑150 (MQ=255)
                                                                                                                                 gacCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTGTATGGCCTTCCGGCAGCGCGTTGGTGAATAGCGGGTTGATTTGCAGCGGGAACACATTTAGCGCATCGAAACCTTGCGGAACGACAATCGGCATATCGTTGGTGGTACg                    >  2:666628/1‑150 (MQ=255)
                                                                                                                                    cagcagTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCATATCGTTGGTGGTACGAAt                 <  2:1489393/150‑1 (MQ=255)
                                                                                                                                       cagTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCATATCGTTGGTGGTACGAATAGt              >  1:904541/1‑150 (MQ=255)
                                                                                                                                            cgcgAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCATATCGTTGGTGGTACGAATAGTTGGgc         <  1:2616257/150‑1 (MQ=255)
                                                                                                                                            cgcgAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCATATCGTTGGTGGTACGAATAGTTGGgc         >  1:1034277/1‑150 (MQ=255)
                                                                                                                                               gAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGAGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCATATCGTTGGTGGTACGAATAGTTGGgcaag      <  2:3004016/150‑1 (MQ=255)
                                                                                                                                                 aTACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCATATCGTTGGTGGTACGAATAGTTGGgcaagca    >  1:850443/1‑150 (MQ=255)
                                                                                                                                                 aTACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCATATCGTTGGTGGTACGAATAGTTGGgcaagca    >  1:1291444/1‑150 (MQ=255)
                                                                                                                                                  tACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCATATCGTTGGTGGTACGAATAGTTGGgcaagcaa   <  1:83482/150‑1 (MQ=255)
                                                                                                                                                  tACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCATATCGTTGGTGGTACGAATAGTTGGgcaagcaa   <  1:3573558/150‑1 (MQ=255)
                                                                                                                                                   aCGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCATATCGTTGGTGGTACGAATAGTTGGGcaagcaag  <  1:4549057/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCATATCGTTGGTGGTACGAATAGTTGGGCAAGCAAG  >  CP009273/4219402‑4219698

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: