Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 6,151,153 T→C 5.6% V55A (GTC→GCC)  PP_5395 → conserved protein of unknown function

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE0154516,151,1530TC 5.6% 236.4 / 3.6 71V55A (GTC→GCC) PP_5395conserved protein of unknown function
Reads supporting (aligned to +/- strand):  ref base T (41/25);  new base C (2/2);  total (44/27)
Fisher's exact test for biased strand distribution p-value = 6.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.71e-01

GCAGGTAGCCCTTAGCCATTAGCTGTCCATCAACGCTAAGGATCAGCCGCTGACATGGAGGTTGATAGTGAAGTTCGATCTAAGGAATGCGATCAATAGCATATCTGACGTGGTGATGAATCGACCCGCCCCCTTGCGCGAGTTCGACCAGATCTACATGAAGGTCGCGGATATGGTCATCCAGGCCGAGTACGTCGCAAGGGTCTTCGACGGGAAGAAGGTGGTCTTCGTTGGTGACGGCGATGGGATAGCCCTTTCAGCTGCGCATCTCAAAGCCCAGGAGGTAATCAGCTACGGGCCAAAATCGATCACGCTCCTGGATTTCGATGAGCGAATCGTCAACTCAGTACGGCGATTCGGTGAAAAACTAGCGCCCAAGATCGAAATCACTGCCCATCTCTACAACGTGGCTGACCCTCTACCTGCCGAACATTTTGGGGCCTACGACGGCTTTCATATC  >  AE015451/6150929‑6151388
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                           
gCAGGTAGCCCTTAGCCATTAGCTGTCCATCAACGCTAAGGATCAGCCGCTGACATGGAGGTTGATAGTGAAGTTCGATCTAAGGAATGCGATCAATAGCATATCTGACGTGGTGATGAATCGACCCGCCCCCTTGCGCGAGTTCGACCAGATCTACATGAAGGTCGCGGATATGGTCATCCAGGCCGAGTACGTCGCAAGGGTCTTCGACGGGAAGAAGGTGGTCTTCGTTGGTg                                                                                                                                                                                                                                  >  1:824309/1‑236 (MQ=255)
          cTTAGCCATTAGCTGTCCATCAACGCGAAGGATCAGCCGCTGACATGGAGGTTGATAGTGAAGTTCGATCTAAGGAATGCGATCAATAGCATATCTGACGTGGTGATGAATCGACCCGCCCCCTTGCGCGAGTTCGACCAGATCTACATGAAGGTCGCGGATATGGTCATCCAGGCCGAGTACGTCGCAAGGGTCTTCGACGGGAAGAAGGTGGTCTTCGTTGGTGACGGCGAt                                                                                                                                                                                                                          <  1:34875/234‑1 (MQ=255)
           ttAGCCATTAGCTGTCCATCAACGCTAAGGATCAGCCGCTGACATGGAGGTTGATAGTGAAGTTCGATCTAAGGAATGCGATCAATAGCATATCTGACGTGGTGATGAATCGACCCGCCCCCTTGCGCGAGTTCGACCAGATCTACATGAAGGTCGCGGATATGGTCATCCAGGCCGAGTACGTCGCAAGGGTCTTCGACGGGAGGAAGGTGGTCTTCGTTGGTGACGGCGATgg                                                                                                                                                                                                                        >  1:638764/1‑235 (MQ=255)
                      cTGTCCATCAACGCTAAGGATCAGCCGCTGACATGGAGGTTGATAGTGAAGTTCGATCTAAGGAATGCGATCAATAGCATATCTGACGTGGTGATGAATCGACCCGCCCCCTTGCGCGAGTTCGACCAGATCTACATGAAGGTCGCGGATATGGTCATCCAGGCCGAGTACGTCGCAAGGGTCTTCGACGGGAAGAAGGTGGTCTTCGTTGGTGACGGCGATGGGATAGCCCttta                                                                                                                                                                                                            >  1:361597/1‑235 (MQ=255)
                      cTGTCCATCAACGCTAAGGATCAGCCGCTGACATGGAGGTTGATAGTGAAGTTCGATCTAAGGAATGCGATCAATAGCATATCTGACGTGGTGATGAATCGACCCGCCCCCTTGCGCGAGTTCGACCAGATCTACATGAAGGTCGCGGATATGGTCATCCAGGCCGAGTACGTCGCAAGGGTCTTCGACGGGAAGAAGGTGGCCTTCGTTGGTGACGGCGATGGGATAGCCCTTTc                                                                                                                                                                                                            >  1:130513/1‑236 (MQ=255)
                      cTGTCCATCAACGCTAAGGATCAGCCGCTGACATGGAGGTTGATAGTGAAGTTCGATCTAAGGAATGCGATCAATAGCATATCTGACGTGGTGATGAATCGACCCGCCCCCTTGCGCGAGTTCGACCAGATCTACATGAAGGTCGCGGATATGGTCATCCAGGCCGAGTACGTCGCAAGGGTCTTCGACGGGAAGAAGGTGGCCTTCGTTGGTGACGGCGATGGGATAGCCCTTTc                                                                                                                                                                                                            >  1:130504/1‑236 (MQ=255)
                       tGTCCATCAACGCTAAGGATCAGCCGCTGACATGGAGGTTGATAGTGAAGTTCGATCTAAGGAATGCGATCAATAGCATATCTGACGTGGTGATGAATCGACCCGCCCCCTTGCGCGAGTTCGACCAGATCTACATGAAGGTCGCGGATATGGTCATCCAGGCCGAGTACGTCGCAAGGGTCTTCGACGGGAAGAAGGTGGTCTTCGTTGGTGACGGCGATGGGATAGCCCTTTCa                                                                                                                                                                                                           >  1:920110/1‑236 (MQ=255)
                          ccATCAACGCTAAGGATCAGCCGCTGACATGGAGGTTGATAGTGAAGTTCGATCTAAGGAATGCGATCAATAGCATATCTGACGTGGTGATGAATCGACCCGCCCCCTTGCGCGAGTTCGACCAGATCTACATGAAGGTCGCGGATATGGTCATCCAGGCCGAGTACGTCGCAAGGGTCTTCGACGGGAAGAAGGTGGTCTTCGTTGGTGACGGCGATggg                                                                                                                                                                                                                       <  1:293371/221‑1 (MQ=255)
                            aTCAACGCTAAGGATCAGCCGCTGACATGGAGGTTGATAGTGAAGTTCGATCTAAGGAATGCGATCAATAGCATATCTGACGTGGTGATGAATCGACCCGCCCCCTTGCGCGAGTTCGACCAGATCTACATGAAGGTCGCGGATATGGTCATCCAGGCCGAGTACGTCGCAAGGGTCTTCGACGGGAAGAAGGTGGTCTTCGTTGGTGACGGCGATGGGATAGCCCtg                                                                                                                                                                                                              >  2:655129/1‑227 (MQ=255)
                            aTCAACGCTAAGGATCAGCCGCTGACATGGAGGTTGATAGTGAAGTTCGATCTAAGGAATGCGATCAATAGCATATCTGACGTGGTGATGAATCGACCCGCCCCCTTGCGCGAGTTCGACCAGATCTACATGAAGGTCGCGGATATGGTCATCCAGGCCGAGTACGTCGCAAGGGTCTTCGACGGGAAGAAGGTGGTCTTCGTTGGTGACGGCGATGGGATAGCCCTTTCAGCt                                                                                                                                                                                                        >  2:381355/1‑234 (MQ=255)
                             tCAACGCTAAGGATCAGCCGCTGACATGGAGGTTGATAGTGAAGTTCGATCTAAGGAATGCGATCAATAGCATATCTGACGTGGTGATGAATCGACCCGCCCCCTTGCGCGAGTTCGACCAGATCTACATGAAGGTCGCGGATATGGTCATCCAGGCCGAGTACGTCGCAAGGGTCTTCGACGGGAGGAAGGTGGTCTTCGTTGGTGACGGCGATGGGATAGCCCTTTCAGCt                                                                                                                                                                                                        <  2:638761/233‑1 (MQ=255)
                                   cTAAGGATCAGCCGCTGACATGGAGGTTGATAGTGAAGTTCGATCTAAGGAATGCGATCAATAGCATATCTGACGTGGTGATGAATCGACCCGCCCCCTTGCGCGAGTTCGACCAGATCTACATGAAGGTCGCGGATATGGTCATCCAGGCCGAGTACGTCGCAAGGGTCTTCGACGGGAAGAAGGTGGTCTTCGTTGGTGACGGCGATGGGATAGCCCTTTCAGCt                                                                                                                                                                                                        <  2:326657/227‑1 (MQ=255)
                                      aGGATCAGCCGCTGACATGGAGGTTGATAGTGAAGTTCGATCTAAGGAATGCGATCAATAGCATATCTGACGTGGTGATGAATCGACCCGCCCCCTTGCGCGAGTTCGACCAGATCTACATGAAGGTCGCGGATATGGTCATCCAGGCCGAGTACGTCGCAAGGGTCTTCGACGGGAAGAAGGTGGTCTTCGTTGGTGACGGCGATGGGATAGCCCTTTCAGCTGCGCATCTCaa                                                                                                                                                                                             >  1:917228/1‑235 (MQ=255)
                                            aGCCGCTGACATGGAGGTTGATAGTGAAGTTCGATCTAAGGAATGCGATCAATAGCATATCTGACGTGGTGATGAATCGACCCGCCCCCTTGCGCGAGTTCGACCAGATCTACATGAAGGTCGCGGATATGGTCATCCAGGCCGAGTACGTCGCAAGGGTCTTCGACGGGAAGAAGGTGGTCTTCGTTGGTGACGGCGATGGGATAGCCCTTTCAGCTGCGCATCTCAAAGCCCag                                                                                                                                                                                      <  1:330357/236‑1 (MQ=255)
                                             gCCGCTGACATGGAGGTTGATAGTGAAGTTCGATCTAAGGAATGCGATCAATAGCATATCTGACGTGGTGATGAATCGACCCGCCCCCTTGCGCGAGTTCGACCAGATCTACATGAAGGTCGCGGATATGGTCATCCAGTCCGAGTACGTCGCAAGGGTCTTCGAACGGGAAGAAGGTGGTCTTCGTTGGTGACGGCGATGGGATAGCCCTTTCAGCTGCGCATCTCAAAGCCCag                                                                                                                                                                                      <  1:391370/236‑1 (MQ=255)
                                               cGCTGACATGGAGGTTGATAGTGAAGTTCGATCTAAGGAATGCGATCAATAGCATATCTGACGTGGTGATGAATCGACCCGCCCCCTTGCGCGAGTTCGACCAGATCTACATGAAGGTCGCGGATATGGTCATCCAGGCCGAGTACGTCGCAAGGGTCTTCGACGGGAAGAAGGTGGTCTTCGTTGGTGACGGCGATGGGATAGCCCTTTCAGCTGCGCATCTCAAAGCCCaggag                                                                                                                                                                                   >  1:562072/1‑236 (MQ=255)
                                                     cATGGAGGTTGATAGTGAAGTTCGATCTAAGGAATGCGATCAATAGCATATCTGACGTGGTGATGAATCGACCCGCCCCCTTGCGCGAGTTCGACCAGATCTACATGAAGGTCGCGGATATGGTCATCCAGGCCGAGTACGTCGCAAGGGTCTTCGACGGGAAGAAGGTGGTCTTCGTTGGTGACGGCGATGGGATAGCCCTTTCAGCTGCGCa                                                                                                                                                                                                   >  2:363311/1‑214 (MQ=255)
                                                            gTTGATAGTGAAGTTCGATCTAAGGAATGCGATCAATAGCATATCTGACGTGGTGATGAATCGACCCGCCCCCTTGCGCGAGTTCGACCAGATCTACATGAAGGTCGCGGATATGGTCATCCAGGCCGAGTACGTCGCAAGGGTCTTCGACGGGAAGAAGGTGGGCTTCGTTGGTGACGGCGATGGGATAGCCCTTTCAGCTGCGCATCTCAAAGCCCAGGAGGTAATCAGCTACg                                                                                                                                                                      >  1:167325/1‑236 (MQ=255)
                                                                   gTGAAGTTCGATCTAAGGAATGCGATCAATAGCATATCTGACGTGGTGATGAATCGACCCGCCCCCTTGCGCGAGTTCGACCAGATCTACATGAAGGTCGCGGATATGGTCATCCAGGCCGAGTACGTCGCAAGGGTCTTCGACGGGAAGAAGGTGGTCTTCGTTGGTGACGGCGATGGGATAGCACTTTCAGCTGCGCATCTCAAAGCCCAGGAGGTAATCAGCTACGGGcacaa                                                                                                                                                               >  1:80193/1‑232 (MQ=255)
                                                                              tCTAAGGAATGCGATCAATAGCATATCTGACGTGGTGATGAATCGACCCGCCCCCTTGCGCGAGTTCGACCAGATCTACATGAAGGTCGCGGATATGGTCATCCAGGCCGAGTACGTCGCAAGGGTCTTCGACGGGAAGAAGGTGGTCTTCGTTGGTGACGGCGATGGGATAGCCCTTTCAGCTGCGCATCTCAAAGCCGAGGAGGTAATCAGCTACGGGCCaaa                                                                                                                                                               >  1:1065575/1‑225 (MQ=255)
                                                                              tCTAAGGAATGCGATCAATAGCATATCTGACGTGGTGATGAATCGACCCGCCCCCTTGCGCGAGTTCGACCAGATCTACATGAAGGTCGCGGATATGGTCATCCAGGCCGAGTACGTCGCAAGGGTCTTCGACGGGAAGAAGGTGGTCTTCGTTGGTGACGGCGATGGGATAGCCCTTTCAGCTGCGCATCTCAAAGCCCAGGAGGTAATCAGCTACGGGCCAAAATCg                                                                                                                                                           >  2:742136/1‑229 (MQ=255)
                                                                                 aaGGAATGCGATCAATAGCATATCTGACGTGGTGATGAATCGACCCGCCCCCTTGCGCGAGTTCGACCAGATCTACATGAAGGTCGCGGATATGGTCATCCAGGCCGAGTACGTCGCAAGGGTCTTCGACGGGAAGAAGGTGGTCTTCGTTGGTGACGGCGATGGGATAGCCCTTTCAGCTGCGCATCTCAAAGCCCAGGAGGTAATCAGCTACGGGCCAAAATCGATCACg                                                                                                                                                     >  2:818656/1‑232 (MQ=255)
                                                                                       tGCGATCAATAGCATATCTGACGTGGTGATGAATCGACCCGCCCCCTTGCGCGAGTTCGACCAGATCTACATGAAGGTCGCGGATATGGTCATCCAGGCCGAGTACGTCGCAAGGGTCTTCGACGGGAAGAAGGTGGCCTTCGTTGGTGACGGCGATGGGATAGCCCTTTCAGCTGCGCATCTCAAAGCCCAGGAGGTAATCAGCTACGGGCCAAAATCGATc                                                                                                                                                        <  2:130502/223‑1 (MQ=255)
                                                                                          gATCAATAGCATATCTGACGTGGTGATGAATCGACCCGCCCCCTTGCGCGAGTTCGACCAGATCTACATGAAGGTCGCGGATATGGTCATCCAGGCCGAGTACGTCGCAAGGGTCTTCGACGGGAAGAAGGTGGTCTTCGTTGGTGACGGCgc                                                                                                                                                                                                                           >  1:1031278/1‑152 (MQ=255)
                                                                                          gATCAATAGCATATCTGACGTGGTGATGAATCGACCCGCCCCCTTGCGCGAGTTCGACCAGATCTACATGAAGGTCGCGGATATGGTCATCCAGGCCGAGTACGTCGCAAGGGTCTTCGACGGGAAGAAGGTGGTCTTCGTTGGTGACGGCGATGGGATAGCCCTTTCAGCTGCGCATCTCAAAGCCCAGGAGGTAATCAGCTACggg                                                                                                                                                                    >  1:185457/1‑208 (MQ=255)
                                                                                          gATCAATAGCATATCTGACGTGGTGATGAATCGACCCGCCCCCTTGCGCGAGTTCGACCAGATCTACATGAAGGTCGCGGATATGGTCATCCAGGCCGAGTACGTCGCAAGGGTCTTCGACGGGAAGAAGGTGGTCTTCGTTGGTGACGGCGATGGGATAGCCCTTTCAGCTGCGCATCTCAAAGCCCAGGAGGTAATCAGCTACGGGCCAAACTCGATCACGCTCCTGGATTTc                                                                                                                                         >  1:461241/1‑235 (MQ=255)
                                                                                          gATCAATAGCATATCTGACGTGGTGATGAATCGACCCGCCCCCTTGCGCGAGTTCGACCAGATCTACATGAAGGTCGCGGATATGGTCATCCAGGCCGAGTACGTCGCAAGGGTCTTCGACGGGAAGAAGGTGGTCTTCGTTGGTGACGGCGATGGGATAGCCCTTTCAGCTGCGCATCTCAAAGCCCAGGAGGTAATCAGCTACGGGCCAAAATCGATCACGCTCCTGGATTTCg                                                                                                                                        >  1:856693/1‑236 (MQ=255)
                                                                                              aaTAGCATATCTGACGTGGTGATGAATCGACCCGCCCCCTTGCGCGAGTTCGACCAGATCTACATGAAGGTCGCGGATATGGTCATCCAGGCCGAGTACGTCGCAAGGGTCTTCGACGGGAAGAAGGTGGTCTTCGTTGGTGACGGCGATGGGATAGCCCTTTCAGCTGCGCATCTCAAAGCCCAGGAGGTAATCAGCTACGGGCCAAAATCGATCACGCTCCTGGATTTCGATg                                                                                                                                     <  1:16213/235‑1 (MQ=255)
                                                                                               aTAGCATATCTGACGTGGTGATGAATCGACCCGCCCCCTTGCGCGAGTTCGACCAGATCTACATGAAGGTCGCGGATATGGTCATCCAGGCCGAGTACGTCGCAAGGGTCTTCGACGGGAAGAAGGTGGTCTTCGTTGGTGACGGCGATGGGATAGCCCTTTCAGCTGCGCATCTCAAAGCCCAGGCGGTa                                                                                                                                                                                >  2:311548/1‑191 (MQ=255)
                                                                                               aTAGCATATCTGACGTGGTGATGAATCGACCCGCCCCCTTGCGCGAGTTCGACCAGATCTACATGAAGGTCGCGGATATGGTCATCCAGGCCGAGTACGTCGCAAGGGTCTTCGACGGGAAGAAGGTGGTCTTCGTTGGTGACGGCGATGGGATAGCCCTTTCAGCTGCGCATCTCAAAGCCCAGGAGGTAATCAGCTACGGGCCAAAATCGATCACGCTCCTGGATTTCGATGa                                                                                                                                    >  1:856090/1‑235 (MQ=255)
                                                                                                tAGCATATCTGACGTGGTGATGAATCGACCCGCCCCCTTGCGCGAGTTCGACCAGATCTACATGAAGGTCGCGGATATGGTCATCCAGGCCGAGTACGTCGCAAGGGTCTTCGACGGGAAGAAGGTGGTCTTCGTTg                                                                                                                                                                                                                                     <  1:951085/137‑1 (MQ=255)
                                                                                                           aCGTGGTGATGAATCGACCCGCCCCCTTGCGCGAGTTCGACCAGATCTACATGAAGGTCGCGGATATGGTCATCCAGGCCGAGTACGTCGCAAGGGTCTTCGACGGGAAGAAGGTGGTCTTCGTTGGTGACGGCGATGGGATAGCCCTTTCAGCTGCGCATCTCAAGGCCCag                                                                                                                                                                                      >  2:503125/1‑173 (MQ=255)
                                                                                                           aCGTGGTGATGAATCGACCCGCCCCCTTGCGCGAGTTCGACCAGATCTACATGAAGGTCGCGGATATGGTCATCCAGGCCGAGTACGTCGCAAGGGTCTTCGACGGGAAGAAGGTGGTCTTCGTTGGTGACGGCGATGGGATAGCCCTTTCAGCTGCGCATCTCAAAGCCCAGGAGGTAATCAGCTACGGGCCAAAATCGATCACGCTCCTGGATTTCGATGAGCGAATCGTCAAc                                                                                                                       >  1:624668/1‑236 (MQ=255)
                                                                                                                         cGACCCGCCCCCTTGCGCGAGTTTGACCATATCTACATGAAGGTCGCGGATATGGTCATCCAGGCCGAGTACGTCGCAAGGGTCTTCGACGTGAAGAAGGTGGTCTTCGTTGGTGACGGCGATGGGATAGCCCTTTCAGCTGCGCATCTCAAAg                                                                                                                                                                                           <  2:51394/154‑1 (MQ=255)
                                                                                                                         cGACCCGCCCCCTTGCGCGAGTTCGACCAGATCTACATGAAGGTCGCGGATATGGTCATCCAGGCCGAGTACGTCGCAAGGGTCTTCGACGGGAAGAAGGTGGTCTTCGTTGGTGACGGCGATGGGATAGCCCTTTCAGCTGCGCATCTCAAAGCCCAGGAGGTAATCAGCTACGGGCCAAAATCGATCACGCTCCTGGATTTCGATg                                                                                                                                     <  1:926045/208‑1 (MQ=255)
                                                                                                                              cGCCCCCTTGCGCGAGTTCGACCAGATCTACATGAAGGTCGCGGATATGGTCATCCAGGCCGAGTACGTCGCAAGGGTCTTCGACGGGAAGAAGGTGGTCTTCGTTGGTGACGGCGCTGGGATAGCCCTTTCAGCTGCGCATCTCAAAGCCCAGGCGGTAATCAGCTACGGGcc                                                                                                                                                                  >  2:18367/1‑174 (MQ=255)
                                                                                                                                   ccTTGCGCGAGTTCGACCAGATCTACATGAAGGTCGCGGATATGGTCATCCAGGCCGAGTACGTCGCAAGGGTCTTCGACGGGAAGAAGGTGGTCTTCGTTGGTGACGGCGATGGGATAGCCCTTTCAGCTGCGCATCTCAAAGCCCAGGAGGTAATCAGCTACGGGCCAAAATCGATCACGCTCCTGGATTTCGATg                                                                                                                                     >  2:407229/1‑198 (MQ=255)
                                                                                                                                          ccaGTTCGACCAGATCTACATGAAGGTCGCGGATATGGTCATCCAGGCCGAGTACGTCGCAAGGGTCTTCGCCGGGAAGAAGGTGGTCTTCGTTGGTGACg                                                                                                                                                                                                                               <  2:657838/99‑1 (MQ=255)
                                                                                                                                                cGACCAGATCTACATGAAGGTCGCGGATATGGTCATCCAGGCCGAGTACGTCGCAAGGGTCTTCGACGGGAAGAAGGTGGTCTTCGTTGGTGACGGCGATGGGATAGCCCTTTCAGCTg                                                                                                                                                                                                       >  2:399650/1‑119 (MQ=255)
                                                                                                                                                    cAGATCTACATGAAGGTCGCGGATATGGTCATCCAGGCCGAGTACGTCGCAAGGGTCTTCGACGGGAAGAAGGTGGTCTTCGTTGGTGACGGCGATGGGATAGCCCTTTCAGCTGCGCATCTCAAAGCCCAGGAGGTAATCAGCTACGGGCCAAAATCGATCACGCTCCTGGATTTCGATGAGCGAATCGTCAACTCAGTACGGCGATTCGGTGAAAAACTAGCGCCCAAGATCGa                                                                              <  1:35949/236‑1 (MQ=255)
                                                                                                                                                        tCTACATGAAGGTCGCGGATATGGTCATCCAGGCCGAGTACGTCGCAAGGGTCTTCGACGGGAAGAAGGTGGTCTTCGTTGGTGACGGCGATGGGATAGCCCTTTCAGCTGCGCATCTCAAAGCCCAGGAGGTAATCAGCTACGGGCCAAAATCGATCACGCTCCTGGATTTCGATGAGCGAATCGTCAACTCAGTACGGCGATTCGGTGAAAAACTAGCGCCCAAGATCGaaa                                                                            <  1:772013/234‑1 (MQ=255)
                                                                                                                                                         cTACATGAAGGTCGCGGATATGGTCATCCAGGCCGAGTACGTCGCAAGGGTCTTCGACGGGAAGAAGGTGGTCTTCGTTGGTGACGGCGATGGGATAGCCCTTTCAGCTGCGCATCTCAAAGCCCAGGAGGTAATCAGCTACGGGCCAAAATCGATCACGCTCCTGGATTTCGATGAGCGAATCGTCAACTCAGTACGGCGATTCGGTGAAAAACTAGCGCCCAAGATCGAAATCa                                                                         <  1:283748/236‑1 (MQ=255)
                                                                                                                                                          tACATGAAGGTCGCGGATATGGTCATCCAGGCCGAGTACGTCGCAAGGGTCTTCGACGGGACGAAGGTGGTCTTCGTTGGTGACGGCGATGGGATAGCCCTTTCAGCTGCGCATCTCAAAGCCCAGGAGGTAATCAGCTACGGGCCAAAGTCGATCACgacc                                                                                                                                                  >  2:294677/1‑159 (MQ=255)
                                                                                                                                                           aCATGAAGGTCGCGGATATGGTCATCCAGGCCGAGTACGTCGCAAGGGTCTTCGACGGGAAGAAGGTGGTCTTCGTTGGTGACGGCGATGGGATAGCCCTTTCAGCTGCGCATCTCAAAGCCCAGGAGGTAATCAGCTACGGGCCAAAATCGATCACGCTCCTGGATTTCGATGAGCGAATCGTCAACTCAGTACGGCGATTCGGTGAAAAACTAGCGCCCAAGATCGAAATCACt                                                                       <  1:440958/236‑1 (MQ=255)
                                                                                                                                                            cATGAAGGTCGCGGATATGGTCATCCAGGCCGAGTACGTCGCAAGGGTCTTCGACGGGAAGAAGGTGGTCTTCGTTGGTGACGGCGATGGGATAGCCCTTTCAGCTGCGCATCTCAAAGCCCAGGAGGTAATCAGCTACGGGCCAAAATCGATCACGCTCCTGGATTTCGATGAGCGAATCGTCAACTCAGTACGGCGATTCGGTGAAAAACTAGCGCCCAAGATCGAAATCACt                                                                       <  1:623933/235‑1 (MQ=255)
                                                                                                                                                             acgAAGGTCGCGGATATGGTCATCCAGGCCGAGTACGTCTCAAGGGTCTTCGACGGGAAGAAGGTGGTCTTCGTTGGTGACGGCGATGGGATAGCCCTTTCAGCTGCGCATCTCAAAGCCCAGGAGGTAATCAGCTACGCGCCAAAATCGATCACGCTCCTGGATTTCGATg                                                                                                                                     <  2:728387/170‑1 (MQ=255)
                                                                                                                                                               gAAGGTCGCGGATATGGTCATCCAGGCCGAGTACGTCGCAAGGGTCTTCGACGGGAAGAAGGTGGTCTTCGTTGGTGACGGCGATGGGATAGCCCTTTCAGCTGCGCATCTCAAAGCCCAGGAGGTAATCAGCTACGGGCCAAAATCGATCACGCTCCTGGATTTCGATGAGCGAATCGTCAACTCAGTACGGCGATTCGGTGAAAAACTAgcgc                                                                                        <  2:407534/215‑1 (MQ=255)
                                                                                                                                                               gAAGGTCGCGGATATGGTCATCCAGGCCGAGTACGTCGCAAGGGTCTTCGACGGGAAGAAGGTGGTCTTCGTTGGTGACGGCGATGGGATAGCCCTTTCAGCTGCGCATCTCAAAGCCCAGGAGGTAATCAGCTACGGGCCAAAATCGATCACGCTCCTGGATTTCGATGAGCGAATCGTCAACTCAGTACGGCGATTCGGTGAAAAACTAGCGCCCAAGATCGAAATCACt                                                                       <  2:920107/232‑1 (MQ=255)
                                                                                                                                                                  ggTCGCGGATATGGTCATCCAGGCCGAGTACGTCGCAAGGGTCTTCGACGGGAAGAAGGTGGTCTTCGTTGGTGACGGCGATGGGATAGCCCTTTCAGCTGCGCATCTCAAAGCCCAGGAGGTAATCAGCTACGGGCCAAAATCGATCACGCTCCTGGATTTCGATGAGCGAATCGTCAACTCAGTACGGCGATTCGGTGAAAAACTAGCGCCCAAGATCGAAATCACTGcc                                                                    <  2:167323/232‑1 (MQ=255)
                                                                                                                                                                    tCGCGGATATGGTCATCCAGGCCGAGTACGTCGCAAGGGTCTTCGACGGGAAGAAGGTGGTCTTCGTTGGTGACGGCGATGGGATAGCCCTTTCAGCTGCGCATCTCAAAGCCCAGGAGGTAATCAGCTACGGGCCAAAATCGATCACGCTCCTGGCTTTCGATGCGCGAATCGTcca                                                                                                                        >  2:61921/1‑176 (MQ=255)
                                                                                                                                                                    tCGCGGATATGGTCATCCAGGCCGAGTACGTCGCAAGGGTCTTCGACGGGAAGAAGGTGGTCTTCGTTGGTGACGGCGATGGGATAGCCCTTTCAGCTGCGCATCTCAAAGCCCAGGAGGTAATCAGCTACGGGCCAAAATCGATCACGCTCCTGGATTTCGATGAGCGAATc                                                                                                                             >  2:1074497/1‑173 (MQ=255)
                                                                                                                                                                    tCGCGGATATGGTCATCCAGGCCGAGTACGTCGCAAGGGTCTTCGACGGGAAGAAGGTGGTCTTCGTTGGTGACGGCGATGGGATAGCCCTTTCAGCTGCGCATCTCAAAGCCCAGGAGGTAATCAGCTACGGGCCAAAATCGATCACGCTCCTGGATTTCGATGAGCGAATCGTCAACTCAGTAc                                                                                                                >  2:699155/1‑186 (MQ=255)
                                                                                                                                                                    tCGCGGATATGGTCATCCAGGCCGAGTACGTCGCAAGGGTCTTCGACGGGAAGAAGGTGGTCTTCGTTGGTGACGGCGATGGGATAGCCCTTTCAGCTGCGCATCTCAAAGCCCAGGAGGTAATCAGCTACGGGCCAAAATCGATCACGCTCCTGGATTTCGATGAGCGAATCGTCAACTAAGTACGGCGATTCGGTGAAAAACTAgca                                                                                         >  1:72942/1‑208 (MQ=255)
                                                                                                                                                                    tCGCGGATATGGTCATCCAGGCCGAGTACGTCGCAAGGGTCTTCGACGGGAAGAAGGTGGTCTTCGTTGGTGACGGCGATGGGATAGCCCTTTCAGCTGCGCATCTCAAAGCCCAGGAGGTAATCAGCTACGGGCCAAAATCGATCACGCCCCTGGATTTCGATGAGCGAATCGTCAACTCAGTACGGCGATGCGGTGCAAAACTAGCGCCCAAGATCGAAATCACTGCCCAt                                                                 >  2:101426/1‑233 (MQ=255)
                                                                                                                                                                     cgcgGATATGGTCATCCAGGCCGAGTACGTCGCAAGGGTCTTCGACGGGAAGAAGGTGGTCTTCGTTGGTGACGGCGATGGGATAGCCCTTTCAGCTGCGCATCTCAAAGCCCAGGAGGTAATCAGCTACGGGCCAAAATCGATCACGCTCCTGGATTTc                                                                                                                                         >  2:1016490/1‑160 (MQ=255)
                                                                                                                                                                       cgGATATGGTCATCCAGGCCGAGTACGTCGCAAGGGTCTTCGACGGGAAGAAGGTGGTCTTCGTTGGTGACGGCGATGGGATAGCCCTTTCAGCTGCGCATCTCAAAGCCCAGGAGGTAATCAGCTACGGGCCAAAATCGATCACGCTCCTGGATTTCGATGAGCGAATCGTCAACTCAGTACGGCGATTCGGTGAAAAACTAGCGCCCAAGATCGAAATCACTGCCCATCTCTAc                                                           <  1:752138/236‑1 (MQ=255)
                                                                                                                                                                         gATATGGTCATCCAGGCCGAGTACGTCGCAAGGGTCTTCGACGGGAAGAAGGTGGTCTTCGTTGGTGACGGCGATGGGATAGCCCTTTCAGCTGCGCATCTCAAAGCCCAGGAGGTAATCAGCTACGGGCCAAAATCGa                                                                                                                                                          >  2:412104/1‑139 (MQ=255)
                                                                                                                                                                         gATATGGTCATCCAGGCCGAGTACGTCGCAAGGGTCTTCGACGGGAAGAAGGTGGTCTTCGTTGGTGACGGCGATGGGATAGCCCTTTCAGCTGCGCATCTCAAAGCCCAGGAGGTAATCAGCTACGGGCCAAAATCGATCACGCTCCTGGATTTCGATGAGCGAATCGTCAACTCAGTACGGCGATTCGGTGAAAAACTAGCGCCCAAGATCGAAATCACTGCCCATCTCTAc                                                           <  2:856690/234‑1 (MQ=255)
                                                                                                                                                                         gATATGGTCATCCAGGCCGAGTACGTCGCAAGGGTCTTCGACGGGAAGAAGGTGGTCTTCGTTGGTGACGGCGATGGGATAGCCCTTTCAGCTGCGCATCTCAAAGCCCAGGAGGTAATCAGCTACGGGCCAAAATCGATCACGCTCCTGGATTTCGATGAGCGAATCGTCAACTCAGTACGGCGATTCGGTGAAAAACTAGCGCCCAAGATCGAAATCACTGCCCATCTCTACa                                                          >  1:934302/1‑235 (MQ=255)
                                                                                                                                                                         gATATGGTCATCCAGGCCGAGTACGTCGCAAGGGTCTTCGACGGGAAGAAGGTGGTCTTCGTTGGTGACGGCGATGGGATAGCCCTTTCAGCTGCGCATCTCAAAGCCCAGGAGGTAATCAGCTACGGGCCAAAATCGATCACGCTCCTGGATTTCGATGAGCGAATCGTCAACTCAGTACGGCGATTCGGTGAAAAACTAGCGCCCAAGATCGAAATCACTGCCCATCTCTACa                                                          >  1:323140/1‑235 (MQ=255)
                                                                                                                                                                                  aTCCAGGCCGAGTACGTCGCAAGGGTCTTCGACGGGAAGAAGGTGGCCTTCGTTGGTGACGGCGATGGGATAGCCCTTTCAGCTGCGCATCTCAAAGCCCAGGAGGTAATCAGCTACGGGCCAAAATCGATc                                                                                                                                                        <  2:130511/132‑1 (MQ=255)
                                                                                                                                                                                      aGGCCGAGTACGTCGCAAGGGTCTTCGACGGGAAGAAGGTGGTCTTCGTTGGTGACGGCGATGGGATAGCCCTTTCAGCTGCGCATCTCAAAGCCCAGGAGGTAATCAGCTACGGGCCAAAATCGATCACg                                                                                                                                                     <  2:657047/131‑1 (MQ=255)
                                                                                                                                                                                      aGGCCGAGTACGTCGCAAGGGTCTTCGACGGGAAGAAGGTGGTCTTCGTTGGTGACGGCGATGGGATAGCCCTTTCAGCTGCGCATCTCAAAGCCCAGGAGGTAAACAGCTACGGGCCAAAATCGATCACg                                                                                                                                                     >  1:657050/1‑131 (MQ=255)
                                                                                                                                                                                           gAGTACGTCGCAAGGGACTTCGACGGGAAGAAGGTGGTCTTCGTTGGTGACGGCGATGGGATAGCCCTGTCAGCTGCGCATCTCAAAGCCCAGGAGGTAATCAGCTACGGGCCAAAATCGATCACGCTCCTGGATTTCGa                                                                                                                                       >  1:855527/1‑140 (MQ=255)
                                                                                                                                                                                                cGTCGCAAGGGTCTTCGACGGGAAGAAGGTGGTCTTCGTTGGTGACGGCGATGGGATAGCCCTTTCAGCTGCGCATCTCAAAGCCCAGGAGGTAATCAGCTACGGGCCAAAATCGATCACGCTCCTGGATTTCGATGAGCGAATCGTCAACTCAGTACGGCGATTCGGTGAAAAACTAGCGCCCAAGATCGAAATCACt                                                                       <  2:1053163/199‑1 (MQ=255)
                                                                                                                                                                                                    gCAAGGGTCTTCGACGGGAAGAAGGTGGTCTTCGTTGGTGACGGCGATGGGATAGCCCTTTCAGCTGCGCATCTCAAAGCCCAGGAGGTAATCAGCTACGGGCCAAAATCGATCACGCTCCTGGATTTCGATGAGCGAATCGTCAACTCAGTACGGCGATTCGGTGAAAAACTa                                                                                            >  2:733032/1‑174 (MQ=255)
                                                                                                                                                                                                          gTCTTCGACGGGAAGAAGGTGGTCTTCGTTGGTGACGGCGATGGGATAGCCCTTTCAGCTGCGCATCTCAAAGCCCAGGAGGTAATCAGCTACGGGCCAAAATCGATCACGCTCCTGGATTTCGATGAGCGAATCGTCAACTCAGTACGGCGATTCGGTGAAAAACTAGCGCCCAAGATCGAAATCACTGCCCATCTCTACAACGTGGCTGACCCTCTACCTGCCGAACATTTTgg                        >  1:307203/1‑236 (MQ=255)
                                                                                                                                                                                                               cGACGGGAAGAAGGTGGTCTTCGTTGGTGACGGCGATGGGATAGCCCTTTCAGCTGCGCATCTCAAAGCCCAGGAGGTAATCAGCTACGGGCCAAAATCGATCACGCTCCTGGATTTCGATGAGCGAATCGTCAACTCAGTACGGCGATTCGGTGAAAAACTAGCGCCCAAGATCGAAATCACTGCCCATCTCTACAACGTGGCTGACCCTCTACCTGcc                                   >  2:539733/1‑220 (MQ=255)
                                                                                                                                                                                                                gACGGGAAGAAGGTGGTCTTCGTTGGTGACGGCGATGGGATAGCCCTTTCAGCTGCGCATCTCAAAGCCCAGGAGGTAATCAGCTACGGGCCAAAATCGATCACGCTCCTGGATTTCGATGAGCGAATCGTCAACTCAGTACGGCGATTCGGTGAAAAACTAGCGCCCAAGATCGAAATCACTGCCCATCTCTACAACGTg                                                     >  2:893637/1‑201 (MQ=255)
                                                                                                                                                                                                                  aCGGGAAGAAGGTGGTCTTCGTTGGTGACGGCGATGGGATAGCCCTTTCAGCTGCGCATCTCAAAGCCCAGGAGGTAATCAGCTACGGGCCAAAATCGATCACGCTCCTGGATTTCGATGAGCGAATCGTCAACTCAGTACGGCGATTCGGTGAAAAACTAGCGCCCAAGATCGAAATCACTGCCCATCTCTACAACGTGGCTGACCCTCTACCTGCCGAACATTTTGGGGCCTa                  <  1:407232/235‑1 (MQ=255)
                                                                                                                                                                                                                       aagaagGTGGTCTTCGTTGGTGACGGCGATGGGATAGCCCTTTCAGCTGCGCATCTCAAAGCCCAGGAGGTAATCAGCTACGGGCCAAAATCGATCACGCTCCTGGATTTCGATGAGCGAATCGTCAACTCAGTACGGCGATTCGGTGAAAAACTAGCGCCCAATATCGAAATCACTGCCCATCTCTACAACGTGGCTGACCCTCTACCTGCCGAACAtt                            >  2:591125/1‑220 (MQ=255)
                                                                                                                                                                                                                                gtCTTCGTTGGTGACGGCGATGGGATAGCCCTTTCAGCTGCGCATCTCAAAGCCCAGGAGGTAATCAGCTACGGGCCAAAATCGATCCCGATCTTGGATGTCGATGAGCGAATCGTCACCTCAATACGGCGCTTCGGTGAAAAGATAGCGCCCag                                                                                    >  2:526162/1‑154 (MQ=255)
                                                                                                                                                                                                                                 tCTTCGTTGGTGACGGCGATGGGATAGCCCTTTCAGCTGCGCATCTCAAAGCCCAGGAGGTAATCAGCTACGGGCCAAAATCGATCACGCTCCTGGATTTCGATGAGCGAATCGTCAACTCAGTACGGCGATTCGGTGAAAAACTAGCGCCCAAGATCGAAATCACTGCCCATCTCTACAACGTGGCTGACCCTCTACCTGCCGAACATTTTGGGGCCTACGACGGCTTTCATATc  <  1:381358/236‑1 (MQ=255)
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                           
GCAGGTAGCCCTTAGCCATTAGCTGTCCATCAACGCTAAGGATCAGCCGCTGACATGGAGGTTGATAGTGAAGTTCGATCTAAGGAATGCGATCAATAGCATATCTGACGTGGTGATGAATCGACCCGCCCCCTTGCGCGAGTTCGACCAGATCTACATGAAGGTCGCGGATATGGTCATCCAGGCCGAGTACGTCGCAAGGGTCTTCGACGGGAAGAAGGTGGTCTTCGTTGGTGACGGCGATGGGATAGCCCTTTCAGCTGCGCATCTCAAAGCCCAGGAGGTAATCAGCTACGGGCCAAAATCGATCACGCTCCTGGATTTCGATGAGCGAATCGTCAACTCAGTACGGCGATTCGGTGAAAAACTAGCGCCCAAGATCGAAATCACTGCCCATCTCTACAACGTGGCTGACCCTCTACCTGCCGAACATTTTGGGGCCTACGACGGCTTTCATATC  >  AE015451/6150929‑6151388

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