Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,198,675 A→G 5.5% intergenic (‑438/‑20) ymfE ← / → lit e14 prophage; putative inner membrane protein/T4 phage exclusion protein; cell death peptidase, e14 prophage

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,198,6750AG 5.5% 230.2 / 3.5 74intergenic (‑438/‑20)ymfE/lite14 prophage; putative inner membrane protein/T4 phage exclusion protein; cell death peptidase, e14 prophage
Reads supporting (aligned to +/- strand):  ref base A (38/31);  new base G (2/2);  total (41/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.64e-01

GTCTTACGACATGCTTACCTTAATCGTTCCCAGTGTGCCTATAGCAGATATTCTAAACATGTCGATAATTCATTACGCATATAGTATCGAACATAGAAAAAACTGAAGATTCATCTTATTTTGTATATACTACCTAGCCCAACAATGTAGAGGTTAACGAAAAATGCGCTCACCAATTTGTCATCTTTTCTCAGCAATTAATTCATCACCATTTAAGATTGCACCAGAGAAGGAGCAAGATCTTAAAACGATAGTTGACGACAAAAAAATTATAATTTCAGTTGTGAGT  >  NC_000913/1198532‑1198820
                                                                                                                                               |                                                                                                                                                 
gtCTTACGACATGCTTACCTTAATCGTTCCCAGTGTGCCTATAGCAGATATTCTAAACATGTCGATAATTCATTACGCATATAGTATCGAACATAGAAAAAACTGAAGATTCATCTTATTTTGTATATACTACCTAGCCCAACGATGTa                                                                                                                                              >  3:3944/1‑149 (MQ=255)
gtCTTACGACATGCTTACCTTAATCGTTCCCAGTGTGCCTATAGCAGATATTCTAAACATGTCGATAATTCATTACGCATATAGTATCGAACATAGAAAAAACTGAAGATTCATCTTATTTTGTATATACTACCTAGCCCAACGATGTa                                                                                                                                              >  1:123876/1‑149 (MQ=255)
  cTTACGACATGCTTACCTTAATCGTTCCCAGTGTGCCTATAGCAGATATTCTAAACATGTCGATAATTCATTACGCATATAGTATCGAACATAGAAAAAACTGAAGATTCATCTTATTTTGTATATACTACCTAGCCCAACAATGTaga                                                                                                                                            <  1:34078/149‑1 (MQ=255)
  cTTACGACATGCTTACCTTAATCGTTCCCAGTGTGCCTATAGCAGATATTCTAAACATGTCGATAATTCATTACGCATATAGTATCGAACATAGAAAAAACTGAAGATTCATCTTATTTTGTATATACTACCTAGCCCAACAACGTaga                                                                                                                                            <  4:436449/149‑1 (MQ=255)
   ttACGACATGCTTACCTTAATCGTTCCCAGTGTGCCTATAGCAGATATTCTAAACATGTCGATAATTCATTACGCATATAGTATCGAACATAGAAAAAACTGAAGATTCATCTTATTTTGTATATACTACCTAGCCCAACAATGTagag                                                                                                                                           <  1:108944/149‑1 (MQ=255)
   ttACGACATGCTTACCTTAATCGTTCCCAGTGTGCCTATAGCAGATATTCTAAACATGTCGATAATTCATTACGCATATAGTATCGAACATAGAAAAAACTGAAGATTCATCTTATTTTGTATATACTACCTAGCCCAACAATGTagag                                                                                                                                           <  4:334147/149‑1 (MQ=255)
   ttACGACATGCTTACCTTAATCGTTCCCAGTGTGCCTATAGCAGATATTCTAAACATGTCGATAATTCATTACGCATATAGTATCGAACATAGAAAAAACTGAAGATTCATCTTATTTTGTATATACTACCTAGCCCAACAATGTagag                                                                                                                                           >  1:291855/1‑149 (MQ=255)
     aCGACATGCTTACCTTAATCGTTCCCAGTGTGCCTATAGCAGATATTCTAAACATGTCGATAATTCATTACGCATATAGTATCGAACATAGAAAAAACTGAAGATTCATCTTATTTTGTATATACTACCTAGCCCAACGATGTAGAGGt                                                                                                                                         <  2:123876/149‑1 (MQ=255)
     aCGACATGCTTACCTTAATCGTTCCCAGTGTGCCTATAGCAGATATTCTAAACATGTCGATAATTCATTACGCATATAGTATCGAACATAGAAAAAACTGAAGATTCATCTTATTTTGTATATACTACCTAGCCCAACGATGTAGAGGt                                                                                                                                         <  4:3944/149‑1 (MQ=255)
         cATGCTTACCTTAATCGTTCCCAGTGTGCCTATAGCAGATATTCTAAACATGTCGATAATTCATTACGCATATAGTATCGAACATAGAAAAAACTGAAGATTCATCTTATTTTGTATATACTACCTAGCCCAACAATGTAGAGGTTAAc                                                                                                                                     <  4:555663/149‑1 (MQ=255)
           tGCTTACCTTAATCGTTCCCAGTGTGCCTATAGCAGATATTCTAAACATGTCGATAATTCATTACGCATATAGTATCGAACATAGAAAAAACTGAAGATTCATCTTATTTTGTATATACTACCTAGCCCAACAATGTAGAGGTTAACGa                                                                                                                                   >  2:248357/1‑149 (MQ=255)
            gCTTACCTTAATCGTTCCCAGTGTGCCTATAGCAGATATTCTAAACATGTCGATAATTCATTACGCATATAGTATCGAACATAGAAAAAACTGAAGATTCATCTTATTTTGTATATACTACCTAGCCCAACAATGTagag                                                                                                                                           <  2:436017/140‑1 (MQ=255)
            gCTTACCTTAATCGTTCCCAGTGTGCCTATAGCAGATATTCTAAACATGTCGATAATTCATTACGCATATAGTATCGAACATAGAAAAAACTGAAGATTCATCTTATTTTGTATATACTACCTAGCCCAACAATGTagag                                                                                                                                           >  1:436017/1‑140 (MQ=255)
                  cTTAATCGTTCCCAGTGTGCCTATAGCAGATATTCTAAACATGTCGATAATTCATTACGCATATAGTATCGAACATAGAAAAAACTGAAGATTCATCTTATTTTGTATATACTACCTAGCCCAACAATGTAGAGGTTAACGAAAAATgc                                                                                                                            >  2:295798/1‑149 (MQ=255)
                  cTTAATCGTTCCCAGTGTGCCTATAGCAGATATTCTAAACATGTCGATAATTCATTACGCATATAGTATCGAACATAGAAAAAACTGAAGATTCATCTTATTTTGTATATACTACCTAGCCCAACAATGTAGAGGTTAACGAAAAATgc                                                                                                                            >  4:496623/1‑149 (MQ=255)
                          ttCCCAGTGTGCCTATAGCAGATATTCTAAACATGTCGATAATTCATTACGCATATAGTATCGAACATAGAAAAAACTGAAGATTCATCTTATTTTGTATATACTACCTAGCCCAACAATGTAGAGGTTAACGAAAAATGCGCTCACCa                                                                                                                    >  4:134442/1‑149 (MQ=255)
                            cccAGTGTGCCTATAGCAGATATTCTAAACATGTCGATAATTCATTACGCATATAGTATCGAACATAGAAAAAACTGAAGATTCATCTTATTTTGTATATACTACCTAGCCCAACAATGTAGAGGTTAACGAAAAATGCGCTCACCAAt                                                                                                                  <  3:479211/149‑1 (MQ=255)
                              cAGTGTGCCTATAGCAGATATTCTAAACATGTCGATAATTCATTACGCATATAGTATCGAACATAGAAAAAACTGAAGATTCATCTTATTTTGTATATACTACCTAGCCCAACAATGTAGAGGTTAACGAAAAATGCGCTCACCAAttt                                                                                                                <  4:50354/149‑1 (MQ=255)
                              cAGTGTGCCTATAGCAGATATTCTAAACATGTCGATAATTCATTACGCATATAGTATCGAACATAGAAAAAACTGAAGATTCATCTTATTTTGTATATACTACCTAGCCCAACAATGTAGAGGTTAACGAAAAATGCGCTCACCAAttt                                                                                                                <  1:490113/149‑1 (MQ=255)
                                 tgtgCCTATAGCAGATATTCTAAACATGTCGATAATTCATTACGCATATAGTATCGAACATAGAAAAAACTGAAGATTCATCTTATTTTGTATATACTACCTAGCCCAACAATGTAGAGGTTAACGAAAAATGCGCTCACCAATTTGTc                                                                                                             <  1:497818/149‑1 (MQ=255)
                                   tgCCTATAGCAGATATTCTAAACATGTCGATAATTCATTACGCATATAGTATCGAACATAGAAAAAACTGAAGATTCATCTTATTTTGTATATACTACCTAGCCCAACAATGTAGAGGTTAACGAAAAATGCGCTCACCAATTTGTCAt                                                                                                           <  1:93942/149‑1 (MQ=255)
                                    gCCTATAGCAGATATTCTAAACATGTCGATAATTCATTACGCATATAGTATCGAACATAGAAAAAACTGAAGATTCATCTTATTTTGTATATACTACCTAGCCCAACAATGTAGAGGTTAACGAAAAATGCGCTCACCAATTTGTCATc                                                                                                          <  2:437668/149‑1 (MQ=255)
                                     ccTATAGCAGATATTCTAAACATGTCGATAATTCATTACGCATATAGTATCGAACATAGAAAAAACTGAAGATTCATCTTATTTTGTATATACTACCTAGCCCAACAATGTAGAGGTTAACGAAAAATGCGCTCACCAATTTGTCATCt                                                                                                         >  2:379054/1‑149 (MQ=255)
                                       tataGCAGATATTCTAAACATGTCGATAATTCATTACGCATATAGTATCGAACATAGAAAAAACTGAAGATTCATCTTATTTTGTATATACTACCTAGCCCAACAATGTAGAGGTTAACGAAAAATGCGCTCACCAATTTGTCATCttt                                                                                                       >  3:152569/1‑149 (MQ=255)
                                         taGCAGATATTCTAAACATGTCGATAATTCATTACGCATATAGTATCGAACATAGAAAAAACTGAAGATTCATCTTATTTTGTATATACTACCTAGCCCAACAATGTAGAGGTTAACGAAAAATGCGCTCACCAATTTGTCATCTTTtc                                                                                                     <  3:451201/149‑1 (MQ=255)
                                         taGCAGATATTCTAAACATGTCGATAATTCATTACGCATATAGTATCGAACATAGAAAAAACTGAAGATTCATCTTATTTTGTATATACTACCTAGCCCAACAATGTAGAGGTTAACGAAAAATGCGCTCACCAATTTGTCATCTTTtc                                                                                                     >  3:347319/1‑149 (MQ=255)
                                          aGCAGATATTCTAAACATGTCGATAATTCATTACGCATATAGTATCGAACATAGAAAAAACTGAAGATTCATCTTATTTTGTATATACTACCTAGCCCAACAATGTAGAGGTTAACGAAAAATGCGCTCACCAATTTGTCATCTTTtct                                                                                                    >  2:332074/1‑149 (MQ=255)
                                           gCAGATATTCTAAACATGTCGATAATTCATTACGCATATAGTATCGAACATAGAAAAAACTGAAGATTCATCTTATTTTGTATATACTACCTAGCCCAACAATGTAGAGGTTAACGAAAAATGCGCTCACCAATTTGTCATCTTTtctc                                                                                                   <  2:491887/149‑1 (MQ=255)
                                            cAGATATTCTAAACATGTCGATAATTCATTACGCATATAGTATCGAACATAGAAAAAACTGAAGATTCATCTTATTTTGTATATACTACCTAGCCCAACAATGTAGAGGTTAACGAAAAATGCGCTCACCAATTTGTCATCTTTTCTCa                                                                                                  >  1:62106/1‑149 (MQ=255)
                                             aGATATTCTAAACATGTCGATAATTCATCACGCATATAGTATCGAACATAGAAAAAACTGAAGATTCATCTTATTTTGTAGATACCACCGAGCCCAACAATGTAGAGGTTAACGAAAAATGCGCTCACCAATTTGTCATCTTTTCTCAg                                                                                                 >  2:417312/1‑149 (MQ=255)
                                              gATATTCTAAACATGTCGATAATTCATTACGCATATAGTATCGAACATAGAAAAAACTGAAGATTCATCTTATTTTGTATATACTACCTAGCCCAACAACGTAGAGGTTAACGAAAAATGCGCTCACCAATTTGTCATCTATTCTCAGc                                                                                                >  4:274301/1‑149 (MQ=255)
                                               atatTCTAAACATGTCGATAATTCATTACGCATATAGTATCGAACATAGAAAAAACTGAAGATTCATCTTATTTTGTATATACTACCTAGCCCAACAATGTAGAGGTTAACGAAAAATGCGCTCACCAATTTGTCATCTTTTCTCAGCa                                                                                               >  1:107499/1‑149 (MQ=255)
                                                    cTAAACATGTCGATAATTCATTACGCATATAGTATCGAACATAGAAAAAACTGAAGATTCATCTTATTTTGTATATACTACCTAGCCCAACAATGTAGAGGTTAACGAAAAATGCGCTCACCAATTTGTCATCTTTTCTCAGCaattaa                                                                                          <  3:164098/149‑1 (MQ=255)
                                                       aaCATGTCGATAATTCATTACGCATATAGTATCGAACATAGAAAAAACTGAAGATTCATCTTATTTTGTATATACTACCTAGCCCAACAATGTAGAGGTTAACGaaaaa                                                                                                                               <  3:85171/109‑1 (MQ=255)
                                                       aaCATGTCGATAATTCATTACGCATATAGTATCGAACATAGAAAAAACTGAAGATTCATCTTATTTTGTATATACTACCTAGCCCAACAATGTAGAGGTTAACGaaaaa                                                                                                                               >  4:85171/1‑109 (MQ=255)
                                                       aaCATGTCGATAATTCATTACGCATATAGTATCGAACATAGAAAAAACTGAAGATTCATCTTATTTTGTATATACTACCTAGCCCAACAATGTAGAGGTTAACGaaaaa                                                                                                                               >  2:309635/1‑109 (MQ=255)
                                                       aaCATGTCGATAATTCATTACGCATATAGTATCGAAAATAGAAAAAACTGAAGATTCATCTTATTTTGTATATACTACCTAGCCCACCAATGTAGAGGTTAACGaaaaa                                                                                                                               <  1:309635/109‑1 (MQ=255)
                                                              ccaTAATTCATTACGCATATAGTATCGAACATAGAAAAAACTGAAGATTCATCTTATTTTGTATATACTACCTAGCCCAACAATGTAGAGGTTAACGAAAAATGCGCTCACCAATTTGTCATCTTTTCTCAGCAATTAATTCATCACCa                                                                                <  3:461676/147‑1 (MQ=255)
                                                                 tAATTCATTACGCATATAGTATCGAACATAGAAAAAACTGAAGATTCATCTTATTTTGTATATACTACCTAGCCCAACAATGTAGAGGTTAACGAAAAATGCGCTCACCAATTTGTCATCTTTTCTCAGCAATTAATTCATCACCAttt                                                                             >  1:12284/1‑149 (MQ=255)
                                                                   aTTCATTACGCATATAGTATCGAACATAGAAAAAACTGAAGATTCATCTTATTTTGTATATACTACCTAGCCCAACAATGTAGAGGTTAACGAAAAATGCGCTCACCAATTTGTCATCTTTTCTCAGCAATTAATTCATCACCATTTaa                                                                           <  4:515088/149‑1 (MQ=255)
                                                                     tCATTACGCATATAGTATCGAACATAGAAAAAACTGAAGATTCATCTTATTTTGTATATACTACCTAGCCCAACAATGTAGAGGTTAACGAAAAATGCGCTCACCAATTTGTCATCTTTTCTCAGCAATTAATTCATCACCATTTAAGa                                                                         <  3:134442/149‑1 (MQ=255)
                                                                       aTTACGCATATAGTATCGAACATAGAAAAAACTGAAGATTCATCTTATTTTGTATATACTACCTAGCCCAACAATGTAGAGGTTAACGAAAAATGCGCTCACCAATTTGTCATCTTTTCTCAGCAATTAATTCATCACCATTTAAGAtt                                                                       >  3:372242/1‑149 (MQ=255)
                                                                       aTTACGCATATAGTATCGAACATAGAAAAAACTGAAGATTCATCTTATTTTGTATATACTACCTAGCCCAACAATGTAGAGGTTAACGAAAAATGCGCTCACCAATTTGTCATCTTTTCTCAGCAATTAATTCATCACCATTTAAGAtt                                                                       >  4:383137/1‑149 (MQ=255)
                                                                        ttACGCATATAGTATCGAACATAGAAAAAACTGAAGATTCATCTTATTTTGTATATACTACCTAGCCCAACAATGTAGAGGTTAACGAAAAATGCGCTCACCAATTTGTCATCTTTTCTCAGCAATTAATTCATCACCATTTAAGATTg                                                                      <  1:248687/149‑1 (MQ=255)
                                                                         tACGCATATAGTATCGAACATAGAAAAAACTGAAGATTCATCTTATTTTGTATATACTACCTAGCCCAACAATGTAGAGGTTAACGAAAAATGCGCTCACCAATTTGTCATCTTTTCTCAGCAATTAATTCATCACCATTTAAGATTGc                                                                     >  1:225776/1‑149 (MQ=255)
                                                                              atatAGTATCGAACATAGAAAAAACTGAAGATTCATCTTATTTTGTATATACTACCTAGCCCAACTATGTAGAGGTTAACGAAAAATGCGCTCACCAATTTGTCATCTTTTCTCAGCAATTAATTCATCACCATTTAAGATTGCACCag                                                                >  4:326554/1‑149 (MQ=255)
                                                                               tataGTATCGAACATAGAAAAAACTGAAGATTCATCTTATTTTGTATATACTACCTAGCCCAACAATGTAGAGGTTAACGAAAAATGCGCTCACCAATTTGTCATCTTTTCTCAGCAATTAATTCATCACCAtt                                                                              <  2:514862/134‑1 (MQ=255)
                                                                               tataGTATCGAACATAGAAAAAACTGAAGATTCATCTTATTTTGTATATACTACCTAGCCCAACAATGTAGAGGTTAACGAAAAATGCGCTCACCAATTTGTCATCTTTTCTCAGCAATTAATTCATCACCAtt                                                                              >  1:338039/1‑134 (MQ=255)
                                                                               tataGTATCGAACATAGAAAAAACTGAAGATTCATCTTATTTTGTATATACTACCTAGCCCAACAATGTAGAGGTTAACGAAAAATGCGCTCACCAATTTGTCATCTTTTCTCAGCAATTAATTCATCACCAtt                                                                              >  1:514862/1‑134 (MQ=255)
                                                                               tataGTATCGAACATAGAAAAAACTGAAGATTCATCTTATTTTGTATATACTACCTAGCCCAACAATGTAGAGGTTAACGAAAAATGCGCTCACCAATTTGTCATCTTTTCTCAGCAATTAATTCATCACCAtt                                                                              <  2:338039/134‑1 (MQ=255)
                                                                                ataGTATCGAACATAGAAAAAACTGAAGATTCATCTTATTTTGTATATACTACCTAGCCCAACAATGTAGAGGTTAACGAAAAATGCGCTCACCAATTTGTCATCTTTTCTCAGCAATTAATTCATCACCATTTAAGATTGCACCagag                                                              >  2:229494/1‑149 (MQ=255)
                                                                                  aGTATCGAACATAGAAAAAACTGAAGATTCATCTTATTTTGTATATACTACCTAGCCCAACAATGTAGAGGTTAACGAAAAATGCGCTCACCAATTTGTCATCTTTTCTCAGCAATTAATTCATCACCATTTAAGATTGCACCAGAGaa                                                            <  2:562541/149‑1 (MQ=255)
                                                                                    tATCGAACATAGAAAAAACTGAAGATTCATCTTATTTTGTATATACTACCTAGCCCAACAATGTAGAGGTTAACGAAAAATGCGCTCACCAATTTGTCATCTTTTCTCAGCAATTAATTCATCACCATTTAAGATTGCACCAGAGAAgg                                                          >  1:2119/1‑149 (MQ=255)
                                                                                     aTCGAACATAGAAAAAACTGAAGATTCATCTTATTTTGTATATACTACCTAGCCCAACAATGTAGAGGTTAACGAAAAATGCGCTCACCAATTTGTCATCTTTTCTCAGCAATTAATTCATCACCATTTAAGATTGCACCAGAGAAGGa                                                         >  1:546136/1‑149 (MQ=255)
                                                                                         aaCATAGAAAAAACTGAAGATTCATCTTATTTTGTATATACTACCTAGCCCAACAATGTAGAGGTTAACGAAAAATGCGCTCACCAATTTGTCATCTTTTCTCAGCAATTAATTCATCACCATTTAAGATTGCACCAGAGAAGGAGCaa                                                     >  3:416299/1‑149 (MQ=255)
                                                                                         aaCATAGAAAAAACTGAAGATTCATCTTATTTTGTATATACTACCTAGCCCAACAATGTAGAGGTTAACGAAAAATGCGCTCACCAATTTGTCATCTTTTCTCAGCAATTAATTCATCACCATTTAAGATTGCACCAGAGAAGGAGCaa                                                     <  4:65297/149‑1 (MQ=255)
                                                                                         aaCATAGAAAAAACTGAAGATTCATCTTATTTTGTATATACTACCTAGCCCAACAATGTAGAGGTTAACGAAAAATGCGCTCACCAATTTGTCATCTTTTCTCAGCAATTAATTCATCACCATTTAAGATTGCACCAGAGAAGGAGCaa                                                     >  4:495111/1‑149 (MQ=255)
                                                                                                aaaaaaCTGAAGATTCATCTTATTTTGTATATACTACCTAGCCCAACAATGTAGAGGTTAACGAAAAATGCGCTCACCAATTTGTCATCTTTTCTCAGCAATTAATTCATCACCATTTAAGATTGCACCAGAGAAGGAGCAAGATCTTa                                              >  4:136600/1‑149 (MQ=255)
                                                                                                    aaCTGAAGATTCATCTTATTTTGTATATACTACCTAGCCCAACAATGTAGAGGTTAACGAAAAATGCGCTCACCAATTTGTCATCTTTTCTCAGCAATTAATTCATCACCATTTAAGATTGCACCAGAGAAGGa                                                         <  3:287641/134‑1 (MQ=255)
                                                                                                    aaCTGAAGATTCATCTTATTTTGTATATACTACCTAGCCCAACAATGTAGAGGTTAACGAAAAATGCGCTCACCAATTTGTCATCTTTTCTCAGCAATTAATTCATCACCATTTAAGATTGCACCAGAGAAGGa                                                         >  4:287641/1‑134 (MQ=255)
                                                                                                      cTGAAGATTCATCTTATTTTGTATATACTACCTAGCCCAACAATGTAGAGGTTAACGAAAAATGCGCTCACCAATTTGTCATCTTTTCTCAGCAATTAATTCATCACCATTTAAGATTGCACCAGAGAAGGAGCAAGATCTTAAAACGa                                        <  3:459867/149‑1 (MQ=255)
                                                                                                      cTGAAGATTCATCTTATTTTGTATATACTACCTAGCCCAACAATGTAGAGGTTAACGAAAAATGCGCTCACCAATTTGTCATCTTTTCTCAGCAATTAATTCATCACCATTTAAGATTGCACCAGAGAAGGAGCAAGATCTTAAAACGa                                        <  1:379054/149‑1 (MQ=255)
                                                                                                       tGAAGATTCATCTTATTTTGTATATACTACCTAGCCCAACAATGTAGAGGTTAACGAAAAATGCGCTCACCAATTTGTCATCTTTTCTCAGCAATTAATTCATCACCATTTAAGATTGCACCAGAGAAGGAGCAAGATCTTAAAACGAt                                       >  1:507037/1‑149 (MQ=255)
                                                                                                       tGAAGATTCATCTTATTTTGTATATACTACCTAGCCCAACAATGTAGAGGTTAACGAAAAATGCGCTCACCAATTTGTCATCTTTTCTCAGCAATTAATTCATCACCATTTAAGATTGCACCAGAGAAGGAGCAAGATCTTAAAACGAt                                       >  3:491116/1‑149 (MQ=255)
                                                                                                                 tCTTATTTTGTATATACTACCTAGCCCAACAATGTAGAGGTTAACGAAAAATGCGCTCACCAATTTGTCATCTTTTCTCAGCAATTAATTCATCACCATTTAAGATTGCACCAGAGAAGGAGCAAGATCTTAAAACGATAGTTgacgac                             >  2:332584/1‑149 (MQ=255)
                                                                                                                 tCTTATTTTGTATATACTACCTAGCCCAACAATGTAACGGTTAACGAAAAATGCGCTCACCAATGTGTCATCTTTTCTCAGAAATTAATTCATCACCATTTAAGTTTGAACCAGACAAGTAGCAAGATCTTAAAAAGATAGTTaacgac                             >  4:412344/1‑149 (MQ=255)
                                                                                                                   ttATTTTGTATATACTACCTAGCCCAACAATGTAGAGGTTAACGAAAAATGCGCTCACCAATTTGTCATCTTTTCTCAGCAATTAATTCATCACCATTTAAGATTGCACCAGAGAAGGAGCAAGATCTTAAAACGATAGTTGACGACaa                           >  2:576215/1‑149 (MQ=255)
                                                                                                                          gTATATACTACCTAGCCCAACAATGTAGAGGTTAACGAAAAATGCGCTCACCAATTTGTCATCTTTTCTCAGCAATTAATTCATCACCATTTAAGATTGCACCAGAGAAGGAGCAAGATCTTAAAACGATAGTTGACGACAAAAAAAtt                    <  4:73228/149‑1 (MQ=255)
                                                                                                                           tatataCTACCTAGCCCAACAATGTAGAGGTTAACGAAAAATGCGCTCACCAATTTGTCATCTTTTCTCAGCAATTAATTCATCACCATTTAAGATTGCACCAGAGAAGGAGCAAGATCTTAAAACGATAGTTGACGACAAAAAAATta                   >  1:6308/1‑149 (MQ=255)
                                                                                                                           tatataCTACCTAGCCCAACAATGTAGAGGTTAACGAAAAATGCGCTCACCAATTTGTCATCTTTTCTCAGCAATTAATTCATCACCATTTAAGATTGCACCAGAGAAGGAGCAAGATCTTAAAACGATAGTTGACGACAAAAAAATta                   >  4:31052/1‑149 (MQ=255)
                                                                                                                                actacCTAGCCCAACAATGTAGAGGTTAACGAAAAATGCGCTCACCAATTTGTCATCTTTTCTCAGCAATTAATTCATCACCATTTAAGATTGCACCAGAGAAGGAGCAAGATCTTAAAACGATAGTTGACGACAAAAAAATTATAAtt              >  4:91665/1‑149 (MQ=255)
                                                                                                                                 ctacCTAGCCCAACAATGTAGAGGTTAACGAAAAATGCGCTCACCAATTTGTCATCTTTTCTCAGCAATTAATTCATCACCATTTAAGATTGCACCAGAGAAGGAGCAAGATCTTAAAACGATAGTTGACGACAAAAAAATTATAAttt             <  2:535346/149‑1 (MQ=255)
                                                                                                                                  tacCTAGCCCAACAATGTAGAGGTTAACGAAAAATGCGCTCACCAATTTGTCATCTTTTCTCAGCAATTAATTCATCACCATTTAAGATTGCACCAGAGAAGGAGCAAGATCTTAAAACGATAGTTGACGACAAAAAAATTATAATTTc            >  1:33024/1‑149 (MQ=255)
                                                                                                                                          ccAACAATGTAGAGGTTAACGAAAAATGCGCTCACCAATTTGTCATCTTTTCTCAGCAATTAATTCATCACCATTTAAGATTGCACCAGAGAAGGAGCAAGATCTTAAAACGATAGTTGACGACAAAAAAATTATAATTTCAGTTgtga    <  3:91665/149‑1 (MQ=255)
                                                                                                                                           caacaaTGTAGAGGTTAACGAAAAATGCGCTCACCAATTTGTCATCTTTTCTCAGCAATTAATTCATCACCATTTAAGATTGCACCAGAGAAGGAGCAAGATCTTAAAACGATAGTTGACGACAAAAAAATTATAATTTCAGTTgtgag   >  4:17404/1‑149 (MQ=255)
                                                                                                                                           caacaaTGTAGAGGTTAACGAAAAATGCGCTCACCAATTTGTCATCTTTTCTCAGCAATTAATTCATCACCATTTAAGATTGCACCAGAGAAGGAGCAAGATCTTAAAACGATAGTTGACGACAAAAAAATTATAATTTCAGTTgtgag   >  4:141572/1‑149 (MQ=255)
                                                                                                                                            aacaaTGTAGAGGTTAACGAAAAATGCGCTCACCAATTTGTCATCTTTTCTCAGCAATTAATTCATCACCATTTAAGATTGCACCAGAGAAGGAGCAAGATCTTAAAACGATAGTTGACGACAAAAAAATTATAATTTCAGTTgtgagt  >  1:519301/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                 
GTCTTACGACATGCTTACCTTAATCGTTCCCAGTGTGCCTATAGCAGATATTCTAAACATGTCGATAATTCATTACGCATATAGTATCGAACATAGAAAAAACTGAAGATTCATCTTATTTTGTATATACTACCTAGCCCAACAATGTAGAGGTTAACGAAAAATGCGCTCACCAATTTGTCATCTTTTCTCAGCAATTAATTCATCACCATTTAAGATTGCACCAGAGAAGGAGCAAGATCTTAAAACGATAGTTGACGACAAAAAAATTATAATTTCAGTTGTGAGT  >  NC_000913/1198532‑1198820

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