Predicted mutation
evidence seq id position mutation annotation gene description
RA AE015451_phaM1 2,586,629 Δ1 bp intergenic (+535/‑7) PP_t39 tRNA (tRNA‑Leu) → / → PP_2266 CDS (DNA‑directed RNA polymerase_ bacteriophage and mitochondrial type) –/–

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451_phaM12,586,6280T.99.9% 4209.8 / ‑3.9 797intergenic (+534/‑8)PP_t39 tRNA (tRNA‑Leu)/PP_2266 CDS (DNA‑directed RNA polymerase_ bacteriophage and mitochondrial type)–/–
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base . (372/424);  minor base G (1/0);  total (373/424)
Fisher's exact test for biased strand distribution p-value = 4.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 798 total aligned reads displayed.

CCGGGGCCGCGCCAGGGCGCTGTGGGAGACCACATCCAGAGTCGCGGACCCACAATTGGATTCCTTAGCGATGCCCTCTGAGATGACCCTTAGAGACATCCTTAGAAATGATTACTTAATGATGATTACTAAGAGATTCCCTTTATGACTACCTCTATGACTCCTCTAGTTGGCCGCCTGGATGGTGAAGACCTCATCGCCACCCAAGTCCGTCTGGAAGAGGAGATGACCCAGCGCGGTGTGAGCCGTTTCCAGCGCAACGTAGAGGC  >  AE015451_phaM1/2586492‑2586760
                                                                                                                                        |                                                                                                                                    
        gcgcCAGGGCGCTGTGGGAGACCACATCCAGAGTCGCGGACCCACAATTGGATTCCTTAGCGATGCCCTCTGAGATGACCCTTAGAGACATCCTTAGAAATGATTACTTAATGATGATTACTAAGAGATCCCTTTtt                                                                                                                             >  2:3903862‑M1/1‑135 (MQ=255)
         cgcCAGGGCGCTGTGGGAGACCACATCCAGAGGCGCGGACCCACAATTGGATTCCTCAGCGATGCCCTCTGAGAGGACCCTTAAAGACATCCTTAGAAATGATTACTTAATGATGATTACTAAGAGAGCCCTTTttg                                                                                                                            <  2:17655515‑M1/137‑4 (MQ=255)
          gcCAGGGCGCTGTGGGAGACCACATCCAGAGTCGCGGACCCACAATTGGATTCCTTAGCGATGCCCTCTGAGATGACCCTTAGAGACATCCTTAGAAATGATTACTTAATGATGATTACTAAGAGATCCCTTTttgg                                                                                                                           >  1:6239037‑M1/1‑133 (MQ=255)
          gcCAGGGCGCTGTGGGAGACCACATCCAGAGTCGCGGACCCACAATTGGATTCCTTAGCGATGCCCTCTGAGATGACCCTTAGAGACATCCTTAGAAATGATTACTTAATGATGATTACTAAGAGATCCCTTTttgg                                                                                                                           >  2:10534859‑M1/1‑133 (MQ=255)
            cAGGGCGCTGTGGGAGACCACATCCAGAGTCGCGGACCCACAATTGGATTCCTTAGCGATGCCCTCTGAGATGACCCTTAGAGACATCCTTAGAAATGCTTACTTAATGATGGTTACTAAGAGATCCCTTTttggga                                                                                                                         >  1:1416353‑M1/1‑131 (MQ=255)
            cAGGGCGCTGTGGGAGACCACATCCAGAGTCGCGGACCCACAATTGGATTCCTTAGCGATGCCCTCTGAGATGACCCTTAGAGACATCCTTAGAAATGATTACTTAATGATGATTACTAAGAGATCCCTTTttggga                                                                                                                         >  2:7037446‑M1/1‑131 (MQ=255)
            cAGGGCGCTGTGGGAGACCACATCCAGAGTCGCGGACCCACAATTGGATTCCTTAGCGATGCCCTCTGAGATGACCCTTAGAGACATCCTTAGAAATGATTACTTAATGATGATTACTAAGAGATCCCTTTttggga                                                                                                                         <  2:847595‑M1/137‑7 (MQ=255)
            cAGGGCGCTGTGGGAGACCACATCCAGAGTCGCGGACCCACAATTGGATTCCTTAGCGATGCCCTCTGAGATGACCCTTAGAGACATCCTTAGAAATGATTACTTAATGATGATTACTAAGAGATCCCTTTttggga                                                                                                                         <  1:16350415‑M1/137‑7 (MQ=255)
            cAGGGCGCTGTGGGAGACCACATCCAGAGTCGCGGACCCACAATTGGATTCCTTAGCGATGCCCTCTGAGATGACCCTTAGAGACATCCTTAGAAATGATTACTTAATGATGATTACTAAGAGATCCCTTTttggga                                                                                                                         >  2:8725203‑M1/1‑131 (MQ=255)
            cAGGGCGCTGTGGGAGACCACATCCAGAGTCGCGGACCCACAATTGGATTCCTTAGCGATGCCCTCTGAGATGACCCTTAGAGACATCCTTAGAAATGATTACTTAATGATGATTACTAAGAGATCCCTTTttggga                                                                                                                         >  2:770668‑M1/1‑131 (MQ=255)
            cAGGGCGCTGTGGGAGACCACATCCAGAGTCGCGGACCCACAATTGGATTCCTTAGCGATGCCCTCTGAGATGACCCTTAGAGACATCCTTAGAAATGATTACTTAATGATGATTACTAAGAGATCCCTTTttggga                                                                                                                         >  2:17073360‑M1/1‑131 (MQ=255)
             aGGGCGCTGTGGGAGACCACATCCAGAGTCGCGGACCCACAATTGGATTCCTTAGCGATGCCCTCTGAGATGACCCTTAGAGACATCCTTAGAAATGATTACTTAATGATGATTACTAAGAGATCCCTTTttgggac                                                                                                                        >  2:13784577‑M1/1‑130 (MQ=255)
                  gcTGTGGGAGACCACATCCAGAGTCGCGGACCCACAATTGGATTCCTTAGCGATGCCCTCTGAGATGACCCTTAGAGACATCCTTAGAAATGATTACTTAATGATGATTACTAAGAGATCCCTTTttg                                                                                                                            >  2:17356853‑M1/1‑125 (MQ=255)
                             ccACATCCAGAGTCGCGGACCCACAATTGGATTCCTTAGCGATGCCCTCTGAGATGACCCTTAGAGACATCCTTAGAAATGATTACTTAATGATGATTACTAAGAGATCCCTTTttgggacacattttctgattcgt                                                                                                        <  2:5342715‑M1/137‑24 (MQ=255)
                             ccACATCCAGAGTCGCGGACCCACAATTGGATTCCTTAGCGATGCCCTCTGAGATGACCCTTAGAGACATCCTTAGAAATGATTACTTAATGATGATTACTAAGAGATCCCTTTttgggacacattttctgattcgt                                                                                                        <  1:15938471‑M1/137‑24 (MQ=255)
                             ccACATCCAGAGTCGCGGACCCACAATTGGATTCCTTAGCGATGCCCTCTGAGATGACCCTTAGAGACATCCTTAGAAATGATTACTTAATGATGATTACTAAGAGATCCCTTTttgggacacattttctgattcgt                                                                                                        <  1:2151813‑M1/137‑24 (MQ=255)
                             ccACATCCAGAGTCGCGGACCCACAATTGGATTCCTTAGCGATGCCCTCTGAGATGACCCTTAGAGACATCCTTAGAAATGATTACTTAATGATGATTACTAAGAGATCCCTTTttgggacacattttctgattcgt                                                                                                        <  2:3555850‑M1/137‑24 (MQ=255)
                                 aTCCAGAGTCGCGGACCCACAATTGGATTCCTTAGCGATGCCCTCTGAGATGACCCTTAGAGACATCCTTAGAAATGATTACTTAATGATGATTACTAAGAGATCCCTTTttgggacacatttctgattcgtggctc                                                                                                    <  2:15136172‑M1/137‑28 (MQ=255)
                                  tCCAGAGTCGCGGACCCACAATTGGATTCCTTAGCGATGCCCTCTGAGATGACCCTTAGAGACATCCTTAGAAATGATTACTTAATGATGATTACTAAGAGATCCCTTTttgggacacattttctgattcgtggctc                                                                                                   <  1:15639782‑M1/137‑29 (MQ=255)
                                  tCCAGAGTCGCGGACCCACAATTGGATTCCTTAGCGATGCCCTCTGAGATGACCCTTAGAGACATCCTTAGAAATGATTACTTAATGATGATTACTAAGAGATCCATTTttgggacacattttctgattcgtggctc                                                                                                   <  1:10065866‑M1/137‑29 (MQ=255)
                                   ccATAGTCGCGGACCCACAAGTGGATTCCTTAGCTATGCCCCCTGAGATGACCCTTATAGACATCCTTAGAAATGATTACTTAATTATGATTACTAAGATATCCCTTTttgggacacattttctgattcgtggctcg                                                                                                  <  2:15715778‑M1/137‑30 (MQ=255)
                                   ccAGAGTCGCGGACCCACAATTGGATTCCTTAGCGATGCCCTCTGAGATGACCCTTAGAGACATCCTTAGAAATGATTACTTAATGATGATTACTAAGAGATCCCTTTttgggacacattttctgattcgtggctcg                                                                                                  >  2:2982981‑M1/1‑108 (MQ=255)
                                   ccAGAGTCGCGGACCCACAATTGGATTCCTTAGCGATGCCCTCTGAGATGACCCTTAGAGACATCCTTAGAAATGATTACTTAATGATGATTACTAAGAGATCCCTTTttgggacacattttctgattcgtggctcg                                                                                                  >  1:7310016‑M1/1‑108 (MQ=255)
                                   ccAGAGTCGCGGACCCACAATTGGATTCCTTAGCGATGCCCTCTGAGATGACCCTTAGAGACATCCTTAGAAATGATTACTTAATGATGATTACTAAGAGATCCCTTTttgggacacattttctgattcgtggctcg                                                                                                  >  1:11956901‑M1/1‑108 (MQ=255)
                                   ccAGAGTCGCGGACCCACAATTGGATTCCTTAGCGATGCCCTCTGAGATGACCCTTAGAGACATCCTTAGAAATGATTACTTAATGATGATTACTAAGAGATCCCTTTttgggacacattttctgattcgtggctcg                                                                                                  >  2:14608422‑M1/1‑108 (MQ=255)
                                   ccAGAGTCGCGGACCCACAATTGGATTCCTTAGCGATGCCCTCTGAGATGACCCTTAGAGACATCCTTAGAAATGATTACTTAATGATGATTACTAAGAGATCCCTTTttgggacacattttctgattcgtggctcg                                                                                                  >  1:9252241‑M1/1‑108 (MQ=255)
                                   ccAGAGTCGCGGACCCACAATTGGATTCCTTAGCGATGCCCTCTGAGATGACCCTTAGAGACATCCTTAGAAATGATTACTTAATGATGATTACTAAGAGATCCCTTTttgggacacattttctgattcgtggctcg                                                                                                  >  2:5587917‑M1/1‑108 (MQ=255)
                                   ccAGAGTCGCGGACCCACAATTGGATTCCTTAGCGATGCCCTCTGAGATGACCATTAGAGACATCCTTAAAAATGATTACTTAATGATGATTACTAAGAGATCCCTTTttgggacacattttctgattcgtggctcg                                                                                                  >  2:9587186‑M1/1‑108 (MQ=255)
                                    cAGAGTCGCGGACCCACAATTGGATTCCTTAGCGATGCCCTCTGAGATGACCCTTAGAGACATCCTTAGAAATGATTACTTAATGATGATTACTAAGAGATCCCTTTttgggacacattttctgattcgtggctcgg                                                                                                 >  1:14891361‑M1/1‑107 (MQ=255)
                                       agTCGCGGACCCACAATTGGATTCCTTAGCGATGCCCTCTGAGATGACCCTTAGAGACATCCTTAGAAATGATTACTTAATGATGATTACTAAGAGATCCCTTTttgggacacattttctgattcgtggctcgggca                                                                                              <  1:5861910‑M1/137‑34 (MQ=255)
                                          cgcgGACCCACAATTGGATTCCTTAGCGATGCCCTCTGAGATGACCCTTAGAGACATCCTTAGCAATGATTACTTAATGATGATTACTAAGAGATCCCTTTttgggacacattttctgattcgtggctcgggcagag                                                                                           >  1:66790‑M1/1‑101 (MQ=255)
                                          cgcgGACCCACAATTGGATTCCTTAGCGATGCCCTCTGAGATGACCCTTAGAGACATCCTTAGAAATGATTACTTAATGATGATTACTAAGAGATCCCTTTttgggacacattttctgattcgtggctcgggcagag                                                                                           >  2:15424813‑M1/1‑101 (MQ=255)
                                          cgcgGACCCAAAATTGGATTCCTTAGCGACGCCCTCTGAGATGACCCTTAGAAACATCCGTAGCCATGATTACTCAATGATGATTAATAAGATATCCCTTTttgggacacattttctgattcgtggctcggacagag                                                                                           >  2:12640333‑M1/1‑101 (MQ=255)
                                            cgGACCCACAATTGGATTCCTTAGCGATGCCCTCTGAGATGACCCTTAGAGACATCCTTAGAAATGATTACTTAATGATGATTACTAAGAGATCCCTTTttgggacacattttctgattcgtggctcgggcagagta                                                                                         >  2:10679866‑M1/1‑99 (MQ=255)
                                              gACCCACAATTGGATTCCTTAGCGATGCCCTCTGAGATGACCCTTAGAGACATCCTTAGAAATGATTACTTAATGATGATTACTAAGAGATCCCTTTttgggacacattttctgattcgtggctcgggcagagtagg                                                                                       <  1:8407677‑M1/137‑41 (MQ=255)
                                              gACCCACAATTGGATTCCTTAGCGATGCCCTCTGAGATGACCCTTAGAGACATCCTTAGAAATGATTACTTAATGATGATTACTAAGAGATCCCTTTttgggacacattttctgattcgtggctcgggcagagtagg                                                                                       >  2:3270097‑M1/1‑97 (MQ=255)
                                               aCCCACAATTGGATTCCTTAGCGATGCCCTCTGAGATGACCCTTAGAGACATCCTTAGAAATGATTACTTAATGATGATTACTAAGAGATCCCTTTttgggacacattttctgattcgtggctcgggcagagtaggc                                                                                      >  1:12457284‑M1/1‑96 (MQ=255)
                                                     aaTTGGATTCCTTAGCGATGCCCTCTGAGATGACCCTTAGAGACATCCTTAGAAATGATTACTTAATGATGATTACTAAGAGATCCCTTTttgggacacattttctgattcgtggctcgggcagagtaggccctgag                                                                                <  2:5164759‑M1/137‑48 (MQ=255)
                                                     aaTTGGATTCCTTAGCGATGCCCTCTGAGATGACCCTTAGAGACATCCTTAGAAATGATTACTTAATGATGATTACTAAGAGATCCCTTTttgggacacattttctgattcgtggctcgggcagagtaggccctgag                                                                                <  2:238730‑M1/137‑48 (MQ=255)
                                                      aTTGGATTCCTTAGCGATGCCCTCTGAGATGACCCTTATAGACATCCTTAGAAATGATTACTTAATGATGATTACTAAGAGATCCCTTTttgggacacattttctgattcgtggctcgggcagagtaggccctgagc                                                                               <  1:16892266‑M1/137‑49 (MQ=255)
                                                      aTTGGATTCCTTAGCGATGCCCTCTGAGATGACCCTTAGAGACATCCTTAGAAATGATTACTTAATGATGATTACTAAGAGATCCCTTTttgggacacattttctgattcgtggctcgggcagagtaggccctgagc                                                                               <  2:17144074‑M1/137‑49 (MQ=255)
                                                      aTTGGATTCCTTAGCGATGCCCTCTGAGATGACCCTTAGAGACATCCTTAGAAATGATTACTTAATGATGATTACTAAGAGATCCCTTTttgggacacattttctgattcgtggctcgggcagagtaggccctgagc                                                                               <  1:171716‑M1/137‑49 (MQ=255)
                                                      aTTGGATTCCTTAGCGATGCCCTCTGAGATGACCCTTAGAGACATCCTTAGAAATGATTACTTAATGATGATTACGAAGAGATCCCTTTtttggacacattttctgattcgtggctcgggcagagtaggccctgagc                                                                               <  1:3215046‑M1/137‑49 (MQ=255)
                                                              ccTTAGCGATGCCCTCTGAGATGACCCTTAGAGACATCCTTAGAAATGATTACTTAATGATGATTACTAAGAGATCCCTTTttgggacacattttctgattcgtggctcgggcagagtaggccctgagcggcgaccg                                                                       >  1:5663633‑M1/1‑81 (MQ=255)
                                                                ttAGCGATGCCCTCTGAGATGACCCTTAGAGACATCCTTAGAAATGATTACTTAATGATGATTACTAAGAGATCCCTTTttgggacacattttctgattcgtggctcgggcagagtaggccctgagcggcgaccgtg                                                                     <  2:16909513‑M1/137‑59 (MQ=255)
                                                                ttAGCGATGCCCTCTGAGATGACCCTTAGAGACATCCTTAGAAATGATTACTTAATGATGATTACTAAGAGATCCCTTTttgggacacattttctgattcgtggctcgggcagagtaggccctgagcggcgaccgtg                                                                     <  1:8104362‑M1/137‑59 (MQ=255)
                                                                ttAGCGATGCCCTCTGAGATGACCCTTAGAGACATCCTTAGAAATGATTACTTAATGATGATTACTAAGAGATCCCTTTttgggacacattttctgattcgtggctcgggcagagtaggccctgagcggcgaccgtg                                                                     <  1:11306495‑M1/137‑59 (MQ=255)
                                                                 tAGCGATGCCCTCTGAGATGACCCTTAGAGACATCCTTAGAAATGATTACTTAATGATGATTACTAAGAGATCCCTTTttgggacacattttctgattcgtggctcgggcagagtaggccctgagcggcgaccgtgg                                                                    <  1:637294‑M1/137‑60 (MQ=255)
                                                                 tAGCGATGCCCTCTGAGATGACCCTTAGAGACATCCTTAGAAATGATTACTTAATGATGATTACTAAGAGATCCCTTTttgggacacattttctgattcgtggctcgggcagagtaggccctgagcggcgaccgtgg                                                                    <  2:14872447‑M1/137‑60 (MQ=255)
                                                                 tAGCGATGCCCTCTGAGATGACCCTTAGAGACATCCTTAGAAATGATTACTTAATGATGATTACTAAGAGATCCCTTTttgggacacattttctgattcgtggctcgggcagagtaggccctgagcggcgaccgtgg                                                                    <  2:14385633‑M1/137‑60 (MQ=255)
                                                                  aGCGATGCCCTCTGAGATGACCCTTAGAGACATCCTTAGAAATGATTACTTAATGATGATTACTAAGAGATCCCTTTttgggacacattttctgattcgtggctcgggcagagtaggccctgagcggcgaccgtggg                                                                   >  2:7927897‑M1/1‑77 (MQ=255)
                                                                    cGATGCCCTCTGAGATGACCCTTAGAGACATCCTTAGAAATGATTACTTAATGATGATTACTAAGAGATCCCTTTttgggacacattttctgattcgtggctcgggcagagtaggccctgagcggcgaccgtgggcc                                                                 >  1:479991‑M1/1‑75 (MQ=255)
                                                                          ccTCTGAGATGACCCTTAGAGACATCCTTAGAAATGATTACTTAATGATGATTACTAAGAGATCCCTTTttgggacacattttctgattcgtggctcgggcagagtaggccctgagcggcgaccgtgggcccaccca                                                           <  2:11065800‑M1/137‑69 (MQ=255)
                                                                          ccTCTGAGATGACCCTTAGAGACATCCTTAGAAATGATTACTTAATGATGATTACTAAGAGATCCCTTTttgggacacattttctgattcgtggctcgggcagagtaggccctgagcggcgaccgtgggcccaccca                                                           <  1:13784577‑M1/137‑69 (MQ=255)
                                                                          ccTCTGAGATGACCCTTAGAGACATCCTTAGAAATGATTACTTAATGATGATTACTAAGAGATCCCTTTttgggacacattttctgattcgtggctcgggcagagtaggccctgagcggcgaccgtgggcccaccca                                                           <  1:1177903‑M1/137‑69 (MQ=255)
                                                                            tctGAGATGACCCTTAGAGACATCCTTAGAAATGATTACTTAATGATGATTACTAAGAGATCCCTTTttgggacacattttctgattcgtggctcgggcagagtaggccctgagcggcgaccgtgggcccacccaaa                                                         >  1:5350198‑M1/1‑67 (MQ=255)
                                                                            tctGAGATGACCCTTAGAGACATCCTTAGAAATGATTACTTAATGATGATTACTAAGAGATCCCTTTttgggacacattttctgattcgtggctcgggcagagtaggccctgagcggcgaccgtgggcccacccaaa                                                         >  1:1048118‑M1/1‑67 (MQ=255)
                                                                               gagaTGACCCTTAGAGACATCCTTAGAAATGATTACTTAATGATGATTACTAAGAGATCCCTTTttgggacacattttctgattcgtggctcgggcagagtaggccctgagcggcgaccgtgggcccacccaaaacg                                                      <  2:11776371‑M1/137‑74 (MQ=255)
                                                                               gagaTGACCCTTAGAGACATCCTTAGAAATGATTACTTAATGATGATTACTAAGAGATCCCTTTttgggacacattttctgattcgtggctcgggcagagtaggccctgagcggcgaccgtgggcccacccaaaacg                                                      <  1:8660527‑M1/137‑74 (MQ=255)
                                                                               gagaTGACCCTTAGAGACATCCTTAGAAATGATTACTTAATGATGATTACTAAGAGATCCCTTTttgggacacattttctgattcgtggctcgggcagagtaggccctgagcggcgaccgtgggcccacccaaaacg                                                      <  1:403045‑M1/137‑74 (MQ=255)
                                                                                 gaTGACCCTTAGAGACATCCTTAGAAATGATTACTTAATGATGATTACTAAGAGATCCCTTTttgggacacattttctgattcgtggctcgggcagagtaggccctgagcggcgaccgtgggcccacccaaaacggg                                                    >  1:10828165‑M1/1‑62 (MQ=255)
                                                                                    gACCCTTAGAGACATCCTTAGAAATGATTACTTAATGATGATTACTAAGAGATCCCTTTttgggacacattttctgattcgtggctcgggcagagtaggccctgagcggcgaccgtgggcccacccaaaacggggcg                                                 <  2:4509394‑M1/137‑79 (MQ=255)
                                                                                    gACCCTTAGAGACATCCTTAGAAATGATTACTTAATGATGATTACTAAGAGATCCCTTTttgggacacattttctgattcgtggctcgggcagagtaggccctgagcggcgaccgtgggcccacccaaaacggggcg                                                 <  1:4989977‑M1/137‑79 (MQ=255)
                                                                                    gACCCTTAGAGACATCCTTAGAAATGATTACTTAATGATGATTACTAAGAGATCCCTTTttgggacacattttctgattcgtggctcgggcagagtaggccctgagcggcgaccgtgggcccacccaaaacggggcg                                                 <  2:13381512‑M1/137‑79 (MQ=255)
                                                                                    gACCCTTAGAGACATCCTTAGAAATGATTACTTAATGATGATTACTAAGAGATCCCTTTttgggacacattttctgattcgtggctcgggcagagtaggccctgagcggcgaccgtgggcccacccaaaacggggcg                                                 <  2:12057602‑M1/137‑79 (MQ=255)
                                                                                      cccTTAGAGACATCCTTAGAAATGATTACTTAATGATGATTACTAAGAGATCCTTTTttgggacacattttctgattcgtggctcgggcagagtaggccctgagcggcgaccgtgggcccacccaaaacggggcggg                                               <  1:4544562‑M1/137‑81 (MQ=255)
                                                                                      cccTTAGAGACATCCTTAGAAATGATTACTTAATGATGATTACTAAGAGATCCCTTTttgggacacattttctgattcgtggctcgggcagagtaggccctgagcggcgaccgtgggcccacccaaaacggggcggg                                               <  2:12339600‑M1/137‑81 (MQ=255)
                                                                                       ccTTAGAGACATCCTTAGAAATGATTACTTAATGATGATTACTAAGAGATCCCTTTttgggacacattttctgattcgtggctcgggcagagtaggccctgagcggcgaccgtgggcccacccaaaacggggcggga                                              >  2:1994207‑M1/1‑56 (MQ=255)
                                                                                            gagaCATCCTTAGAAATGATTACTTAATGATGATTACTAAGAGATCCCTTTttgggacacattttctgattcgtggctcgggcagagtaggccctgagcggcgaccgtgggcccacccaaaacggggcgggacacac                                         <  1:12789196‑M1/137‑87 (MQ=255)
                                                                                              gaCATCCTTAGAAATGATTACTTAATGATGATTACTAAGAGATCCCTTTttgggacacattttctgattcgtggctcgggcagagtaggccctgagcggcgaccgtgggcccacccaaaacggggcgggacacaccg                                       <  1:1871879‑M1/137‑89 (MQ=255)
                                                                                                cATCCTTAGAAATGATTACTTAATGATGATTACTAAGAGATCCCTTTttgggacacattttctgattcgtggcgcgggcagagtaggccctgagcggcgaccgtgggcccacccaaaacggggcgggacacaccggg                                     <  1:7769853‑M1/137‑91 (MQ=255)
                                                                                                 aTCCTTAGAAATGATTACTTAATGATGATTACTAAGAGATCCCTTTttgggacacattttctgattcgtggctcgggcagagtaggccctgagcggcgaccgtgggcccacccaaaacggggcgggacacaccggga                                    <  2:16523795‑M1/137‑92 (MQ=255)
                                                                                                  tCCTTAGAAATGATTACTTAATGATGATTACTAAGAGATCCCTTTttgggacacattttctgattcgtggctcgggcagagtaggccctgagcggcgaccgtgggcccacccaaaacggggcgggacacaccgggac                                   <  2:6231176‑M1/137‑93 (MQ=255)
                                                                                                   ccTTAGAAATGATTACTTAATGATGATTACTAAGAGATCCCTTTttgggacacattttctgattcgtggctcgggcagagtaggccctgagcggcgaccgtgggcccacccaaaacggggcgggacacaccgggaca                                  >  1:3646095‑M1/1‑44 (MQ=255)
                                                                                                    cTTAGAAATGATTACTTAATGATGATTACTAAGAGATCCCTTTttgggacacattttctgattcgtggctcgggcagagtaggccctgagcggcgaccgtgggcccacccaaaacggggcgggacacaccgggacac                                 >  2:1303692‑M1/1‑43 (MQ=255)
                                                                                                    cTTAGAAATGATTACTTAATGATGATTACTAAGAGATCCCTTTttgggacacattttctgattcgtggctcgggcagagtaggccctgagcggcgaccgtgggcccacccaaaacggggcgggacacaccgggacac                                 >  2:10049709‑M1/1‑43 (MQ=255)
                                                                                                    cTTAGAAATGATTACTTAATGATGATTACTAAGAGATCCCTTTttgggacacattttctgattcgtggctcgggcagagtaggccctgagcggcgaccgtgggcccacccaaaacggggcgggacacaccgggacac                                 >  2:9259818‑M1/1‑43 (MQ=255)
                                                                                                      tAGAAATGATTACTTAATGATGATTACTAAGAGATCCCTTTttgggacacattttctgattcgtggctcgggcagagtaggccctgagcggcgaccgtgggcccacccaaaacggggcgggacacaccgggacacat                               <  1:2759520‑M1/137‑97 (MQ=255)
                                                                                                      tAGAAATGATTACTTAATGATGATTACTAAGAGATCCCTTTttgggacacattttctgattcgtggctcgggcagagtaggccctgagcggcgaccgtgggcccacccaaaacggggcgggacacaccgggacacat                               >  1:166910‑M1/1‑41 (MQ=255)
                                                                                                       aGAAATGATTACTTAATGATGATTACTAAGAGATCCCTTTttgggacacattttctgattcgtggctcgggcagagtaggccctgagcggcgaccgtgggcccacccaaaacggggcgggacacaccgggacacatg                              >  2:8206257‑M1/1‑40 (MQ=255)
                                                                                                        gAAATGATTACTTAATGATGATTACTAAGAGATCCCTTTttgggacacattttctgattcgtggctcgggcagagtaggccctgagcggcgaccgtgggcccacccaaaacggggcgggacacaccgggacacatgg                             >  2:9818378‑M1/1‑39 (MQ=255)
                                                                                                               ttACTTAATGATGATTACTAAGAGATCCCTTTttgggacacattttctgattcgtggctcgggcagagtaggccctgagcggcgaccgtgggcccacccaaaacggggcgggacacaccgggacacatggggacact                      >  1:791911‑M1/1‑32 (MQ=255)
                                                                                                                 aCTTAATGATGATTACTAAGAGATCCCTTTttgggacacattttctgattcgtggctcgggcagagtaggccctgagcggcgaccgtgggcccacccaaaacggggcgggacacaccgggacacatggggacacttg                    >  2:4343553‑M1/1‑30 (MQ=255)
                                                                                                                         atgatTACTAAGAGATCCCTTTttgggacacattttctgattcgtggctcgggcagagtaggccctgagcggcgaccgtgggcccacccaaaacggggcgggacacaccgggacacatggggacacttgggacacga            <  1:15313643‑M1/137‑116 (MQ=255)
                                                                                                                         atgatTACTAAGAGATCCCTTTttgggacacattttctgattcgtggctcgggcagagtaggccctgagcggcgaccgtgggcccacccaaaacggggcgggacacaccgggacacatggggacacttgggacacga            <  1:891318‑M1/137‑116 (MQ=255)
                                                                                                                           gatTACTAAGAGATCCCTTTttgggacacattttctgattcgtggctcgggcagagtaggccctgagcggcgaccgtgggcccacccaaaacggggcgggacacaccgggacacatggggacacttgggacacgaat          <  1:13214440‑M1/137‑118 (MQ=255)
                                                                                                                             ttACTAAGAGATCCCTTTttgggacacattttctgattcgtggctcgggcagagtaggccctgagcggcgaccgtgggcccacccaaaacggggcgggacacaccgggacacatggggacacttgggacacgaattt        <  2:66790‑M1/137‑120 (MQ=255)
                                                                                                                             ttACTAAGAGATCCCTTTttgggacacattttctgattcgtggctcgggcagagtaggccctgagcggcgaccgtgggcccacccaaaacggggcgggacacaccgggacacatggggacacttgggacacgaattt        <  1:14374077‑M1/137‑120 (MQ=255)
                                                                                                                             ttACTAAGAGATCCCTTTttgggacacattttctgattcgtggctcgggcagagtaggccctgagcggcgaccgtgggcccacccaaaacggggcgggacacaccgggacacatggggacacttgggacacgaattt        <  2:11326640‑M1/137‑120 (MQ=255)
                                                                                                                             ctaCTAAGAGATCCCTTTttgggacacattttctgattcttggctcggacagagtaggccctgagcggcgaccgtgggcccacccaaaacggggcgggacacaccgggacacatggggacacttgggacacgaattt        <  2:4879819‑M1/136‑120 (MQ=255)
                                                                                                                              tACTAAGAGATCCCTTTttgggacacattttctgattcgtggctcgggcagagtaggccctgagcggcgaccgtgggcccacccaaaacggggcgggacacaccgggacacatggggacacttgggacacgaattta       <  1:2546962‑M1/137‑121 (MQ=255)
                                                                                                                               aCTAAGAGATCCCTTTttgggacacattttctgattcgtggctcgggcagagtaggccctgagcggcgaccgtgggcccacccaaaacggggcgggacacaccgggacacatggggacacttgggacacgaatttag      >  1:982388‑M1/1‑16 (MQ=255)
                                                                                                                               aCTAAGAGATCCCTTTttgggacacattttctgattcgtggctcgggcagagtaggccctgagcggcgaccgtg                                                                     >  2:5803666‑M1/1‑16 (MQ=255)
                                                                                                                                 tAAGAGATCCCTTTttgggacacattttctgattcgtggctcgggcagagtaggccctgagcggcgaccgtgggcccacccaaaacggggcgggacacaccgggacacatggggacacttgggacacgaatttagcg    <  1:466553‑M1/137‑124 (MQ=255)
                                                                                                                                  aaGAGATCCCTTTttgggacacattttctgattcgtggctcgggcagagtaggccctgagcggcgaccgtgggcccacccaaaacggggcgggacacaccgggacacatggggacacttgggacacgaatttagcgt   <  2:14007609‑M1/137‑125 (MQ=255)
                                                                                                                                  aaGAGATCCCTTTttgggacacattctctgattcgtggctcgggcagagtaggccctgagcggcgaccgtgggcccacccaaaacggggcgggacacaccgggacacatggggacacttgggacacgaatttagcgt   <  2:16168247‑M1/137‑125 (MQ=255)
                                                                                                                                   agagATCCCTTTttgggacacattttctgattcgtggctcgggcagagtaggccctgagcggcgaccgtgggcccacccaaaacggggcgggacacaccgggacacatggggacacttgggacacgaatttagcgtt  <  2:2759095‑M1/137‑126 (MQ=255)
                                                                                                                                   agagATCCCTTTttgggacacattttctgattcgtggctcgggcagagtaggccctgagcggcgaccgtgggcccacccaaaacggggcgggacacaccgggacacatggggacacttgggacacgaatttagcgtt  <  1:420777‑M1/137‑126 (MQ=255)
                                                                                                                                   agagATCCCTTTttgggacacattttctgattcgtggctcgggcagagtaggccctgagcggcgaccgtgggcccacccaaaacggggcgggacacaccgggacacatggggacacttgggacacgaatttagcgtt  <  1:17358073‑M1/137‑126 (MQ=255)
                                                                                                                                   agagATCCCTTTttgggacacattttctgattcgtggctcgggcagagtaggccctgagcggcgaccgtgggcccacccaaaacggggcgggacacaccgggacacatggggacacttgggacacgaatttagcgtt  <  1:10273735‑M1/137‑126 (MQ=255)
                                                                                                                                        |                                                                                                                                    
CCGGGGCCGCGCCAGGGCGCTGTGGGAGACCACATCCAGAGTCGCGGACCCACAATTGGATTCCTTAGCGATGCCCTCTGAGATGACCCTTAGAGACATCCTTAGAAATGATTACTTAATGATGATTACTAAGAGATTCCCTTTATGACTACCTCTATGACTCCTCTAGTTGGCCGCCTGGATGGTGAAGACCTCATCGCCACCCAAGTCCGTCTGGAAGAGGAGATGACCCAGCGCGGTGTGAGCCGTTTCCAGCGCAACGTAGAGGC  >  AE015451_phaM1/2586492‑2586760

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