Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,761,279 A→G intergenic (+106/+72) yfjJ → / ← abpB CP4‑57 prophage; protein YfjJ/CP4‑57 prophage; putative helicase YfjK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,761,2790AG98.5% 237.0 / ‑4.7 68intergenic (+106/+72)yfjJ/abpBCP4‑57 prophage; protein YfjJ/CP4‑57 prophage; putative helicase YfjK
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (30/37);  total (31/37)
Fisher's exact test for biased strand distribution p-value = 4.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

GAGACCGTAATTTCGGCTGCAGTCGCGGGTAATTTATTTTTTACAACAATGGCCTTTCTTATGGCTGGAAAGGCCTGATTTGTTTGTATTAAGATATCGGTAATCTTTTGATGGTTGGAAGCCCGTCCTGGTAGACGACCTAGCCTACCGAAGATGTGAGTTATTCAAAAATGTCGCTTTCCCGCGTTCCGTAGACAAACGTACCTAAATCAATGTCGTCTAATACCGAGCGCGCGCCTGATAAACATCTGATCAACGTACCCAATCGATC  >  NC_000913/2761142‑2761412
                                                                                                                                         |                                                                                                                                     
gagaCCGTAATTTCGGCTGCAGTCGCGGGTAATTTATTTTTTACAACAATGGCCTTTCTTATGGCTGGAAAGGCCTGATTTGTTTGTATTAAGATATCGGTAATCTTTTGATGGTTGGAAGCCCGTCCTGGTAGACGGcc                                                                                                                                     <  1:635829/140‑1 (MQ=255)
     cGTAATTTCGGCTGCAGTCGCGGGTAATTTATTTTTTACAACAATGGCCTTTCTTATGGCTGGAAAGGCCTGATTTGTTTGTATTAAGATATCGGTAATCTTTTGATGGTTGGAAGCCCGTCCTGGTAGACGGCCTa                                                                                                                                   >  1:385997/1‑137 (MQ=255)
     cGTAATTTCGGCTGCAGTCGCGGGTAATTTATTTTTTACAACAATGGCCTTTCTTATGGCTGGAAAGGCCTGATTTGTTTGTATTAAGATATCGGTAATCTTTTGATGGTTGGAAGCCCGTCCTGGTAGACGGCCTa                                                                                                                                   <  2:385997/137‑1 (MQ=255)
               gCTGCAGTCGCGGGTAATTTATTTTTTACAACAATGGCCTTTCTTATGGCTGGAAAGGCCTGATTTGTTTGTATTAAGATATCGGTAATCTTTTGATGGTTGGAAGCCCGTCCTGGTAGACGGCCTAGCCTACCGAAGAt                                                                                                                      >  1:159700/1‑140 (MQ=255)
                    aGTCGCGGGTAATTTATTTTTTACAACAATGGCCTTTCTTATGGCTGGAAAGGCCTGATTTGTTTGTATTAAGATATCGGTAATCTTTTGATGGTTGGAAGCCCGTCCTGGTAGACGGCCTAGCCTACCGAAGATGTGAg                                                                                                                 <  1:516642/140‑1 (MQ=255)
                    aGTCGCGGGTAATTTATTTTTTACAACAATGGCCTTTCTTATGGCTGGAAAGGCCTGATTTGTTTGTATTAAGATATAAGTAATCTTTTGATGGTTGGAAGCCCGTCCTGGTAGACGGCCTAGCCTACCGAAGATGTGAg                                                                                                                 <  2:773126/140‑1 (MQ=255)
                             tAATTTATTTTTTACAACAATGGCCTTTCTTATGGCTGGAAAGGCCTGATTTGTTTGTATTAAGATATCGGTAATCTTTTGATGGTTGGAAGCCCGTCCTGGTAGACGGCCTAGCCTACCGAAGATGTGAGTTATTCaaa                                                                                                        >  1:112960/1‑140 (MQ=255)
                                  tatttTTTACAACAATGGCCTTTCTTATGGCTGAAAAGGCCTGATTTGTTTGTATTAAGATATCGGTAATCTTTTGATGGTTGGAAGCCCGTCCTGGTAGACGGCCTAGCCTACCGAAGATGTGAGTTATTCAAAAATGt                                                                                                   <  1:764363/140‑1 (MQ=255)
                                      ttttACAACAATGGCCTTTCTTATGGCTGGAAAGGCCTGATTTGTTTGTATTAAGATATCGGTAATCTTTTGATGGTTGGAAGCCCGTCCTGGTAGACGGCCTAGCCTACCGAAGATGTGAGTTATTCAAAAATGTCGCt                                                                                               <  1:220453/140‑1 (MQ=255)
                                                aTGGCCTTTCTTATGGCTGGAAAGGCCTGATTTGTTTGTATTAAGATATCGGTAATCTTTTGATGGTTGGAAGCCCGTCCTGGTAGACGGCCTAGCCTACCGAAGATGTGAGTTATTCAAAAATGTCGCTTTCCCGCGtt                                                                                     <  2:112960/140‑1 (MQ=255)
                                                 tGGCCTTTCTTATGGCTGGAAAGGCCTGATTTGTTTGTATTAAGATATCGGTAATCTTTTGATGGTTGGAAGCCCGTCCTGGTAGACGGCCTAGCCTACCGAAGATGTGAGTTATTCAAAAATGTCGCTTTCCCGCGTTc                                                                                    >  2:7548/1‑140 (MQ=255)
                                                        tCTTATGGCTGGAAAGGCCTGATTTGTTTGTATTAAGATATCGGTAATCTTTTGATGGTTGGAAGCCCGTCCTGGTAGACGGCCTAGCCTACCGAAGATGTGAGTTATTCAAAAATGTCGCTTTCCCGCGTTCCGTa                                                                                <  2:334091/137‑1 (MQ=255)
                                                        tCTTATGGCTGGAAAGGCCTGATTTGTTTGTATTAAGATATCGGTAATCTTTTGATGGTTGGAAGCCCGTCCTGGTAGACGGCCTAGCCTACCGAAGATGTGAGTTATTCAAAAATGTCGCTTTCCCGCGTTCCGTa                                                                                >  1:334091/1‑137 (MQ=255)
                                                        tCTTATGGCTGGAAAGGCCTGATTTGTTTGTATTAAGATATCGGTAATCTTTTGATGGTTGGAAGCCCGTCCTGGTAGACGGCCTAGCCTACCGAAGATGTGAGTTATTCAAAAATGTCGCTTTCCCGCGTTCCGTAGAc                                                                             >  1:259116/1‑140 (MQ=255)
                                                            aTGGCTGGAAAGGCCTGATTTGTTTGTATTAAGATATCGGTAATCTTTTGATGGTTGGAAGCCCGTCCTGGTAGACGGCCTAGCCTACCGAAGATGTGAGTTATTCAAAAATGTCGCTTTCCCGCGTTCCGTAGACAAAc                                                                         <  2:331723/140‑1 (MQ=255)
                                                            aTGGCTGGAAAGGCCTGATTTGTTTGTATTAAGATATCGGTAATCTTTTGATGGTTGGAAGCCCGTCCTGGTAGACGGCCTAGCCTACCGAAGATGTGAGTTATTCAAAAATGTCGCTTTCCCGCGTTCCGTAGACAAAc                                                                         <  2:159700/140‑1 (MQ=255)
                                                             tGGCTGGAAAGGCCTGATTTGTTTGTATTAAGATATCGGTAATCTTTTGATGGTTGGAAGCCCGTCCTGGTAGACGGCCTAGCCTACCGAAGATGTGAGTTATTCAAAAATGTCGCTTTCCCGCGTTCCGTAGACAAAc                                                                         <  1:284185/139‑1 (MQ=255)
                                                             tGGCTGGAAAGGCCTGATTTGTTTGTATTAAGATATCGGTAATCTTTTGATGGTTGGAAGCCCGTCCTGGTAGACGGCCTAGCCTACCGAAGATGTGAGTTATTCAAAAATGTCGCTTTCCCGCGTTCCGTAGACAAAc                                                                         >  2:284185/1‑139 (MQ=255)
                                                              ggCTGGAAAGGCCTGATTTGTTTGTATTAAGATATCGGTAATCTTTTGATGGTTGGAAGCCCGTCCTGGTAGACGGCCTAGCCTACCGAAGATGTGAGTTATTCAAAAATGTCGCTTTCCCGCGTTCCGTAGACAAACGt                                                                       >  2:863772/1‑140 (MQ=255)
                                                                  ggAAAGGCCTGATTTGTTTGTATTAAGATATCGGTAATCTTTTGATGGTTGGAAGCCCGTCCTGGTAGACGGCCTAGCCTACCGAAGATGTGAGTTATTCAAAAATGTCGCTTTCCCGCGTTCCGTAGACAAACGTACCt                                                                   <  2:727738/140‑1 (MQ=255)
                                                                   gAAAGGCCCGATTTGTTTGTATTAAGATATCGGTAATCTTTTGATGGTTGGAAGCCCGTCCTGGTAGACGACCTAGCCTACCGAAGATGTGAGTTATTCAAAAATGTCGCTTTCCCGCGTTCCGTAGACAAACGTACCTa                                                                  >  2:77173/1‑140 (MQ=255)
                                                                         ccTGATTTGTTTGTATTAAGATATCGGTAATCTTTTGATGGTTGGAAGCCCGTCCTGGTAGACGGCCTAGCCTACCGAAGATGTGAGTTATTCAAAAATGTCGCTTTCCCGCGTTCCGTAGACAAACGTACCTAAATCaa                                                            >  2:27049/1‑140 (MQ=255)
                                                                          cTGATTTGTTTGTATTAAGATATCGGTAATCTTTTGATGGTTGGAAGCCCGTCCTGGTAGACGGCCTAGCCTACCGAAGATGTGAGTTATTCAAAAATGTCGCTTTCCCGCGTTCCGTAGACAAACGTACCTAAATCAAt                                                           >  2:610150/1‑140 (MQ=255)
                                                                           tGATTTGTTTGTATTAAGATATCGGTAATCTTTTGATGGTTGGAAGCCCGTCCTGGTAGACGGCCTAGCCTACCGAAGATGTGAGTTATTCAAAAATGTCGCTTTCCCGCGTTCCGTAGACAAACGTACCTAAATCAATg                                                          <  1:61312/140‑1 (MQ=255)
                                                                           tGATTTGTTTGTATTAAGATATCGGTAATCTTTTGATGGTTGGAAGCCCGTCCTGGTAGACGGCCTAGCCTACCGAAGATGTGAGTTATTCAAAAATGTCGCTTTCCCGCGTTCCGTAGACAAACGTACCTAAATCAATg                                                          <  1:863772/140‑1 (MQ=255)
                                                                           tGATTTGTTTGTATTAAGATATCGGTAATCTTTTGATGGTTGGAAGCCCGTCCTGGTAGACGGCCTAGCCTACCGAAGATGTGAGTTATTCAAAAATGTCGCTTTCCCGCGTTCCGTAGACAAACGTACCTAAATCAATg                                                          <  2:461302/140‑1 (MQ=255)
                                                                           tGATTTGTTTGTATTAAGATATCGGTAATCTTTTGATGGTTGGAAGCCCGTCCTGGTAGACGGCCTAGCCTACCGAAGATGTGAGTTATTCAAAAATGTCGCTTTCCCGCGTTCCGTAGACAAACGTACCTAAATCAATg                                                          <  2:522715/140‑1 (MQ=255)
                                                                                  tttgtATTAAGATATCGGTAATCTTTTGATGGTTGGAAGCCCGTCCTGGTAGACGGCCTAGCCTACCGAAGATGTGAGTTATTCAAAAATGTCGCTTTCCCGCGTTCCGTAGACAAACGTACCTAAATCAATGTCGTCTa                                                   >  1:583471/1‑140 (MQ=255)
                                                                                  tttgtATTAAGATATCGGTAATCTTTTGATGGTTGGAAGCCCGTCCTGGTAGACGGCCTAGCCTACCGAAGATGTGAGTTATTCAAAAATGTCGCTTTCCCGCGTTCCGTAGACAAACGTACCTAAATCAATGTCGTCTa                                                   >  2:123889/1‑140 (MQ=255)
                                                                                  tttgtATTAAGATATCGGTAATCTTTTGATGGTTGGAAGCCCGTCCTGGTAGACGGCCTAGCCTACCGAAGATGTGAGTTATTCAAAAATGTCGCTTTCCCGCGTTCCGTAGACAAACGTACCTAAATCAATGTCGTCTa                                                   >  1:806314/1‑140 (MQ=255)
                                                                                   ttgtATTAAGATATCGGTAATCTTTTGATGGTTGGAAGCCCGTCCTGGTAGACGGCCTAGCCTACCGAAGATGTGAGTTATTCAAAAATGTCGCTTTCCCGCGTTCCGTAGACAAACGTACCTAAATCAATg                                                          >  2:340689/1‑132 (MQ=255)
                                                                                   ttgtATTAAGATATCGGTAATCTTTTGATGGTTGGAAGCCCGTCCTGGTAGACGGCCTAGCCTACCGAAGATGTGAGTTATTCAAAAATGTCGCTTTCCCGCGTTCCGTAGACAAACGTACCTAAATCAATg                                                          <  1:146552/132‑1 (MQ=255)
                                                                                   ttgtATTAAGATATCGGTAATCTTTTGATGGTTGGAAGCCCGTCCTGGTAGACGGCCTAGCCTACCGAAGATGTGAGTTATTCAAAAATGTCGCTTTCCCGCGTTCCGTAGACAAACGTACCTAAATCAATg                                                          <  1:340689/132‑1 (MQ=255)
                                                                                   ttgtATTAAGATATCGGTAATCTTTTGATGGTTGGAAGCCCGTCCTGGTAGACGGCCTAGCCTACCGAAGATGTGAGTTATTCAAAAATGTCGCTTTCCCGCGTTCCGTAGACAAACGTACCTAAATCAATg                                                          >  2:146552/1‑132 (MQ=255)
                                                                                   ttgtATTAAGATATCGGTAATCTTTTGATGGTTGGAAGCCCGTCCTGGTAGACGGCCTAGCCTACCGAAGATGTGAGTTATTCAAAAATGTCGCTTTCCCGCGTTCAGTGGACAAACGTACCTAAATCAATg                                                          >  2:655807/1‑132 (MQ=255)
                                                                                   ttgtATTAAGATATCGGTAATCTTTTGATGGTTGGAAGCCCGTCCTGGTAGACGGCCTAGCCTACCGAAGATGTGAGTTATTCAAAAATGTCGCTTTCCCGCGTTCAGTGGACAAACGTACCTAAATCAATg                                                          <  1:655807/132‑1 (MQ=255)
                                                                                               atCGGTAATCTTTTGATGGTTGGAAGCCCGTCCTGGTAGACGGCCTAGCCTACCGAAGATGTGAGTTATTCAAAAATGTCGCTTTCCCGCGTTCCGTAGACAAACGTACCTAAATCAATGTCGTCTAATACCGAgcgcgc                                      <  2:228281/140‑1 (MQ=255)
                                                                                                tCGGTAATCTTTTGATGGTTGGAAGCCCGTCCTGGTAGACGGCCTAGCCTACCGAAGATGTGAGTTATTCAAAAATGTCGCTTTCCCGCGTTCCGTAGACAAACGTACCTAAATCAATGTCGTCTAATACCGAgcgcgcg                                     >  2:668978/1‑140 (MQ=255)
                                                                                                tCGGTAATCTTTTGATGGTTGGAAGCCCGTCCTGGTAGACGGCCTAGCCTACCGAAGATGTGAGTTATTCAAAAATGTCGCTTTCCCGCGTTCCGTAGACAAACGTACCTAAATCAATGTCGTCTAATACCGAgcgcgcg                                     >  2:289119/1‑140 (MQ=255)
                                                                                                tCGGTAATCTTTTGATGGTTGGAAGCCCGTCCTGGTAGACGGCCTAGCCTACCGAAGATGTGAGTTATTCAAAAATGTCGCTTTCCCGCGTTCCGTAGACAAACGTACCTAAATCAATGTCGTCTAATACCGAgcgcgcg                                     >  1:702963/1‑140 (MQ=255)
                                                                                                 cGGTAATCTTTTGATGGTTGGAAGCCCGTCCTGGTAGACGGCCTAGCCTACCGAAGATGTGAGTTATTCAAAAATGTCGCTTTCCCGCGTTCCGTAGACAAACGTACCTAAATCAATGTCGTCTAATACCGAgcgcgcgc                                    >  1:480776/1‑140 (MQ=255)
                                                                                                  ggTAATCTTTTGATGGTTGGAAGCCCGTCCTGGTAGACGGCCTAGCCTACCGAAGATGTGAGTTATTCAAAAATGTCGCTTTCCCGCGTTCCGTAGACAAACGTACCTAAATCAATGTCGTCTAATACCGAGCGCGCGcc                                   <  2:651764/140‑1 (MQ=255)
                                                                                                  ggTAATCTTTTGATGGTTGGAAGCCCGTCCTGGTAGACGGCCTAGCCTACCGAAGATGTGAGTTATTCAAAAATGTCGCTTTCCCGCGTTCCGTAGACAAACGTACCTAAATCAATGTCGTCTAATACCGAGCGCGCGcc                                   <  1:289119/140‑1 (MQ=255)
                                                                                                   gTAATCTTTTGATGGTTGGAAGCCCGTCCTGGTAGACGGCCTAGCCTACCGAAGATGTGAGTTATTCAAAAATGTCGCTTTCCCGCGTTCCGTAGACAAACGTACCTAAATCAATGTCGTCTAATACCGAGCGCGCGCCt                                  >  2:391316/1‑140 (MQ=255)
                                                                                                     aaTCTTTTGATGGTTGGAAGCCCGTCCTGGTAGACGGCCTAGCCTACCGAAGATGTGAGTTATTCAAAAATGTCGCTTTCCCGCGTTCCGTAGACAAACGTACCTAAATCAATGTCGTCTAATACCGAgcgcg                                       >  2:782978/1‑133 (MQ=255)
                                                                                                     aaTCTTTTGATGGTTGGAAGCCCGTCCTGGTAGACGGCCTAGCCTACCGAAGATGTGAGTTATTCAAAAATGTCGCTTTCCCGCGTTCCGTAGACAAACGTACCTAAATCAATGTCGTCTAATACCGAgcgcg                                       <  1:782978/133‑1 (MQ=255)
                                                                                                       tCTTTTGATGGTTGGAAGCCCGTCCTGGTAGACGGCCTAGCCTACCGAAGATGTGAGTTATTCAAAAATGTCGCTTTCCCGCGTTCCGTAGACAAACGTACCTAAATCAATGTCGTCTAATACCGAGCGCGCGCCTGATa                              <  2:806314/140‑1 (MQ=255)
                                                                                                       tCTTTTGATGGTTGGAAGCCCGTCCTGGTAGACGGCCTAGCCTACCGAAGATGTGAGTTATTCAAAAATGTCGCTTTCCCGCGTTCCGTAGACAAACGTACCTAAATCAATGTCGTCTAATACCGAGCGCGCGCCTGATa                              <  1:391316/140‑1 (MQ=255)
                                                                                                           ttGATGGTTGGAAGCCCGTCCTGGTAGACGGCCTAGCCTACCGAAGATGTGAGTTATTCAAAAATGTCGCTTTCCCGCGTTCCGTAGACAAACGTACCTAAATCAATGTCGTCTAATACCGAGCGCGCGCCTGATAAACa                          <  1:27049/140‑1 (MQ=255)
                                                                                                           ttGATGGTTGGAAGCCCGTCCTGGTAGACGGCCTAGCCTACCGAAGATGTGAGTTATTCAAAAATGTCGCTTTCCCGCGTTCCGTAGACAAACGTACCTAAATCAATGTCGTCTAATACCGAGCGCGCGCCTGATAAACa                          <  2:583471/140‑1 (MQ=255)
                                                                                                            tGATGGTTGGAAGCCCGTCCTGGTAGACGGCCTAGCCTACCGAAGATGTGAGTTATTCAAAAATGTCGCTTTCCCGCGTTCCGTAGACAAACGTACCTAAATCAATGTCGTCTAATACCGAGCGCGCGCCTGATAAACAt                         >  2:217303/1‑140 (MQ=255)
                                                                                                            tGATGGTTGGAAGCCCGTCCTGGTAGACGGCCTAGCCTACCGAAGATGTGAGTTATTCAAAAATGTCGCTTTCCCGCGTTCCGTAGACAAACGTACCTAAATCAATGTCGTCTAATACCGAGCGCGCGCCTGATAAACAt                         <  2:636381/140‑1 (MQ=255)
                                                                                                            tGATGGTTGGAAGCCCGTCCTGGTAGACGGCCTAGCCTACCGAAGATGTGAGTTATTCAAAAATGTCGCTTTCCCGCGTTCCGTAGACAAACGTACCTAAATCAATGTCGTCTAATACCGAGCGCGCGCCTGATAAACAt                         >  1:382453/1‑140 (MQ=255)
                                                                                                            tGATGGTTGGAAGCCCGTCCTGGTAGACGGCCTAGCCTACCGAAGATGTGAGTTATTCAAAAATGTCGCTTTCCCGCGTTCCGTAGACAAACGTACCTAAATCAATGTCGTCTAATACCGAGCGCGCGCCTGATAAACAt                         <  1:668978/140‑1 (MQ=255)
                                                                                                               tGGTTGGAAGCCCGTCCTGGTAGACGGCCTAGCCTACCGAAGATGTGAGTTATTCAAAAATGTCGCTTTCCCGCGTTCCGTAGACAAACGTACCTAAATCAATGTCGTCTAATACCGAGCGCGCGCCTGATAAACATCTg                      <  2:388054/140‑1 (MQ=255)
                                                                                                                ggTTGGAAGCCCGTCCTGGTAGACGGCCTAGCCTACCGAAGATGTGAGTTATTCAAAAATGTCGCTTTCCCGCGTTCCGTAGACAAACGTACCTAAATCAATGTCGTCTAATACCGAGCGCGCGCCTGATAAACATCTGa                     <  2:306611/140‑1 (MQ=255)
                                                                                                                   tGGAAGCCCGTCCTGGTAGACGGCCTAGCCTACCGAAGATGTGAGTTATTCAAAAATGTCGCTTTCCCGCGTTCCGTAGACAAACGTACCTAAATCAATGTCGTCTAATACCGAGCGCGCGCCTGATAAACATCTGATCa                  <  1:7548/140‑1 (MQ=255)
                                                                                                                   tGGAAGCCCGTCCTGGTAGACGGCCTAGCCTACCGAAGATGTGAGTTATTCAAAAATGTCGCTTTCCCGCGTTCCGTAGACAAACGTACCTAAATCAATGTCGTCTAATACCGAGCGCGCGCCTGATAAACATCTGATCa                  >  2:577757/1‑140 (MQ=255)
                                                                                                                   tGGAAGCCCGTCCTGGTAGACGGCCTAGCCTACCGAAGATGTGAGTTATTCAAAAATGTCGCTTTCCCGCGTTCCGTAGACAAACGTACCTAAATCAATGTCGTCTAATACCGAGCGCGCGCCTGATAAACATCTGATCa                  >  1:671329/1‑140 (MQ=255)
                                                                                                                     gAAGCCCGTCCTGGTAGACGGCCTAGCCTACCGAAGATGTGAGTTATTCAAAAATGTCGCTTTCCCGCGTTCCGTAGACAAACGTACCTAAATCAATGTCGTCTAATACCGAGCGCGCGCCTGATAAACATCTGATCAAc                >  1:30611/1‑140 (MQ=255)
                                                                                                                           cGTCCTGGTAGACGGCCTAGCCTACCGAAGATGTGAGTTATTCAAAAATGTCGCTTTCCCGCGTTCCGTAGACAAACGTACCTAAATCAATGTCGTCTAATACCGAGCGCGCGCCTGATAAACATCTGATCAACGTAccc          >  1:251650/1‑140 (MQ=255)
                                                                                                                           cGTCCTGGTAGACGGCCTAGCCTACCGAAGATGTGAGTTATTCAAAAATGTCGCTTTCCCGCGTTCCGTAGACAAACGTACCTAAATCAATGTCGTCTAATACCGAGCGCGCGCCTGATAAACATCTGATCAACGTAccc          <  1:215223/140‑1 (MQ=255)
                                                                                                                           cGTCCTGGTAGACGGCCTAGCCTACCGAAGATGTGAGTTATTCAAAAATGTCGCTTTCCCGCGTTCCGTAGACAAACGTACCTAAATCAATGTCGTCTAATACCGAGCGCGCGCCTGATAAACATCTGATCAACGTAccc          >  1:586216/1‑140 (MQ=255)
                                                                                                                            gTCCTGGTAGACGGCCTAGCCTACCGAAGATGTGAGTTATTCAAAAATGTCGCTTTCCCGCGTTCCGTAGACAAACGTACCTAAATCAATGTCGTCTAATACCGAGCGCGCGCCTGATAAACATCTGa                     <  2:103100/128‑1 (MQ=255)
                                                                                                                            gTCCTGGTAGACGGCCTAGCCTACCGAAGATGTGAGTTATTCAAAAATGTCGCTTTCCCGCGTTCCGTAGACAAACGTACCTAAATCAATGTCGTCTAATACCGAGCGCGCGCCTGATAAACATCTGa                     >  1:103100/1‑128 (MQ=255)
                                                                                                                             tCCTGGTAGACGGCCTAGCCTACCGAAGATGTGAGTTATTCAAAAATGTCGCTTTCCCGCGTTCCGTAGACAAACGTACCTAAATCAATGTCGTCTAATACCGAGCGCGCGCCTGATAAACATCTGATCAACGTACCCaa        <  1:610150/140‑1 (MQ=255)
                                                                                                                              ccTGGTAGACGGCCTAGCCTACCGAAGATGTGAGTTATTCAAAAATGTCGCTTTCCCGCGTTCCGTAGACAAACGTACCTAAATCAATGTCGTCTAATACCGAGCGCGCGCCTGATAAACATCTGATCAACGTACCCAAt       <  2:30611/140‑1 (MQ=255)
                                                                                                                                   tAGACGGCCTAGCCTACCGAAGATGTGAGTTATTCAAAAATGTCGCTTTCCCGCGTTCCGTAGACAAACGTACCTAAATCAATGTCGTCTAATACCGAGCGCGCGCCTGATAAACATCTGATCAACGTACCCAATCGATc  <  1:241248/140‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                     
GAGACCGTAATTTCGGCTGCAGTCGCGGGTAATTTATTTTTTACAACAATGGCCTTTCTTATGGCTGGAAAGGCCTGATTTGTTTGTATTAAGATATCGGTAATCTTTTGATGGTTGGAAGCCCGTCCTGGTAGACGACCTAGCCTACCGAAGATGTGAGTTATTCAAAAATGTCGCTTTCCCGCGTTCCGTAGACAAACGTACCTAAATCAATGTCGTCTAATACCGAGCGCGCGCCTGATAAACATCTGATCAACGTACCCAATCGATC  >  NC_000913/2761142‑2761412

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