Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 6,132,601 T→C 12.8% T418A (ACC→GCC)  copA‑II ← copper resistance protein A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE0154516,132,6010TC12.8% 198.9 / 6.1 81T418A (ACC→GCC) copA‑IIcopper resistance protein A
Reads supporting (aligned to +/- strand):  ref base T (36/34);  new base C (4/7);  total (40/41)
Fisher's exact test for biased strand distribution p-value = 5.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

TTCCATGCTGCTGTGGTCCATATCGGCCATGCCGCCTTGGTCCATCCCAGCCATCTTGCTATGGTCCATATCGGCCATGCCGCCTTGGTCCATCCCAGCCATCTTGCTGTGGTCCATGCCGGCCATGCCGCCGTGGTCCATCTCAGCCATTTTGCTGTGATCCATGCCGGCCATGCTGCCATGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGTCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGCCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATACCGGCCATGTTGCCGTGGTCCATGCCGGCCATGCTGCCGTGATCCATACCCATATCACTCATTGAGATGAGAGGGCGCGGATCAACTGCCGGAACTGCAGCTTGCATACCTTCGCGAACGGCTAAAGTGCCTCT  >  AE015451/6132383‑6132826
                                                                                                                                                                                                                          |                                                                                                                                                                                                                                 
ttCCATGCTGCTGTGGTCCATATCGGCCATGCCGCCTTGGTCCATCCCAGCCATCTTGCTATGGTCCATATCGGCCATGCCGCCTTGGTCCATCCCAGCCATCTTGCTGTGGTCCATGCCGGCCATGCCGCCGTGGTCCATCTCAGCCATTTTGCTGTGATCCATGCCGGCCATGCTGCCATGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGTCATGTTGc                                                                                                                                                                                                                           >  1:89182/1‑227 (MQ=255)
     tgctgcTGTGGTCCATATCGGCCATGCCGCCTTGGTCCATCCCAGCCATCTTGCTATGGTCCATATCGGCCATGCCGCCTTGGTCCATCCCAGCCATCTTGCTGTGGTCCATGCCGGCCATGCCGCCGTGGTCCATCTCAGCCATTTTGCTGTGATCCATGCCGGCCATGCTGCCATGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGTCATGTTGCCGTGGTCCATccc                                                                                                                                                                                                              >  1:1104731/1‑235 (MQ=255)
      gctgctGTGGTCCATATCGGCCATGCCGCCTTGGTCCATCCCAGCCATCTTGCTATGGTCCATATCGGCCATGCCGCCTTGGTCCATCCCAGCCATCTTGCTGTGGTCCATGCCGGCCATGCTGCCGTGGTCCATCCCAGCCATTTTGCTGTGATCCATGCCGGCCATGCTGCCATGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGTCATGTTTcc                                                                                                                                                                                                                          >  1:990094/1‑222 (MQ=255)
      gctgctGTGGTCCATATCGGCCATGCCGCCTTGGTCCATCCCAGCCATCTTGCTATGGTCCATATCGGCCATGCCGCCTTGGTCCATCCCAGCCATCTTGCTGTGGTCCATGCCGGCCATGCCGCCGTGGTCCATCTCAGCCATTTTGCTGTGATCCATGCCGGCCATGCTGCCATGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGTCATGTTGCCGTGGTCCATCCCAg                                                                                                                                                                                                            >  1:653990/1‑236 (MQ=255)
             tgGTCCATATCGGCCATGCCGCCTTGGTCCATCCCAGCCATCTTGCTATGGTCCATATCGGCCATGCCGCCTTGGTCCATCCCAGCCATCTTGCTGTGGTCCATGCCGGCCATGCCGCCGTGGTCCATCTCAGCCATTTTGCTGTGATCCATGCCGGCCATGCTGCCATGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGt                                                                                                                                                                                                                                   >  1:881765/1‑206 (MQ=255)
             tgGTCCATATCGGCCATGCCGCCTTGGTCCATCCCAGCCATCTTGCTATGGTCCATATCGGCCATGCCGCCTTGGTCCATCCCAGCCATCTTGCTGTGGTCCATGCCGGCCATGCCGCCGTGGTCCATCTCAGCCATTTTGCTGTGATCCATGCCGGCCATGCTGCCATGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGTCATGTTGCCGTGGTCCATCCCAGCCAttt                                                                                                                                                                                                      >  1:386574/1‑235 (MQ=255)
                      tCGGCCATGCCGCCTTGGTCCATCCCAGCCATCTTGCTATGGTCCATATCGGCCATGCCGCCTTGGTCCATCCCAGCCATCTTGCTGTGGTCCATGCCGGCCATGCCGCCGTGGTCCATCTCAGCCATTTTGCTGTGATCCATGCCGGCCATGCTGCCATGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGTCATGTTGCCGTGGTCCATcc                                                                                                                                                                                                               <  1:960415/217‑1 (MQ=255)
                       cGGCCATGCCGCCTTGGTCCATCCCAGCCATCTTGCTATGGTCCATATCGGCCATGCCGCCTTGGTCCATCCCAGCCATCTTGCTGTGGTCCATGCCGGCCATGCCGCCGTGGTCCATCTCAGCCATTTTGCTGTGATCCATGCCGGCCATGCTGCCATGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGTCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGt                                                                                                                                                                                             >  1:679591/1‑234 (MQ=255)
                            aTGCCGCCTTGGTCCATCCCAGCCATCTTGCTATGGTCCATATCGGCCATGCCGCCTTGGTCCATCCCAGCCATCTTGCTGTGGTCCATGCCGGCCATGCCGCCGTGGTCCATCTCAGCCATTTTGCTGTGATCCATGCCGGCCATGCTGCCATGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGTCATGTTGc                                                                                                                                                                                                                           >  1:1172017/1‑199 (MQ=255)
                                 gccTTGGTCCATCCCAGCCATCTTGCTATGGTCCATATCGGCCATGCCGCCTTGGTCCATCCCAGCCATCTTGCTGTGGTCCATGCCGGCCATGCCGCCGTGGTCCATCTCAGCCATTTTGCTGTGATCCATGCCGGCCATGCTGCCATGGTCCATCCCAGCCATCTTGCTGTGGTCCATTCCGGTCATGTTGc                                                                                                                                                                                                                           >  1:85196/1‑194 (MQ=255)
                                      ggTCCATCCCAGCCATCTTGCTATGGTCCATATCGGCCATGCCGCCTTGGTCCATCCCAGCCATCTTGCTGTGGTCCATGCCGGCCATGCCGCCGTGGTCCATCTCAGCCATTTTGCTGTGATCCATGCCGGCCATGCTGCCATGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGTc                                                                                                                                                                                                                                  <  1:172654/182‑1 (MQ=255)
                                      ggTCCATCCCAGCCATCTTGCTATGGTCCATATCGGCCATGCCGCCTTGGTCCATCCCAGCCATCTTGCTGTGGTCCATGCCGGCCATGCCGCCGTGGTCCATCTCAGCCATTTTGCTGTGATCCATGCCGGCCATGCTGCCATGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGCCATGTTGCCGTGGTCCa                                                                                                                                                                                                                  >  1:265240/1‑198 (MQ=255)
                                                     tCTTGCTATGGTCCATATCGGCCATGCCGCCTTGGTCCATCCCAGCCATCTTGCTGTGGTCCATGCCGGCCATGCCGCCGTGGTCCATCTCAGCCATTTTGCTGTGATCCATGCCGGCCATGCTGCCATGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGTCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGc                                                                                                                                                                                   >  1:485216/1‑214 (MQ=255)
                                                        tGCTATGGTCCATATCGGCCATGCCGCCTTGGTCCATCCCAGCCATCTTGCTGTGGTCCATGCCGGCCATGCCGCCGTGGTCCATCTCAGCCATTTTGCTGTGATCCATGCCGGCCATGCTGCCATGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGTc                                                                                                                                                                                                                                  <  1:889716/164‑1 (MQ=255)
                                                           ggtGGTCCATATCGGCCATGCCGCCTTGGTCCATCCCAGCCATCTTGCTGTGGTCCATGCCGGCCATGCCGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGCCATGCTGCCATGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGTCATGTTGCCGTGGTCCATCCCAGCCATTTTGCtgt                                                                                                                                                                                                <  2:485211/193‑1 (MQ=255)
                                                                      tCGGCCATGCCGCCTTGGTCCATCCCAGCCATCTTGCTGTGGTCCATGCCGGCCATGCCGCCGTGGTCCATCTCAGCCATTTTGCTGTGATCCATGCCGGCCATGCTGCCATGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGTCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATGc                                                                                                                                                                                       <  1:84048/193‑1 (MQ=255)
                                                                       cGGCCATGCCGCCTTGGTCCATCCCAGCCATCTTGCTGTGGTCCATGCCGGCCATGCCGCCGTGGTCCATCTCAGCCATTTTGCTGTGATCCATGCCGGCCATGCTGCCATGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGTCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGCCATGTTGCCGTGGTCCATCCCAGCCAtttt                                                                                                                                                     >  1:21831/1‑226 (MQ=255)
                                                                            aTGCCGCCTTGGTCCATCCCAGCCATCTTGCTGTGGTCCATGCCGGCCATGCCGCCGTGGTCCATCTCAGCCATTTTGCTGTGATCCATGCCGGCCATGCTGCCATGGTCCATCCCAGCCATTTTGCTGTGGTCCATACCGGCCATGTTGCCGTGGTCCAt                                                                                                                                                                                                                 >  2:977534/1‑161 (MQ=255)
                                                                             tGCCGCCTTGGTCCATCCCAGCCATCTTGCTGTGGTCCATGCCGGCCATGCCGCCGTGGTCCATCTCAGCCATTTTGCTGTGATCCATGCCGGCCATGCTGCCATGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGTCATGTTGCCGTGGTCCATCCCAGCCAtttt                                                                                                                                                                                                     <  1:937769/172‑1 (MQ=255)
                                                                                 gccTTGGTCCATCCCATCCATCTTGCTGTGGTCCATGCCGGCCATGCCGCCGTGGTCCATCTCAGCCATTTTGCTGTGATCCATGCCGGCCATGCTGCCATGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGTc                                                                                                                                                                                                                                  <  2:905101/139‑1 (MQ=255)
                                                                                 gccTTGGTCCATCCCAGCCATCTTGCTGTGGTCCATGCCGGCCATGCTGCCGTGGTCCATCTCAGCCATTTTGCTGTGATCCATGCCGGCCATGCTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGTCATGTTGCCGTGGTc                                                                                                                                                                                                                    >  2:727396/1‑153 (MQ=255)
                                                                                      ggTCCATCCCATCCATCTTGCTGTGGTCCATGCCGGCCATGCCGCCGTGGTCCATCTCAGCCATTTTGCTGTGGTCCATGCCGGCCATGCTGCCATGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGTCATGTTGCCGTGGTcc                                                                                                                                                                                                                   <  2:458435/149‑1 (MQ=255)
                                                                                          cATCCCAGCCATCTTGCTGTGGTCCATGCCGGCCATGCCGCCGTGGTCCATCTCAGCCATTTTGCTGTGATCCATGCCGGCCATGCTGCCATGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGTCATGTTGCCGTGGTCCATCCCAGCCATTTTGCtgt                                                                                                                                                                                                <  1:817345/164‑1 (MQ=255)
                                                                                                aGCCATCTTGCTGTGGTCCATGCCGGCCATGCCGCCGTGGTCCATCCCAGCCATTTTGCTGTGATCCATGCCGGCCATGCTGCCATGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGCCATGTTGCCGTGGTCCATCCCAGCCAttt                                                                                                                                                                                                      >  1:795173/1‑152 (MQ=255)
                                                                                                     tCTTGCTGTGGTCCATGCCGGCCATGCCGCCGTGGTCCATCTCAGCCATTTTGCTGTGATCCATGCCGGCCATGCTGCCATGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGTCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGCCATGTTg                                                                                                                                                                            >  2:131892/1‑173 (MQ=255)
                                                                                                              ggTCCATGCCGGCCATGCCGCCGGGGTCCATCGCAGCCATTTTGCTGTGATCCAGGCCGGCCATGCTGCCATGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGTCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATGc                                                                                                                                                                                       <  2:97611/153‑1 (MQ=255)
                                                                                                               gTCCATGCCGGCCATGCCGCCGTGGTCCATCTCAGCCATTTTGCTGTGATCCATGCCGGCCATGCTGCCATGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGTCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATACCGGc                                                                                                                                                                                   >  2:661810/1‑156 (MQ=255)
                                                                                                                  cATGCCGGCCATGCCGCCGTGGTCCATCTCAGCCATTTTGCTGTGATCCATGCCTGCCATGCTGCCATGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGTCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGCCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATACCgt                                                                                                                                    >  1:906352/1‑199 (MQ=255)
                                                                                                                  cATGCCGGCCATGCCGCCGTGGTCCATCTCAGCCATTTTGCTGTGATCCATGCCGGCCATGCTGCCATGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGTCATGTTGCCGTGGTCCATCCCAGCCa                                                                                                                                                                                                         >  2:180215/1‑131 (MQ=255)
                                                                                                                  cATGCCGGCCATGCCGCCGTGGTCCATCTCAGCCATTTTGCTGTGATCCATGCCGGCCATGCTGCCATGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGTCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGCCATGTTGCCGTCGTCCATCCCAGCCAtt                                                                                                                                                       >  2:359277/1‑181 (MQ=255)
                                                                                                                     gCCGGCCATGCCGCCGTGGTCCATCTCAGCCATTTTGCTGTGATCCATGCCGGCCATGCTGCCATGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGTCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGCCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATACCGGCCATGTTGCCGTGGTCCATGCCGGCCCtgtctcttatac                                                                                             >  1:347128/1‑226 (MQ=255)
                                                                                                                       cGGCCATGCCGCCGTGGTCCATCTCCGCCATTTTGCTGTGATCCATGCCGGCCATGCTGCCATGGTCCATCCCAGCCATTTTGCCGTGGTCCATGCCGGCCATGTTGCCGTGGTCCACACCAg                                                                                                                                                                                                            >  2:1038041/1‑123 (MQ=14)
                                                                                                                       cGGCCATGCCGCCGTGGTCCATCTCAGCCATTTTGCTGTGATCCATGCCGGCCATGCTGCCATGGTCCATCCCAGCCCTTTTGCTGTGGTCCATGCCGGTCATGATGCCGTGGTCCATcc                                                                                                                                                                                                               >  2:302860/1‑120 (MQ=255)
                                                                                                                        ggCCATGCCGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGCCATGCTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGTCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGCCATGTTGCCGTGGTCCATcc                                                                                                                                                               <  2:1096808/167‑1 (MQ=255)
                                                                                                                                cgccgTCGTCCATCCCATCCATTTTGCTGTGGTCCATGCCGGCCATGCTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGTCATGCTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGc                                                                                                                                                                                   <  2:1197612/139‑1 (MQ=11)
                                                                                                                                 gccgTGGTCCATCTCAGCCATTTTGCTGTGATCCATGCCGGCCATGCTGCCATGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGTCATGTTGCCGTGGTCCATcccg                                                                                                                                                                                                             >  2:563338/1‑111 (MQ=255)
                                                                                                                                      cgTCCATCTCAGCCATTTTGCTGTGTTCCATGCCGGCCATGCTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGTCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGCCATGTTGCCGTGGTCCATCCCAGCCAtttt                                                                                                                                                     <  2:984364/162‑1 (MQ=255)
                                                                                                                                        cccATCCCAGCCATTTTGCTGTGATCCATGCCGGCCATGCTGCCATGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGTCATGTTGCCGTGGTCCATCCCAGCCAtttt                                                                                                                                                                                                     <  2:239014/112‑1 (MQ=255)
                                                                                                                                           acctcAGCCATTTTGCTGTGATCCATGCCGGCCATGCTGCCATGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGTCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGCCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATACCGGCCATGtt                                                                                                                             <  1:276896/180‑1 (MQ=255)
                                                                                                                                               cAGCCATTTTGCTGTGATCCATGCCGGCCATGCTGCCATGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGTCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGCCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATACCGGCCATGTTGCCGTGGTCCATGCCGGCCATGCTGCCGTGATCCATACCCATATCACTCATtga                                                                       >  1:42510/1‑232 (MQ=255)
                                                                                                                                                   cATCTTGCTGTGGTCCATGCCGGCCATGCTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGTCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGCCATGTTGCCGTGGTCCATCCCAGcc                                                                                                                                                          <  2:851676/145‑1 (MQ=255)
                                                                                                                                                    aTTTTGCTGTGATCCATGCCGGCCATGCTGCCATGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGTCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGCCa                                                                                                                                                                                 >  2:895589/1‑121 (MQ=255)
                                                                                                                                                            gtgATCCATGCCGGCCATGCTGCCATGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGTCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGCCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATACCGGcc                                                                                                                                  <  1:423509/160‑1 (MQ=255)
                                                                                                                                                               aTCCATGCCGGCCATGCTGCCATGGTCCATCCCGGCCTTTTTGCTGTGGTCCATGCCGGTCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGCCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATACCGGCCATGTTg                                                                                                                            <  1:325533/163‑1 (MQ=255)
                                                                                                                                                                 ccATGCCGGCCATGCTGCCTTGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGCCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGCCATGTTGCCGTGGTCCATCCCAGCCAtttt                                                                                                                                                     <  2:763290/136‑1 (MQ=255)
                                                                                                                                                                  cATGCCGGCCATGCTGCCATGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGTCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGCCATGtt                                                                                                                                                                             >  1:1030197/1‑111 (MQ=255)
                                                                                                                                                                  cATGCCGGCCATGCTGCCATGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGTCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGCCATGtt                                                                                                                                                                             <  2:1030192/111‑1 (MQ=255)
                                                                                                                                                                   aTGTCGGCCATGCTGCCATGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGTCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGCCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATACCGGCCATGTTGCCGTGGTCCATGCCGGCCATGCTGCCGTGATc                                                                                            <  1:785864/191‑1 (MQ=255)
                                                                                                                                                                       cGGCCATGCTGCCATGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGTCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGCCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATACCGGCCATGTTGCCGTGGTCCATGCCGGCCATGCTGCCg                                                                                                 >  2:667832/1‑182 (MQ=255)
                                                                                                                                                                        ggCCATGCTGCCATGGTCCATCCCATCCATTTTGCTGTGGTCCATGCCGGCCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGCCATGTTGCCGTGGTCCATCCCAGCCAtttt                                                                                                                                                     <  2:1104726/129‑1 (MQ=255)
                                                                                                                                                                        ggCCATGCTGCCATGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGTCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGCCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATACCGGCCATGTTGCCGTGGTCCATGCCGGCCATGCTGCCGTGATCCATAcc                                                                                      <  1:440806/192‑1 (MQ=255)
                                                                                                                                                                        ggCCATGCTGCCATGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGCCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGCCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATACCGGCCATGTTGCCGTGGTCCATGCCGGCCAtgctg                                                                                                    <  2:679586/178‑1 (MQ=255)
                                                                                                                                                                         gCCATGCTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGTCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGCCATGTTGCCGTGGTCCATCCCAGCCAtttt                                                                                                                                                     <  2:884329/128‑1 (MQ=255)
                                                                                                                                                                         gCCATGCTGCCATGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGTCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGcc                                                                                                                                                                                  <  2:1225490/99‑1 (MQ=255)
                                                                                                                                                                            acgctgcCATGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGTCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGCCATGTTGCCGTGGCCCATCCCAGCCATTTTGCTGTGGTCCATACCGGCAGTGTTGCCGTGGTGCATGc                                                                                                               >  2:467509/3‑163 (MQ=255)
                                                                                                                                                                            aTGCTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGTCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGCCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATACCGGCCATGTTGCCGTGGTCCATGCCGGCCAtgct                                                                                                     <  1:769046/173‑1 (MQ=255)
                                                                                                                                                                             tgctgcCATGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGTCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGCCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATACCGGCCATGTTGCCGTGGTCCATGCCGGCCATGCTGCCGTGATCCATACCCATATCACTCATtgagatg                                                                   >  1:666642/1‑206 (MQ=255)
                                                                                                                                                                                 gcCATGGTCCCTCCCAGCCATTTTGCTGTGGTCCATGCCGGTCATGTTGCCGTGGTCCCTCCCAGCCa                                                                                                                                                                                                         >  2:234265/1‑68 (MQ=255)
                                                                                                                                                                                  ccATGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGCCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGCCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATACCGGCCATGtt                                                                                                                             <  2:386569/143‑1 (MQ=255)
                                                                                                                                                                                   cATGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGTCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGCCATGTTGCCGTGGTCCATCTCAGcc                                                                                                                                                          <  1:131895/113‑1 (MQ=255)
                                                                                                                                                                                    gtcgTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGCCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGCCATGTTGCCGTGGTCCATCCCAGCCAtttt                                                                                                                                                     <  2:723894/114‑1 (MQ=255)
                                                                                                                                                                                       gTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGTCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGCCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCAATTCCGGccct                                                                                                                                >  2:1280547/1‑133 (MQ=255)
                                                                                                                                                                                          cATCCCTTCCATGTTGCTGTGGTCCATATCGGCCATGCCGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGCCATGTTGCCGTGGTCCATCCCAGCCAtttt                                                                                                                                                     <  2:577456/111‑1 (MQ=11)
                                                                                                                                                                                           aTCCCAGCCATTTTGCTGTGGTCCATGCCGGTCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGCCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATACCGGCCATGTTGCCGTGGTCCATGCCGGCCATGCTGCCGTGATCCATACCCATATCACTCATTGAGATGAGAGGGCGCGGATCAACTGCCGGAACTGCAGCTTGCATACCtt                        >  1:1104671/1‑235 (MQ=255)
                                                                                                                                                                                               cAGCCATTTTGCTGTGGTCCATGCCGGTCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGCCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATACCGGCCATGTTGCCGTGGTCCATGCCGGCCAtgctgc                                                                                                   >  1:884031/1‑156 (MQ=255)
                                                                                                                                                                                               cAGCCATTTTGCTGTGGTCCATGCCGGTCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGCCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATACCGGCCATGTTGCCGTGGTCCATGCCGGCCATGCTGCCGTGATCCATACCCATATCa                                                                               >  2:1202945/1‑176 (MQ=255)
                                                                                                                                                                                               cAGCCATTTTGCTGTGGTCCATGCCGGTCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGCCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATACCGGCCATGTTGCCGTGGTCCATGACGGCCATGCTGCCGTGACCCATACCCATATCAATCATtgagatg                                                                   >  1:7848/1‑188 (MQ=255)
                                                                                                                                                                                                aGCCATTTTGCTGTGGTCCATGCCGGTCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGCCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATACCGGCCATGTTGCCGTGGTCCATGCCGGCCATGCTGCCGTGATCCATACCCATATCACTCAtt                                                                         <  1:918030/181‑1 (MQ=255)
                                                                                                                                                                                                   cATTTTGCTGTGGTCCATGCCGGTCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGCCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATACCGGCCATGTTGCCGTGGTCCATGCCGGCCAtgct                                                                                                     <  1:1270391/150‑1 (MQ=255)
                                                                                                                                                                                                    aTTTTGCTGTGGTCCATGCCGGTCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGCCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATACCGGCCATGTTGCCGTGGTCCATGCCGGCCATGCTGCCGTGATCCATACCCATATCACTCAt                                                                          >  2:885025/1‑176 (MQ=255)
                                                                                                                                                                                                    aTTTTGCTGTGGTCCATGCCGGTCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCAAGCCGGCCATGTTGACGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATACCGGCCATGTTGCCGTGGTCCATGCCGGCCATGCTGCCGTGATCCATACCCATATCACTc                                                                            >  2:693065/1‑174 (MQ=255)
                                                                                                                                                                                                     ttttGCTGTGGTCCATGCCGGTCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGCCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTTGTCCATAcc                                                                                                                                      >  1:1159818/1‑115 (MQ=255)
                                                                                                                                                                                                     ttttGCTGTGGTCCATGCCGGTCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGCCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATACCGGCCATGTTGCCGTGGTCCATg                                                                                                                >  2:658475/1‑137 (MQ=255)
                                                                                                                                                                                                     ttttGCTGTGGTCCATGCCGGTCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGCCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATACCGGCCATGTTGCCGTGGTCCATGCCGGCCATGCTGCCGTGATCCATACCCATATCACTCATTGAGATGAGAGGGCGCGGATCAACTGCCGGAACTGCAGCTTGCATACCTTCGCGAACGGc              <  1:693070/235‑1 (MQ=255)
                                                                                                                                                                                                     ttttGCTGTGGTCCATGCCGGTCATGCCGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGCCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATACCGGcc                                                                                                                                  <  2:879607/119‑1 (MQ=255)
                                                                                                                                                                                                        tGCTGTGGTCCACGCCGGTCATGTGGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGCCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATACCGGCCATGTTGCCGTGGTCCATGCCGGCCATGCTGCCGTGATCCATACCCATATCACTCATTGAGATGAGAGGGCGCGGATCAACTGCCGGAACt                                        <  1:149870/206‑1 (MQ=255)
                                                                                                                                                                                                             tgGTCCATGCCGGCCATGTTGCCGTGGTCCACCCCAGCCATTTTGCTGTGGTCCATGCCGGCCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATACCGGCCATGTTGCCGt                                                                                                                        <  2:347123/121‑1 (MQ=255)
                                                                                                                                                                                                              ggTCCATGCCGGTCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGCCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATACCGGCCATGTTGCCGTGGTCCATGCCGGCCATGCTGCCGTGATCCATACCCATATCACTCATTGAGATGAGAGGGCGCGGATCAACTGCCGGAACTGCCGCt                                  >  1:973783/1‑206 (MQ=255)
                                                                                                                                                                                                               gTCCATGCCGGTCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGCCATGTTg                                                                                                                                                                            <  1:847355/67‑1 (MQ=255)
                                                                                                                                                                                                                 ccATGCCGGTCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGCCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATACCGGCCATGTTGCCGTGGTCCATGCCGGCCATGCTGCCGTGATCCATACCCATATCACACATTGAGATGAGAGGGCGCGGATCAACTGCCGGAACTGCAGCTTGCATACCTTc                       <  1:1094221/214‑1 (MQ=255)
                                                                                                                                                                                                                      ccGGTCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGCCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATACCGGCCATGTTGCCGTGGTCCATGCCGGCCATGCTGCCGTGATCCATACCCATATCACTCATTGAGATGAGAGGGCGCGGATCAACTGCCGGAACTGCAGCTTGCATACCTTCGCGAACGGCTAAAGTGCctct  <  1:79439/230‑1 (MQ=255)
                                                                                                                                                                                                                         gTCATGTTGCCGTGGTCCATACCAGCCATTTTGCTGTGGTCCATGCCGGCCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATACCGGCCATGTTGCCGTGGTCCATGCCGGCCATGCTGCCGTGATCCATACCCatat                                                                                 <  2:1040004/148‑1 (MQ=255)
                                                                                                                                                                                                                          tCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGCCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATACCGGCCATGTTGCCGTGGTCCATGCCGGCCATGCTGCCGTGATCCATACCCATATCACTCATTGAGATGAGAGGGCGCGGATCAACTGCCGGAACTGCAGCTTGCATACCTTCGCGAACGGCTAAAGTg       <  1:295675/221‑1 (MQ=255)
                                                                                                                                                                                                                          |                                                                                                                                                                                                                                 
TTCCATGCTGCTGTGGTCCATATCGGCCATGCCGCCTTGGTCCATCCCAGCCATCTTGCTATGGTCCATATCGGCCATGCCGCCTTGGTCCATCCCAGCCATCTTGCTGTGGTCCATGCCGGCCATGCCGCCGTGGTCCATCTCAGCCATTTTGCTGTGATCCATGCCGGCCATGCTGCCATGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGTCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGCCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATACCGGCCATGTTGCCGTGGTCCATGCCGGCCATGCTGCCGTGATCCATACCCATATCACTCATTGAGATGAGAGGGCGCGGATCAACTGCCGGAACTGCAGCTTGCATACCTTCGCGAACGGCTAAAGTGCCTCT  >  AE015451/6132383‑6132826

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑