Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 336,125:1 +T 6.7% coding (210/237 nt) PP_0278 → conserved protein of unknown function

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451336,1251.T 6.7% 353.0 / 18.0 89coding (210/237 nt)PP_0278conserved protein of unknown function
Reads supporting (aligned to +/- strand):  ref base . (54/29);  new base T (2/3);  total (57/32)
Fisher's exact test for biased strand distribution p-value = 3.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.85e-01

TCTCCTTCCGTGAATAAGGATTTGCCATGAAAAAAGATGTTTCTGAGCCAATCACAACCGCAGGCGTGGAAACCTTCCTCGCAGCCGTTAACCCATCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGCCCCGGCGATTTCAGAGATTGCGCTGAAAACGGGGGCCGCTTTGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCCGCTTGCGCTAACCACAAG  >  AE015451/335890‑336347
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                              
tctcCTTCCGTGAATAAGGATTTGCCATGAAAAAAGATGTTTCTGAGCCAATCACAACCGCAGGCGTGGAAACCTTCCTCGCAGCCGTTAACCCATCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGAt                                                                                                                                                                                                                                 <  1:1147080/236‑1 (MQ=255)
   ccTTCCGTGAATAAGGATTTGCCATGAAAAAAGATGTTTCTGAGCCAATCACAACCGCAGGCGTGGAAACCTTCCTCGCAGCCGTTAACCCATCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGAttac                                                                                                                                                                                                                              >  1:623592/1‑233 (MQ=255)
             aTAAGGATTTGCCATGAAAAAAGATGTTTCTGAGCCAATCACAACCGCAGGCGTGGAAACCTTCCTCGCAGCCGTTAACCCATCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGAc                                                                                                                                                                                                                     >  1:37988/1‑235 (MQ=255)
               aaGGATTTGCCATGAAAAAAGATGTTTCTGAGCCAATCACAACCGCAGGCGTGGAAACCTTCCTCGCAGCCGTTAACCCATCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATg                                                                                                                                                                                                                  >  1:293951/1‑236 (MQ=255)
                aGGATTTGCCATGAAAAAAGATGTTTCTGAGCCAATCACAACCGCAGGCGTGGAAACCTTCCTCGCAGCCGTTAACCCATCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATgg                                                                                                                                                                                                                 <  1:421461/236‑1 (MQ=255)
                     ttGCCATGAAAAAAGATGTTTCTGAGCCAATCACAACCGCAGGCGTGGAAACCTTCCTCGCAGCCGTTAACCCATCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAg                                                                                                                                                                                                            >  1:485835/1‑236 (MQ=255)
                     ttGCCATGAAAAAAGATGTTTCTGAGCCAATCACAACCGCAGGCGTGGAAACCTTCCTCGCAGCCGTTAACCCATCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAg                                                                                                                                                                                                            <  1:24249/236‑1 (MQ=255)
                               aaaaGATGTTTCTGAGCCAATCACAACCGCAGGCGTGGAAACCTTCCTCGCAGCCGTTAACCCATCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCtg                                                                                                                                                                                                                           >  2:253715/1‑211 (MQ=255)
                                             cgCCAATCACAACCGCAGGCGTGGAAACCTTCCTCGCAGCCGTTAACCCATCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATc                                                                                                                                                                                    <  2:623589/235‑1 (MQ=255)
                                                      caACCGCAGGCGTGGAAACCTTCCTCGCAGCCGTTAACCCATCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGt                                                                                                                                                                           >  1:627398/1‑236 (MQ=255)
                                                       aaCCGCAGGCGTGGAAACCTTCCTCGCAGCCGTTAACCCATCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGt                                                                                                                                                                           <  1:1230416/235‑1 (MQ=255)
                                                            cAGGCGTGGAAACCTTCCTCGCAGCCGTTAACCCATCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCtgat                                                                                                                                                                                                    >  2:814011/1‑205 (MQ=255)
                                                            cAGGCGTGGAAACCTTCCTCGCAGCCGTTAACCCATCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCtgat                                                                                                                                                                                                    <  1:814014/205‑1 (MQ=255)
                                                               gCGTGGAAACCTTCCTCGCAGCCGTTAACCCATCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCAt                                                                                                                                                                  >  1:158525/1‑236 (MQ=255)
                                                                  tGGAAACCTTCCTCGCAGCCGTTAACCCATCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGCAGGACGACCATAAATTCCTGGGGGCAGCCGAttaccc                                                                                                                                                                                                                            >  2:57178/1‑170 (MQ=255)
                                                                    gAAACCTTCCTCGCAGCCGTTAACCCATCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGGTGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATGACCTGTGCGACATGGCCAAg                                                                                                                                                                                                            >  1:384412/1‑189 (MQ=255)
                                                                         cTTCCTCGCAGCCGTTAACCCATCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGa                                                                                                                                                                                                                      <  1:1017789/174‑1 (MQ=255)
                                                                                                 ccTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGG                                                                                                                                   >  2:611364‑M1/1‑233 (MQ=255)
                                                                                                  cTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCa                                                                                                                                                                                                              >  2:246016/1‑157 (MQ=255)
                                                                                                  cTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGGAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGC                                                                                                                                  >  1:144660‑M1/1‑233 (MQ=255)
                                                                                                        ctcCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTACGGCGCATTAGGCAt                                                                                                                                                           >  1:585847/1‑202 (MQ=255)
                                                                                                             aCGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGGGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCG                                                                                                                       <  2:672018‑M1/233‑1 (MQ=255)
                                                                                                             aCGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTT                                                                                                                               >  1:346381‑M1/1‑225 (MQ=255)
                                                                                                               gcgTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGCTGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCt                                                                                                                                                                                                       >  2:61682/1‑151 (MQ=255)
                                                                                                               gcgTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATc                                                                                                                                                                                          <  2:37988/164‑1 (MQ=255)
                                                                                                               gcgTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTc                                                                                                                                                 >  1:1170245/1‑205 (MQ=255)
                                                                                                               gcgTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATgg                                                                                                                  >  1:1071976/1‑236 (MQ=255)
                                                                                                                 gTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATggg                                                                                                                 <  1:820299/235‑1 (MQ=255)
                                                                                                                      gCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGt                                                                                                                                                                           >  1:971719/1‑172 (MQ=255)
                                                                                                                      gCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGt                                                                                                                                                                           <  2:971716/172‑1 (MQ=255)
                                                                                                                            cGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAagcagc                                                                                                     >  1:508311/1‑236 (MQ=255)
                                                                                                                             gTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGcc                                                                                                    >  1:1180015/1‑236 (MQ=255)
                                                                                                                              tCCCCGTCGACGCGGCCTACGAGCATGGCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTG                                                                                                                              <  2:1135992‑M1/209‑1 (MQ=255)
                                                                                                                                cacGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATc                                                                                                                                                                                          <  2:338731/145‑1 (MQ=255)
                                                                                                                                  cGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCCAGGCGCTGATGAACGCTATCGAGATCGCCAAGGGTAAGGCGCa                                                                                                                                                                   >  2:346860/1‑168 (MQ=255)
                                                                                                                                   gTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCtgat                                                                                                                                                                                                    <  2:6744/134‑1 (MQ=255)
                                                                                                                                   gTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGG                                                                                                                                         <  2:1162109‑M1/193‑1 (MQ=255)
                                                                                                                                    gcgACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGG                                                                                                                                         <  2:449620‑M1/191‑1 (MQ=255)
                                                                                                                                          cgGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGccc                                                                                                   <  1:44297/224‑1 (MQ=255)
                                                                                                                                                aCGGGCATGTCTCGATCCTGATGGGCTACAGCAAGCGCCGGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCAT                                                                                                                                               <  2:219195‑M1/174‑1 (MQ=255)
                                                                                                                                                    gCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATggg                                                                                                                 <  2:485832/200‑1 (MQ=255)
                                                                                                                                                         tctcGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGGAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCG                                                                                                                            >  2:265363‑M1/1‑184 (MQ=255)
                                                                                                                                                         tctcGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGAGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCg                                                                                                                                                                                         <  1:38090/123‑1 (MQ=255)
                                                                                                                                                          ctcGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCTCCAAGGGTAAGGCGCAGTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGT                                                                                                                           >  2:317875‑M1/1‑184 (MQ=255)
                                                                                                                                                          ctcGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGTCCTTGCGTCGCGAATGGGCTGCAAAGCAGCCCCGGCgg                                                                                             >  2:623306/1‑213 (MQ=255)
                                                                                                                                                          ctcGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGCCCCGGCGAttt                                                                                          >  1:1225450/1‑217 (MQ=255)
                                                                                                                                                          ctcGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGCCCCGGCGAt                                                                                            >  2:737580/1‑215 (MQ=255)
                                                                                                                                                          ctcGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGCCCCGGCGATTTCCGAGATTGCGCTGAAAAc                                                                       >  1:224719/1‑236 (MQ=255)
                                                                                                                                                          ctcGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGCCCCGGCGATTTCAGAGATTGCGCTGAAAAc                                                                       >  1:644261/1‑236 (MQ=255)
                                                                                                                                                          ctcGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTC                                                                                                                          >  1:561918‑M1/1‑185 (MQ=255)
                                                                                                                                                            cGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGAGAGCGGCCTTGCGgc                                                                                                                          >  2:395848‑M1/1‑181 (MQ=255)
                                                                                                                                                              aTCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGCCCCGGCGATTTCAGAGATTGCGCTGAAAACgggg                                                                   <  1:253715/236‑1 (MQ=255)
                                                                                                                                                                 cTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGgcg                                                                                                                                                                     >  2:66904/1‑135 (MQ=255)
                                                                                                                                                                 cTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGGCGCGA                                                                                                                      >  1:997571‑M1/1‑182 (MQ=255)
                                                                                                                                                                  ttatgGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCTGCCTT                                                                                                                               <  1:1180070‑M1/170‑1 (MQ=255)
                                                                                                                                                                  tGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGCCCCGGCGATTTCAGAGAt                                                                                   >  1:1193520/1‑216 (MQ=255)
                                                                                                                                                                       ggCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGCCCCGGCGATTTCAGAGATTGCGCTGAAAACGGGGGcc                                                                >  1:1211318/1‑230 (MQ=255)
                                                                                                                                                                       ggCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTG                                                                                                                                          >  2:577681‑M1/1‑156 (MQ=255)
                                                                                                                                                                            cATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGCTCGCCAAg                                                                                                                                                                              >  2:102030/1‑115 (MQ=255)
                                                                                                                                                                            cATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGA                                                                                                                      >  2:731271‑M1/1‑171 (MQ=255)
                                                                                                                                                                            cATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGG                                                                                                                                         >  2:1024120‑M1/1‑152 (MQ=255)
                                                                                                                                                                            cATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGT                                                                                                                                           >  2:255978‑M1/1‑150 (MQ=255)
                                                                                                                                                                                     ccTGGTGCGCGAGGGCGATATCGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGCCCCGGCGATTTCAGAGATTGCGCTGAAAACGGGGGCCGCTTTGCGGCCCTt                                                  >  2:638845‑M2/1‑229 (MQ=255)
                                                                                                                                                                                      cgggTGCGCGAGGGCGATATGGAGGACGACCAGAAATTCCTGGGGGCAGCCGAGTACCTGTGCGACATGGCTCAGGCGCTGATGAACGGTATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATggg                                                                                                                 <  2:218605/164‑1 (MQ=255)
                                                                                                                                                                                         gTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCG                                                                                                                            >  2:493182‑M1/1‑152 (MQ=255)
                                                                                                                                                                                          tGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAagcagc                                                                                                     >  1:1132701/1‑174 (MQ=255)
                                                                                                                                                                                           gcgcGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCT                                                                                                                                >  2:607590‑M1/1‑146 (MQ=255)
                                                                                                                                                                                            cgcgAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCaaa                                                                                                          >  2:260614/1‑167 (MQ=255)
                                                                                                                                                                                            cgcgAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGccc                                                                                                   >  2:80892/1‑174 (MQ=255)
                                                                                                                                                                                            cgcgAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGccc                                                                                                   <  1:80892/174‑1 (MQ=255)
                                                                                                                                                                                               gAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGCCCCGGCGATTTCAGAGATTGCGCTGaa                                                                          >  2:617455/1‑196 (MQ=255)
                                                                                                                                                                                               gAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGCCCCGGCGATTTCAGAGATTGCGCTGaa                                                                          <  1:617458/196‑1 (MQ=255)
                                                                                                                                                                                               gAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGCCCCGGCGATTTCAGAGATTGCGCTGAAAACGGGGGCCGCGTTGCGGCCCTttcgcgacgcaaggccg                                  <  1:167104‑M2/236‑18 (MQ=255)
                                                                                                                                                                                                 gggCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGCACCGGCg                                                                                              >  2:56528/1‑174 (MQ=255)
                                                                                                                                                                                                      atatGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGt                                                                                                                                                                           <  1:960993/92‑1 (MQ=255)
                                                                                                                                                                                                      atatGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGt                                                                                                                                                                           >  2:960990/1‑92 (MQ=255)
                                                                                                                                                                                                            aggaCGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGCCCCGGCGATTTCAGAGATTGCGCTGAAAACGGGGGCCGCTTTGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGc                     <  1:291892/236‑1 (MQ=255)
                                                                                                                                                                                                                cgacCATAAATTCCTGGGTGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGCGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGCCCCGGCGAttt                                                                                          >  2:1219659/1‑163 (MQ=255)
                                                                                                                                                                                                                cgacCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGg                                                                                                                           >  2:447186‑M1/1‑129 (MQ=255)
                                                                                                                                                                                                                    cATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGTAGCGGCCGTGC                                                                                                                             >  2:606493‑M1/1‑124 (MQ=255)
                                                                                                                                                                                                                    cATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGTAGC                                                                                                                                     >  1:11518‑M1/1‑116 (MQ=255)
                                                                                                                                                                                                                    cATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGCCCCGGCGag                                                                                            >  1:46967/1‑156 (MQ=255)
                                                                                                                                                                                                                    cATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGCCCCGGCGATTTCAGAGATTGCGCTGAAAACGGGGGCCGCTTTGCGGCCCTTTCGCGACGCCAGGCCGCTCCCACaaaa                       >  1:978843/1‑226 (MQ=255)
                                                                                                                                                                                                                    cATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGG                                                                                                                                         <  1:815767‑M1/112‑1 (MQ=255)
                                                                                                                                                                                                                    cATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGG                                                                                                                                         >  2:815764‑M1/1‑112 (MQ=255)
                                                                                                                                                                                                                    cATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGCCCCGGCGATTTCAGAGATTGCGCTGAAAACGGGGGCCGCTTTGCGGCCCTttcgcgacgcaaggccg                                  <  1:701109‑M2/215‑18 (MQ=255)
                                                                                                                                                                                                                    cATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGCCCCGGCGATTTCAGAGATTGCGCTGAAAACGGGGGCCGCTTTGCGGCCCTttcgcgaagcacggcc                                   <  1:105693‑M2/214‑17 (MQ=255)
                                                                                                                                                                                                                      tAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGCCCCGGCGATTTCAGAGATTGCGCTGAAAACGGGGGCCGCTTTGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCCGCTTgc              <  2:1071973/233‑1 (MQ=255)
                                                                                                                                                                                                                        aaTTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGCCCCGGCGATTTCAGAGATTGCGCTGAAAACGGGGGCCGCTTTGCGGCCCTt                                                  >  1:445524‑M2/1‑194 (MQ=255)
                                                                                                                                                                                                                                ggggCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGCGCGGCCTTGCGGCGCGAAGGGGCTg                                                                                                              >  2:213515/1‑127 (MQ=255)
                                                                                                                                                                                                                                ggggCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGCCCCGGCGATTTCAGAGATTGCGCTGAAAACGGGGGCCg                                                               >  1:697624/1‑174 (MQ=255)
                                                                                                                                                                                                                                 gggCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGGATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGCCCCGGCGATTTCAGAGATTGCGCTGAAAACGGGGGCCGCTTTGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCCGCTTGCGCTAACCACAAg  <  2:1211315/234‑1 (MQ=255)
                                                                                                                                                                                                                                        cgattACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGCCCCGGCGATTTCAGAGATTGCGCTGAAAACGGGGGCCGCTTTGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCCGCTTGCGCTAACca      >  1:186972/5‑223 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                              
TCTCCTTCCGTGAATAAGGATTTGCCATGAAAAAAGATGTTTCTGAGCCAATCACAACCGCAGGCGTGGAAACCTTCCTCGCAGCCGTTAACCCATCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGCCCCGGCGATTTCAGAGATTGCGCTGAAAACGGGGGCCGCTTTGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCCGCTTGCGCTAACCACAAG  >  AE015451/335890‑336347

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: