Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 1,530,160 C→T 100% R395C (CGT→TGT) 
noncoding (49/244 nt)		 ftsZ →
PP_mr19 →		 Cell division protein FtsZ
Xoo8 ncRNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE0154511,530,1600CT98.7% 304.0 / ‑4.3 79R395C (CGT→TGT) 
noncoding (49/244 nt)		ftsZ
PP_mr19		Cell division protein FtsZ
Xoo8 ncRNA
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (26/52);  total (26/53)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.78e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

AGAAGCCTGTGAAGGTGGTCGACAATACCCTGCAGACTGCCCAGCAGGCGTACGAAGCGTCCAATCCTGCGCCGGTGCGCCAGGAGCAGCCAGCGGTCAACTACCGTGACCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTAATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCAATTTGCGCTGAAGCGGCCAATGCCATGATTAAACAACGCACCCTGAAGAATACCATCCGTGCCACAGGTGTCGGTCTGCACTCCGGGGAGAAGGTTTACCTGACCCTCAAGCCTGCACCTGTTGACA  >  AE015451/1529939‑1530389
                                                                                                                                                                                                                             |                                                                                                                                                                                                                                     
agaAGCCTGTGAAGGTGGTCGACAATACCCTGCAGACTGCCCAGCAGGCGTACGAAGCGTCCAATCCTGCGCCGGTGCGCCAGGAGCAGCCAGCGGTCAACTACCGTGACCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTTGTCAGGCTGATTaa                                                                                                                                                                                                                         >  1:691627/1‑236 (MQ=255)
         tgAAGGTGGTCGACAATACCCTGCAGACTGCCCAGCAGGCGTACGAAGCGTCCAATCCTGCGCCGGTGCGCCAGGAGCAGCCAGCGGTCAACTACCGTGACCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTTGTCAGGCTGATTAATGGAAt                                                                                                                                                                                                                   >  2:745178/1‑233 (MQ=255)
           aaGGTGGTCGACAATACCCTGCAGACTGCCCAGCAGGCGTACGAAGCGTCCAATCCTGCGCCGGTGCGCCAGGAGCAGCCAGCGGTCAACTACCGTGACCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTTGTCAGGCTGATTAATGGAATTTAt                                                                                                                                                                                                               <  1:962998/235‑1 (MQ=255)
           aaGGTGGTCGACAATACCCTGCAGACTGCCCAGCAGGCGTACGAAGCGTCCAATCCTGCGCCGGTGCGCCAGGAGCAGCCAGCGGTCAACTACCGTGACCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTTGTCAGGCTGATTAATGGAATTTAt                                                                                                                                                                                                               <  1:745179/235‑1 (MQ=255)
           aaGGTGGTCGACAATACCCTGCAGACTGCCCAGCAGGCGTACGAAGCGTCCAATCCTGCGCCGGTGCGCCAGGAGCAGCCAGCGGTCAACTACCGTGACCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTTGTCAGGCTGATTAATGGAATTTAt                                                                                                                                                                                                               <  1:647324/235‑1 (MQ=255)
           aaGGTGGTCGACAATACCCTGCAGACTGCCCAGCAGGCGTACGAAGCGTCCAATCCTGCGCCGGTGCGCCAGGAGCAGCCAGCGGTCAACTACCGTGACCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTTGTCAGGCTGATTAATGGAATTTAt                                                                                                                                                                                                               <  1:47102/235‑1 (MQ=255)
            aGGTGGTCGACAATACCCTGCAGACTGCCCAGCAGGCGTACGAAGCGTCCAATCCTGCGCCGGTGCGCCAGGAGCAGCCAGCGGTCAACTACCGTGACCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTTGTCAGGCTGATTAATGGAATTTAt                                                                                                                                                                                                               <  2:708182/234‑1 (MQ=255)
                 gtCGACAATACCCTGCAGACTGCCCAGCAGGCGTACGAAGCGTCCAATCCTGCGCCGGTGCGCCAGGAGCAGCCAGCGGTCAACTACCGTGACCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTTGTCAGGCTGATTAATGGAATTTATCAgggg                                                                                                                                                                                                         >  2:499725/1‑235 (MQ=255)
                 gtCGACAATACCCTGCAGACTGCCCAGCAGGCGTACGAAGCGTCCAATCCTGCGCCGGTGCGCCAGGAGCAGCCAGCGGTCAACTACCGTGACCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTTGTCAGGCTGATTAATGGAATTTATCAGGGGt                                                                                                                                                                                                        <  1:272804/236‑1 (MQ=255)
                 ggcGACAATACCCTGCGGACTGCCCAGCAGGCGTACGAAGCGTCCAATCCTGCGCCGGTGCGCCAGGAGCAGCCAGCGGTCAACTACCGTGACCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTTGTCAGGCTGAtt                                                                                                                                                                                                                           <  2:999695/215‑1 (MQ=255)
                  tCGACAATACCCTGCAGACTGCCCAGCAGGCGTACGAAGCGTCCAATCCTGCGCCGGTGCGCCAGGAGCAGCCAGCGGTCAACTACCGTGACCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTTGTCAGGCTGATTAATGGAATTTATCAGGGGta                                                                                                                                                                                                       >  1:1018273/1‑236 (MQ=255)
                   cGACAATACCCTGCAGACTGCCCAGCAGGCGTACGAAGCGTCCAATCCTGCGCCGGTGCGCCAGGAGCAGCCAGCGGTCAACTACCGTGACCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTTGTCAGGCTGATTAATGGAATTTATCAGGGGt                                                                                                                                                                                                        >  2:907538/1‑234 (MQ=255)
                              tGCAGACTGCCCAGCAGGCGTACGAAGCGTCCAATCCTGCGCCGGTGCGCCAGGAGCAGCCAGCGGTCAACTACCGTGACCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTTGTCAGGCTGATTCAt                                                                                                                                                                                                                        >  1:1039060/1‑207 (MQ=255)
                                cAGACTGCCCAGCAGGCGTACGAAGCGTCCAATCCTGCGCCGGTGCGCCAGGAGCAGCCAGCGGTCAACTACCGTGACCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTTGTCAGGCTGATTAATGGAATTTATCAGGGGTATAAAGGTGATTgg                                                                                                                                                                                          <  1:1040718/235‑1 (MQ=255)
                                   aCTGCCCAGCAGGCGTACGAAGCGTCCAATCCTGCGCCGGTGCGCCAGGAGCAGCCAGCGGTCAACTACCGTGACCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTTGTCAGGCTGATTACTGGa                                                                                                                                                                                                                     >  1:112064/1‑205 (MQ=255)
                                     tGCCCAGCAGGCGTACGAAGCGTCCAATCCTGCGCCGGTGCGCCAGGAGCAGCCAGCGGTCAACTACCGTGACCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTTGTCAGGCTGATTAATGGAATTTATCAGGGGTATAAAGGt                                                                                                                                                                                                <  1:490131/224‑1 (MQ=255)
                                     tGCCCAGCAGGCGTACGAAGCGTCCAATCCTGCGCCGGTGCGCCAGGAGCAGCCAGCGGTCAACTACCGTGACCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTTGTCAGGCTGATTAATGGAATTTATCAGGGGTATAAAGGTGATTGGTATtca                                                                                                                                                                                    <  1:992464/236‑1 (MQ=255)
                                      gCCCAGCAGGCGTACGAAGCGTCCAATCCTGCGCCGGTGCGCCAGGAGCAGCCAGCGGTCAACTACCGTGACCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTTGTCAGGCTGATTaa                                                                                                                                                                                                                         >  2:1124880/1‑198 (MQ=255)
                                          agcagGCGTACGCAGCGTCCAATCCTGCGCCGGTGCGCCAGGAGCAGCCAGCGGTCAACTACCGTGACCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTTGTCAGGCTGATTAATGGAATTTATCAGGGGTATAAAGGTGATTgg                                                                                                                                                                                          <  2:543989/225‑1 (MQ=255)
                                          agcagGCGTACGAAGCGTCCAATCCTGCGCCGGTGCGCCAGGAGCAGCCAGCGGTCAACTACCGTGACCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTTGTCAGGCTGATTAATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCaaa                                                                                                                                                                               <  1:643204/236‑1 (MQ=255)
                                          agcagGCGTACGAAGCGTCCAATCCTGCGCCGGTGCGCCAGGAGCAGCCAGCGGTCAACTACCGTGACCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTTGTCAGGCTGATTAATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCaaa                                                                                                                                                                               <  1:533233/236‑1 (MQ=255)
                                                    cGAAGCGTCCAATCCTGCGCCGGTGCGCCAGGAGCAGCCAGCGGTCAACTACCGTGACCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTTGTCAGGCTGATTAATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCt                                                                                                                                                                     >  1:71865/1‑236 (MQ=255)
                                                     gcaGCGTCCAATCCGCCGCCGGTGCGCCAGGAGCAGCCAGCGGTCAACTACCGTGACCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTTGTCAGGCTGATTAATGGAATTTATCAGGGGTATAAAGGTGATTg                                                                                                                                                                                           <  2:739819/211‑1 (MQ=255)
                                                           tCCAATCCTGCGCCGGTGCGCCAGGAGCAGCCAGCGGTCAACTACCGTGACCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTTGTCAGGCTGATTAATGGAATTTAt                                                                                                                                                                                                               <  2:932424/187‑1 (MQ=255)
                                                            ccAATCCTGCGCCGGTGCGCCAGGAGCAGCCAGCGGTCAACTACCGTGACCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTTGTCAGGCTGATTAATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTatca                                                                                                                                                              >  1:515991/1‑235 (MQ=255)
                                                               aTCCTGCGCCGGTGCGCCAGGAGCAGCCAGCGGTCAACTACCGTGACCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTTGTCAGGCTGATTAATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATccc                                                                                                                                                          <  1:215380/236‑1 (MQ=255)
                                                                  cTGCGCCGGTGCGCCAGGAGCAGCCAGCGGTCAACTACCGTGACCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTTGTCAGGCTGATTAATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAg                                                                                                                                                       <  1:469189/236‑1 (MQ=255)
                                                                  cTGCGCCGGTGCGCCAGGAGCAGCCAGCGGTCAACTACCGTGACCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTTGTCAGGCTGATTAATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAg                                                                                                                                                       <  1:1003399/236‑1 (MQ=255)
                                                                       ccGGTGCGCCAGGAGCAGCCAGCGGTCAACTACCGTGACCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTAATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGt                                                                                                                                                                       <  1:84498/215‑1 (MQ=255)
                                                                        cGGTGCGCCAGGAGCAGCCAGCGGTCAACTACCGTGACCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTTGTCAGGCTGATTAATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCttt                                                                                                                                                  >  2:574958/1‑235 (MQ=255)
                                                                             cgcCAGGAGCAGCCAGCGGTCAACTACCGTGACCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTTGTCAGGCTGATTAATGGAAt                                                                                                                                                                                                                   >  2:558809/1‑165 (MQ=255)
                                                                             cgcCAGGAGCAGCCAGCGGTCAACTACCGTGACCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTTGTCAGGCTGATTAATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGAt                                                                                                                                            >  1:191203/1‑236 (MQ=255)
                                                                             cgcCAGGAGCAGCCAGCGGTCAACTACCGTGACCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTTGTCAGGCTGATTAATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGAt                                                                                                                                            >  1:139382/1‑236 (MQ=255)
                                                                                cAGGAGCAGCCAGCGGTCAACTACCGTGACCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTTGTCAGGCTGATTAATGGAATTTAt                                                                                                                                                                                                               <  2:1040097/166‑1 (MQ=255)
                                                                                   gAGCAGCCAGCGGTCAACTACCGTGACCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTTGTCAGGCTGATTAATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAgg                                                                                                                                                                             <  2:198065/197‑1 (MQ=255)
                                                                                         ccagcGGTCAACTACCGTGACCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTTGTCAGGCTGATTAATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGt                                                                                                                                                                   <  2:729058/201‑1 (MQ=255)
                                                                                                      aCCGTGACCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTTGTCAGGCTGATTAATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCAATTTGCGCTGAAGc                                                                                                                   <  1:586108/236‑1 (MQ=255)
                                                                                                       ccGTGACCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTTGTCAGGCTGATTAATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCAATTTGCGCTGAAGc                                                                                                                   <  1:1124811/235‑1 (MQ=255)
                                                                                                        cGTGACCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTTGTCAGGCTGATTAATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCAATTTGCGCTGAAGCg                                                                                                                  <  1:485797/235‑1 (MQ=255)
                                                                                                            aCCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTTGTCAGGCTGATTAATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCAATTTGCGCTGAAGCGGCCaa                                                                                                             <  1:1075129/236‑1 (MQ=255)
                                                                                                                     gcCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTTGTCAGGCTGATTAATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCAATTTGCGCTGAAGCGGCCaa                                                                                                             <  2:810787/227‑1 (MQ=255)
                                                                                                                     gcCCCACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCGGCCGCTAAACTGCTCCCTCAGGATGACCTGGACTACCTCTATCTCCCGGCATTCCTGCGTTGTCAGGCTGAt                                                                                                                                                                                                                            <  2:193135/116‑1 (MQ=255)
                                                                                                                         gACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGGTGACCTGGACTACCCCGATATCCCGGCATTCCTGCGTTGTCAGGCTGATTAATGGAATTTATCAGGGGTTTAAAGGTGATTGGGATTCATCAAAGGTTGGGTCTGTTATCATccc                                                                                                                                                          >  2:792622/1‑178 (MQ=255)
                                                                                                                          aCCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTTGTCAGGCTGATTAATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTg                                                                                                                                                                          <  1:910826/161‑1 (MQ=255)
                                                                                                                           ccGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTTGTCAGGCTGATTAATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCAATTTGCGCTGAAGCGGCCAATGCCATGATTAaaca                                                                                              <  1:907539/236‑1 (MQ=255)
                                                                                                                               cATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTTGTCAGGCTGATTAATGGAATTTATGAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCAATTTGCGCTGAAGCGGCCAATGCCATGATTAAACAACGc                                                                                          <  1:1060914/236‑1 (MQ=255)
                                                                                                                               cATGCGCAACCAGGCGCATGCAGGTGCTGACGCAGCCGCTAAACTGAACCCTCAGGATGACCTTGACTAACTCGATATCCCGGCATGCCTGCGTTTTCAGGCTGATTACTGGAATTTATCAGGGGt                                                                                                                                                                                                        <  2:1011175/126‑1 (MQ=255)
                                                                                                                                aTGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGATCCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTTGTCAGGCTGATTAATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCAATTTGCGCTGAAGCGGCCaa                                                                                                             <  2:993095/216‑1 (MQ=255)
                                                                                                                                 tGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCGCAGGATGTCCTGGACTACCTCGATATCCCGGCATTCCTGCGTTGTCAGGCTGATTAATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCAATTTGCGCTGAAGCg                                                                                                                  <  1:122219/210‑1 (MQ=255)
                                                                                                                                   cgcAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTTGTCAGGCTGATTAATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCAATTTGCGCTGAAGCGGCCAATGCCATGATTAAACAACGCAccc                                                                                      <  1:835083/236‑1 (MQ=255)
                                                                                                                                     cAACCAGGCGCATGCAGGTGCTGCCGCAGCCGTTAGACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCtttgtc                                                                                                                                                                                                                                    <  2:663623/92‑1 (MQ=255)
                                                                                                                                        ccAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTTGTCAGGCTGATTAATGGAGTTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTatcatc                                                                                                                                                            <  2:536561/161‑1 (MQ=255)
                                                                                                                                                   cAGGTGCTGCCGCAGCCGCTAAACTGAGCCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTTGTCAGGCTGATTAATGGAATTTATCAGGGGTATAAAGGTGATTgg                                                                                                                                                                                          <  2:294815/120‑1 (MQ=255)
                                                                                                                                                     ggTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTTGTCAGGCTGATTAATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCAATTTGCGCTGAAGCGGCCAATGCCATGATTAAACAACg                                                                                           >  2:488059/1‑213 (MQ=255)
                                                                                                                                                          tgcCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTTGTCAGGCTGATTAATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTttgttg                                                                                                                                              <  1:28617/157‑1 (MQ=255)
                                                                                                                                                           gcCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTTGTCAGGCTGATTAATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCAATTTGCGCTg                                                                                                                       <  2:665585/179‑1 (MQ=255)
                                                                                                                                                              gCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTTGTCAGGCTGATTAATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCAATTTGCGCTGAAGCGGCCAATGCCATGATTAAACAACGCACCCTGAAGAATACCATCCGTGCCACAGgtg                                                           <  1:973565/236‑1 (MQ=255)
                                                                                                                                                                                   cAGGATGACCTGGAATACCTCGCTATCCCGGCATTCCTGCGTTGTCAGGCTGATTAATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCAATTTGCGCTGAAGCGGCCaa                                                                                                             <  1:486782/165‑1 (MQ=255)
                                                                                                                                                                                    aGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTTGTCAGGCTGATTAATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCAATTTGCGCTGAAGCGGCCAATGCCATGATTAAACAACGCACCCTGAAGAATACCATCCGTGCCACAGGTGTCGGTCTGCACTCCGGGgaga                                      <  1:688174/235‑1 (MQ=255)
                                                                                                                                                                                        tGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTTGTCAGGCTGATTAATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCAATTTGCGCTGAAGCGGCCAATGCCATGATTAaacaa                                                                                             <  2:944213/176‑1 (MQ=255)
                                                                                                                                                                                         gACCTGGACTACCTCGATATCCCGGCATTCCTGCGTTGTCAGGCTGATTAATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCAATTTGCGCTGAAGCGGCCAATGCCATGATTAAACAACGCACCCTGAAGAATACCATCCGTGCCACAGGTGTCGGTCTGCACTCCGGGGAGAAGGttt                                <  1:285480/236‑1 (MQ=255)
                                                                                                                                                                                            cTGGACTACCTCGATATCCCGGCATTCCTGCGTTGTCAGGCTGATTAATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCAATTTGCGCTGAAGCGGCCAATGCCATGATTAAACAACGCACCCTGAAGAATACCATCCGTGCCACAGGTGTCGGTCTGCACTCCGGGGAGAAGGTTTAc                              <  1:411430/235‑1 (MQ=255)
                                                                                                                                                                                             tGGACTACCTCGATATCCCGGCATTCCTGCGTTGTCAGGCTGATTAATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGc                                                                                                                                  <  1:989634/134‑1 (MQ=255)
                                                                                                                                                                                             tGGACTACCTCGATATCCCGGCATTCCTGCGTTGTCAGGCTGATTAATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGc                                                                                                                                  >  2:989633/1‑134 (MQ=255)
                                                                                                                                                                                             tGGACTACCTCGATATCCCGGCATTCCTGCGTTGTCAGGCTGATTAATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCAATTTGCGCTGAAGCGGCCAATGCCATGATTAAACAACGCACCCTGAAGAATAc                                                                            >  1:608720/1‑188 (MQ=255)
                                                                                                                                                                                             tGGACTACCTCGATATCCCGGCATTCCTGCGTTGTCAGGCTGATTAATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGCTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTCCGCAATTTGCGCTGAAGCGGCCAATGCCATGATTAAACAACGCACCCTGAAGAATAc                                                                            <  2:608719/188‑1 (MQ=255)
                                                                                                                                                                                                aCTACCTCGATATCCCGGCATTCCTGCGTTGTCAGGCTGATTAATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCAATTTGCGCTGAAGCGGCCAATGCCATGATTAAACAACGCACCCTGAAGAATACCATCCGTGCCACAGGTGTCGGTCTGCACTCCGGGGAGAAGGTTTACCTGAc                         <  1:614640/236‑1 (MQ=255)
                                                                                                                                                                                                 cTACCTCGATATCCCGGCATTCCTGCGTTGTCAGGCTGATTAATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCAATTTGCGCTGAAGCGGCCAATGCCATGATTAAACAACGCACCCTGAAGAATACCATCCGTGCCACAGGTGTCGGTCTGCACTCCGGGGAGAAGGTTTACCTGAcc                        >  1:676998/1‑236 (MQ=255)
                                                                                                                                                                                                       cGATATCCCGGCATTCCTGCGTTGTCAGGCTGATTAATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCAATTTGCGCTGAAGCGGCCAATGCCATGATTAAACAACGCACCCTGAAGAATACCATCCGTGCCACAGGTGTCGGTCTGCACTCCGGGGAGAAGGTTTACCTGACCCTCAAg                  >  1:1056337/1‑236 (MQ=255)
                                                                                                                                                                                                       cGATATCCCGGCATTCCTGCGTTGTCAGGCTGATTAATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCAATTTGCGCTGAAGCGGCCAATGCCATGATTAAACAACGCACCCTGAAGAATACCATCCGTGCCACAGGTGTCGGTCTGCACTCCGGGCAGACGGt                                  >  1:410098/1‑220 (MQ=255)
                                                                                                                                                                                                        gATATCCCGGCATTCCTGCGTTGTCAGGCTGATTAATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCAATTTGCGCTGAAGCGGCCAATGCCATGATTAAACAACGCACCCTGAAGAATCCCATCCGTGCCACCGGTGTCGGTCTGCACTCCGGGga                                        >  2:582979/1‑213 (MQ=255)
                                                                                                                                                                                                           atCCCGTCATTCCTGCGTTGTCAGGCTGATTAATGGAATTTATCAGGGGTATAAAGTTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCAATTTGCGCTg                                                                                                                       <  1:969430/131‑1 (MQ=255)
                                                                                                                                                                                                           atCCCGGCATTCCTGCGTTGTCAGGCTGATTAATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCAATTTGCGCTGAAGCGGCCAATGCCATGATTAAACAACGCACCCTGAAGAATACCATCCGTGCCACAGGTGTCGGTCTGCACTCCGGGGAGAAGGTTTACCTGACCCTCAAGc                 <  2:868169/233‑1 (MQ=255)
                                                                                                                                                                                                            tCCCGGCATTCCTGCGTTGTCAGGCTGATTAATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCAATTTGCGCTGAAGCGGCCAATGCCATGATTAAACAACGCACCCTGAAGAATACCATACGTACcacacg                                                             >  1:1059713/1‑186 (MQ=255)
                                                                                                                                                                                                               cGGCATTCCTGCGTTGTCAGGCTGATTAATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTgg                                                                                                                                                                         <  1:752666/77‑1 (MQ=255)
                                                                                                                                                                                                                ggCATTCCTGCGTTGTCAGGCTGATTAATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCAATTTGCGCTGAAGCGGCCAATGCCATGATTAaac                                                                                               <  2:532867/150‑1 (MQ=255)
                                                                                                                                                                                                                 gCATTCCTGCGTTGTCAGGCTGATTAATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCAATTTGCGCTGAAGCGGCCaa                                                                                                             <  2:987411/135‑1 (MQ=255)
                                                                                                                                                                                                                 gCATTCCTGCGTTGTCAGGCTGATTAATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCAATTTGCGCTGAAGCGGCCAATGCCATGATTAAACAACGCACCCTGAAGAATACCATCCGTGCCACAGGTGTCGGTCTGCACTCCGGGGAGACGGt                                  >  1:230573/1‑210 (MQ=255)
                                                                                                                                                                                                                  cATTCCTGCGTTGTCAGGCTGATTAATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCAATTTGCGCTGAAGCGGCCAATGCCATGATTAAACAACGCACCCTGAAGAATACCATCCGTGCCACAGGTGTCGGTCTGCACTCCGGGGAGAAGGTTTAc                              >  2:358409/1‑213 (MQ=255)
                                                                                                                                                                                                                       cTGCGTTGTCAGGCTGATTAATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCAATTTGCGCTGAAGCGGCCAATGCCATGATTAAACAACGCACCCTGAAGAATACCATCCGTGCCACAGGTGTCGGTCTGCACTCCGGGGAGAAGGTTTACCTGACCCTCAAGCCTGCACCTGTTGaca  >  1:265228/1‑236 (MQ=255)
                                                                                                                                                                                                                             |                                                                                                                                                                                                                                     
AGAAGCCTGTGAAGGTGGTCGACAATACCCTGCAGACTGCCCAGCAGGCGTACGAAGCGTCCAATCCTGCGCCGGTGCGCCAGGAGCAGCCAGCGGTCAACTACCGTGACCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTAATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCAATTTGCGCTGAAGCGGCCAATGCCATGATTAAACAACGCACCCTGAAGAATACCATCCGTGCCACAGGTGTCGGTCTGCACTCCGGGGAGAAGGTTTACCTGACCCTCAAGCCTGCACCTGTTGACA  >  AE015451/1529939‑1530389

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