Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 294,752 (T)7→6 coding (48/855 nt) yagM ← CP4‑6 prophage; protein YagM

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913294,7460T.98.4% 248.7 / ‑2.8 61coding (54/855 nt)yagMCP4‑6 prophage; protein YagM
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base . (29/31);  minor base A (0/1);  total (29/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAAGAGTGAGATTTATTTGTTAGTCTGGCATAATGATCAATTACTCCTAATGAAGATTGATTTTTGAACGTGACTGAGTGTTTTTTTATGGGAGTAGTCTGTACTACCTGTGACTCATTGTTTTTGCCTTGACAAGGTTTTTTTCCTGGTAGAACGCTGCTCATCTCAAAATTAGTAACACTATTAGACATTTTTAAGGAAACCTCTTTGTTATCTCAATATAAAACTACTATCTGAAATTTCTCAGTGTTGCCCCATTTTTCTCGGGGCA  >  NC_000913/294609‑294879
                                                                                                                                         |                                                                                                                                     
gAAGAGTGAGATTTATTTGTTAGTCTGGCATAATGATCAATTACTCCTAATGAAGATTGATTTTTGAACGTGACTGAGTGTTTTTTTATGGGAGTAGTCTGTACTACCTGTGACTCATTGTTTTTGCCTTGACAAGGttt                                                                                                                                     <  2:87839/140‑1 (MQ=255)
  agagTGAGATTTATTTGTTAGTCTGGCATAATGATCAATTACTCCTAATGAAGATTGATTCTTGAACGTGACTGAGTGTTTTTTTATGGGAGTAGTCTGTACTACCTGTGACTCATTGTTTTTGCCTTGACAAGGttttt                                                                                                                                   <  2:502942/140‑1 (MQ=255)
     gTGAGATTTATTTGTTAGTCTGGCATAATGATCAATTACTCCTAATGAAGATTGATTTTTGAACGTGACTGAGTGTTTTTTTATGGGAGTAGTCTGTACTACCTGTGACTCATTGTTTTTGCCTTGACAAGGttttttcc                                                                                                                                >  1:534742/1‑138 (MQ=255)
        agaTTTATTTGTTAGTCTGGCATAATGATCAATTACTCCTAATGAAGATTGATTTTTGAACGTGACTGAGTGTTTTTTTATGGGAGTAGTCTGTACTACCTGTGACTCATTGTTTTTGCCTTGACAAGGTTTTTTCCTgg                                                                                                                            >  2:51210/1‑140 (MQ=255)
          atttatttGTTAGTCTGGCATAATGATCAATTACTCCTAATGAAGATTGATTTTTGAACGTGACTGAGTGTTTTTTTATGGGAGTAGTCTGTACTACCTGTGACTCATTGTTTTTGCCTTGACAAGGTTTTTTCCTGGTa                                                                                                                          <  1:51210/140‑1 (MQ=255)
          atttatttGTTAGTCTGGCATAATGATCAATTACTCCTAATGAAGATTGATTTTTGAACGTGACTGAGTGTTTTTTTATGGGAGTAGTCTGTACTACCTGTGACTCATTGTTTTTGCCTTGACAAGGTTTTTTCCTGGTa                                                                                                                          <  1:186336/140‑1 (MQ=255)
          atttatttGTTAGTCTGGCATAATGATCAATTACTCCTAATGAAGATTGATTTTTGAACGTGACTGAGTGTTTTTTTATGGGAGTAGTCTGTACTACCTGTGACTCATTGTTTTTGCCTTGACAAGGTTTTTTCCTGGTa                                                                                                                          <  2:534742/140‑1 (MQ=255)
           tttatttGTTAGTCTGGCATAATGATCAATTACTCCTAATGAAGATTGATTTTTGAACGTGACTGAGTGTTTTTTTATGGGAGTAGTCTGTACTACCTGTGACTCATTGTTTTTGCCTTGACAAGGTTTTTTCCTGGTAg                                                                                                                         >  2:828234/1‑140 (MQ=255)
             tatttGTTAGTCTGGCATAATGATCAATTACTCCTAATGAAGATTGATTTTTGAACGTGACTGAGTGTTTTTTTATGGGAGTAGTCTGTACTACCTGTGACTCATTGTTTTTGCCTTGACAAGGTTTTTTCCTGGTAGaa                                                                                                                       <  2:474319/140‑1 (MQ=255)
                   ttAGTCTGGCATAATGATCAATTGCTCCTAATGAAGATTGATTTTTGAACGTGACTGAGTGTTTTTTTATGGGAGTAGTCTGTACTACCTGTGACTCATTGTTTTTGCCTTGACAAGGTTTTTTCCTGGTAGAACgctgc                                                                                                                 <  1:183240/140‑1 (MQ=255)
                     aGTCTGGCATAATGATCAATTACTCCTAATGAAGATTGATTTTTGAACGTGACTGAGTGTTTTTTTATGGGAGTAGTCTGTACTACCTGTGACTCATTGTTTTTGCCTTGACAAGGTTTTTTCCTGGTAGAACGCTGCTc                                                                                                               <  1:828234/140‑1 (MQ=255)
                       tCTGGCATAATGATCAATTACTCCTAATGAAGATTGATTTTTGAACGTGACTGAGTGTTTTTTTATGGGAGTAGTCTGTACTACCTGTGACACATTGTTTATGCATTGACAAGGTTTATTCCAGGTAGAACGCTGCTCAt                                                                                                             >  2:264256/1‑140 (MQ=255)
                          ggCATAATGATCAATTACTCCTAATGAAGATTGATTTTTGAACGTGACTGAGTGTTTTTTTATGGGAGTAGTCTGTACTACCTGTGACTCATTGTTTTTGCCTTGACAAGGTTTTTTCCTGGTAGAACGCTGCTCAtctc                                                                                                          >  1:49633/1‑140 (MQ=255)
                           gCATAATGATCAATTACTCCTAATGAAGATTGATTTTTGAACGTGACTGAGTGTTTTTTTATGGGAGTAGTCTGTACTACCTGTGACTCATTGTTTTTGCCTTGACAAGGCTTTCTCGTGGTAGAACGCTGCTCATCTCa                                                                                                         >  1:273557/1‑140 (MQ=255)
                             aTAATGATCAATTACTCGTAATCAAGATTGATTTTTGAACGTGACTGAATGTTTTTTTATGGGAGTAGTCTGTACTGCCTGTGACTCATTGTTTTTGCCTTGACAAGGTTTTTTCCTGGTAGAACGCTGCTCATCTCaaa                                                                                                       <  2:357520/140‑1 (MQ=255)
                             aTAATGATCAATTACTCCTAATGAAGATTGATTTTTGAACGTGACTGAGTGTTTTTTTATGGGAGTAGTCTGTACTACCTGTGACTCATTGTTTTTGCCTTGACAAGGTTTTTTCCTGGTAGAACGCTGCTCATCTCaaa                                                                                                       <  1:163720/140‑1 (MQ=255)
                                aTGATCAATTACTCCTAATGAAGATTGATTTTTGAACGTGACTGAGTGTTTTTTTATGGGAGTAGTCTGTACTACCTGTGACTCATTGTTTTTGCCTTGACAAGGTTTTTTCCTGGTAGAACGCTGCTCATCTCAAAAtt                                                                                                    >  2:664539/1‑140 (MQ=255)
                                   acagaTTACTCCTAATGCAGATAGATTTTTGAACGCGACCGAGTGTTTTTTTATGGGAGTAGTCTGTACTACCTGTGACTCATTGTTTTTGCCTTGACAAGGTTTTTTCCTGGTAGAACGCTGCTCATCTCAAAATTAGt                                                                                                 <  2:270377/136‑1 (MQ=255)
                                     cAATTACTCCTAATGAAGATTGTTTCTTGCACGTGACTGAGTGTTTTTTTATGGGAGTAGTCTGTACTACCTGTGACTCATTGTTTTTGCCTTGACAAGGTTTTTTCCTGGTAGAACGCTGCTCATCTCAAAATTAGTaa                                                                                               <  1:264256/140‑1 (MQ=255)
                                     cAATTACTCCTAATGAAGATTGATTTTTGAACGTGACTGAGTGTTTTTTTATGGGAGTAGTCTGTACTACCTGTGACTCATTGTTTTTGCCTTGACAAGGTTTTTTCCTGGTAGAACGCTGCTCATCTCAAAATTAGTaa                                                                                               <  2:17978/140‑1 (MQ=255)
                                     cAATTACTCCTAATGAAGATTGATTTTTGAACGTGACTGAGTGTTTTTTTATGGGAGTAGTCTGTACTACCTGTGACTCATTGTTTTTGCCTTGACAAGGTTTTTTCCTGGTAGAACGCTGCTCATCTCAAAATTAGTaa                                                                                               <  1:797277/140‑1 (MQ=255)
                                     cAATTACTCCTAATGAAGATTGATTTTTGAACGTGACTGAGTGTTTTTTTATGGGAGTAGTCTGTACTACCTGTGACTCATTGTTTTTGCCTTGACAAGGTTTTTTCCTGGTAGAACGCTGCTCATCTCAAAATTAGTaa                                                                                               <  1:816792/140‑1 (MQ=255)
                                     cAATTACTCCTAATGAAGATTGATTTTTGAACGTGACTGAGTGTTTTTTTATGGGAGTAGTCTGTACTACCTGTGACTCATTGTTTTTGCCTTGACAAGGTTTTTTCCTGGTAGAACGCTGCTCATCTCAAAATTAGTaa                                                                                               <  1:464890/140‑1 (MQ=255)
                                     cAATTACTCCTAATGAAGATTGATTTTTGAACGTGACTGAGTGTTTTTTTATGGGAGTAGTCTGTACTACCTGTGACTCATTGTTTTTGCCTTGACAAGGTTTTTTCCTGGTAGAACGCTGCTCATCTCAAAATTAGTaa                                                                                               <  1:497716/140‑1 (MQ=255)
                                     cAATTACTCCTAATGAAGATTGATTTTTGAACGTGACTGAGTGTTTTTTTATGGGAGTAGTCTGTACTACCTGTGACTCATTGTTTTTGCCTTGACAAGGTTATTTCCTGTTAGAACGCTGCTCGTCTCAAAATTAGTaa                                                                                               >  1:357520/1‑140 (MQ=255)
                                     cAATTACTCCTAATGAAGATTGATTTTTGAACGTGACTGAGTGTTTTTTTATGGGAGTAGTCTGTACTACCTGTGACTCATTGTTTTTGCCCTGACAAGGATTTATCCGGGTAGAACGCTGCTCATCTCAAAATTAGTaa                                                                                               >  2:347805/1‑140 (MQ=255)
                                     cAATTACCCCTAATGAAGATTGATTTTTGAACGTGACTGAGTGTTTTTTTATGGGAGTAGTCTGTACTACCTGTGACTCATTGTTTTTGCCTTGACAAGGTTTTTTCCTGGTAGAACGCTGCTCATCTCAAAATTAGTaa                                                                                               <  1:827604/140‑1 (MQ=255)
                                       aTTACTCCTAATGAAGATTGATTTTTGAACGTGACTGAGTGTTTTTTTATGGGAGTAGTCTGTACTACCTGTGACTCATTGTTTTTGCCTTGACAAGGTTTTTTCCTGGTAGAACGCTGCTCATCTCAAAATTAGTAaca                                                                                             <  2:62078/140‑1 (MQ=255)
                                           cTCCTAATGAAGATTGATTTTTGAACGTGACTGAGTGTTTTTTTATGGGAGTAGTCTGTACTACCTGTGACTCATTGTTTTTGCCTTGACAAGGTTTTTTCCTGGTAGAACGCTGCTCATCTCAAAATTAGTAACACTAt                                                                                         >  1:319753/1‑140 (MQ=255)
                                           cTCCTAATGAAGATTGATTTTTGAACGTGACTGAGTGTTTTTT‑‑TGGGAGAACTGTCTACTACCTGTGAGTCATTGTTTTTTTCTCGAGAAGCAGTTTTCCCGTTATAAACCTGCGCAACTCAACATTATAACCACTATTa                                                                                       >  2:411759/1‑140 (MQ=255)
                                             ccTAATGAAGATTGATTTTTGAACGTGACTGAGTGTTTTTTTATGGGAGTAGTCTGTACTACCTGTGACTCATTGTTTTTGCCTTGACAAGGTTTTTTCCTGGTAGAACGCTGCTCATCTCAAAATTAGTAACACTATTa                                                                                       >  1:516064/1‑140 (MQ=255)
                                              cTAATGAAGATTGATTTTTGAACGTGACTGAGTGTTTTTTTATGGGAGTAGTCTGTACTACCTGTGACTCATTGTTTTTGCCTTGACAAGGCTTTTTCCTTGTAAAACACTGCTCATCTCAAGATTAGTAACACTATTAg                                                                                      >  2:364598/1‑140 (MQ=255)
                                               tAATGAAGATTGATTCTTGAACCTGTCTGCGTGTTTTCTTATGGGAGTAGTCTGTACTACCTGTGACTCATTGTTTTTGCCTTGACAAGGTTTTTTCCTGGTAGAACGCTGCTCATCTCAAAATTAGTAACACTATTAGa                                                                                     <  1:253412/140‑1 (MQ=255)
                                                aaTGAAGATTGATTTTTGAACGTGACTGAGTGTTTTTTTATGGGAGTAGTTTGTACTACCTGTGACTCATTGTTTTTGCCTTGACAAGGTTTTTTCCTGGTAGAACGCTGCTCATCTCAAAATTAGTAACACTATTAGAc                                                                                    >  2:437705/1‑140 (MQ=255)
                                                 aTGAAGATTGATTTTTGAACGTGACTGAGTGTTTTTTTATGGGAGTAGTCTGTACTACCTGTGACTCATTGTTTTTGCCTTGACAAGGTTTTTTCCTGGTAGAACGCTGCTCATCTCAAAATTAGTAACACTATTAGACa                                                                                   >  2:135290/1‑140 (MQ=255)
                                                  tGAAGATTGATTTTTGAACGTGACTGAGTGTTTTTTTATGGGAGTAGTCTGTACTACCTGTGACTCATTGTTTTTGCCTTGACAAGGTTTTTTCCTGGTAGACCGCTGCTAATCACATAATTAGGAACACTATTAgctgt                                                                                  >  1:270377/1‑136 (MQ=255)
                                                  tGAAGATTGATTTTTGAACGTGACTGAGTGTTTTTTTATGGGAGTAGTCTGTACTACCTGTGACTCATTGTTTTTGCCTTGACAAGGTTTTTTCCTGGTAGAACGCTGCTCATCTCAAAATTAGTAACACTATTAGACAt                                                                                  >  2:847169/1‑140 (MQ=255)
                                                           attTTTGAACGTGACTGAGTGTTTTTTTATGGGAGTAGTCTGTACTACCTGTGACTCATTGTTTTTGCCTTGACAAGGTTTTTTCCTGGTAGAACGCTGCTCATCTCAAAATTAGTAACACTATTAGACATTTTTAAGGa                                                                         >  2:677689/1‑140 (MQ=255)
                                                           attTTTGAACGTGACTGAGTGTCTTTTTATGGGAGTAGTCTGTACTACCTGTGACTCATTGTTTTTGCCTTGACAAGGTTTTTTCCTGGTAGAACGCTGCTCATCTCAAAATTAGTAACACTATTAGACATTTTTAAGGa                                                                         <  1:347805/140‑1 (MQ=255)
                                                           attTTTGAACGTGACAGAGTGTCTTGTTACGGGAGTATTCTGTACTACCAGTGACTCATTGTTTTTGCCTTGACAAGGTTTTTTCCTGGTAGAACGCTGCTCATCTCAAAATTAGTAACACTATTAGACATTTTTAAGGa                                                                         <  1:264499/140‑1 (MQ=255)
                                                            tttttGAACGTGACTGAGTGTTTTTTTATGGGAGTAGTCTGTACTACCTGTGACTCATTGTTTTTGCCTTGACAAGGTTTTTTCCTGGTAGAACGCTGCTCATCTCAAAATTAGTAACACTATTAGACATTTTTAAGGaa                                                                        <  1:677689/140‑1 (MQ=255)
                                                               ttGAACGTGACTGAGTGCTTTTTTATGGTAGCAGTCTGTCCTACCTGTTCCTCTTTGTTTTTGCCTTGACAAGGTTTTTTCCTGGTAGAACGCTGCTCATCTCAAAATTAGTAACACTATTAGACATTTTTAAGGAAAcc                                                                     <  2:234857/140‑1 (MQ=255)
                                                                tGAACGTGACTGAGTGTGTTTTTATGGGAGTCGTCAGTACTACCTGTGACTCATTGTTTTTGCCTTGACAAGGTTTTTTCCTGGTAGAACGCTGCTCATCTCAAAATTAGTAACACTATTAGACATTTTTAAGGAAACct                                                                    <  2:273557/140‑1 (MQ=255)
                                                                         cTGAGTGTTTTTTTATGGGAGTGGTCTGTACCACCGGTGGCTCATTGTTTTTGCCTTGACAAGGTTTTTTCCTGGTAGAACGCTGCTCATCTCAAAATTAGTAACACTATTAGACATTTTTAAGGAAACCTCTTTGTTAt                                                           <  2:319753/140‑1 (MQ=255)
                                                                           gAGTGTTTTTTTATGGGAGTAGTCTGTACTACCTGTGACTCATTGTTTTTGCCTTGACAAGGTTTTTTCCTGGTAGAACGCTGCTCATCTCAAAATTAGTAACACTATTAGACATTTTTAAGGAAACCTCTTTGTTAtct                                                         <  1:817536/140‑1 (MQ=255)
                                                                              tgtTTTTTTATGGGAGTAGTCTGTACTACCTGTGACTCATTGTTTTTGCCTTGACAAGGTTTTTTCCTGGTAGAACGCTGCTCATCTCAAAATTAGTAACACTATTAGACATTTTTAAGGAAACCTCTTTGTTATCTCaa                                                      <  1:636219/140‑1 (MQ=255)
                                                                                   ttttATGGGAGTAGTCTGTACTACCTGTGACTCATTGTTTTTGCCTTGACAAGGTTTTTTCCTGGTAGAACGCTGCTCATCTCAAAATTAGTAACACTATTAGACATTTTTAAGGAAACCTCTTTGTTAtctc                                                        >  1:736845/1‑133 (MQ=255)
                                                                                   ttttATGGGAGTAGTCTGTACTACCTGTGACTCATTGTTTTTGCCTTGACAAGGTTTTTTCCTGGTAGAACGCTGCTCATCTCAAAATTAGTAACACTATTAGACATTTTTAAGGAAACCTCTTTGTTAtctc                                                        <  2:736845/133‑1 (MQ=255)
                                                                                              tagtCTGTACTACCTGTGACTCATTGTTTTTGCCTTGACAAGGTTTTTTCCTGGTAGAACGCTGCTCATCTCAAAATTAGTAACACTATTAGACATTTTTAAGGAAACCTCTTTGTTATCTCaa                                                      >  1:117367/1‑124 (MQ=255)
                                                                                              tagtCTGTACTACCTGTGACTCATTGTTTTTGCCTTGACAAGGTTTTTTCCTGGTAGAACGCTGCTCATCTCAAAATTAGTAACACTATTAGACATTTTTAAGGAAACCTCTTTGTTATCTCaa                                                      <  2:117367/124‑1 (MQ=255)
                                                                                                 tCTGTACTACCTGTGACCCATTTTTTTTGCCTTGACAAGATTTTTTCCTGGTAGAACGCTGCTCATCTCAAAATTAGTAACACTATTAGACATTTTTAAGGAAACCTCTTTGTTATCTCAATATAAAACTACTATCTGaa                                   <  2:258503/140‑1 (MQ=255)
                                                                                                  cTGTACTACCTGTGACTCATTGTTTTTGCCTTGACAAGGTTTTTTCCTGGTAGAACGCTGCTCATCTCAAAATTAGTAACACTATTAGACATTTTTAAGGAAACCTCTTTGTTATCTCAATATAAAACTACTATCTGaaa                                  <  1:664539/140‑1 (MQ=255)
                                                                                                      actacCTGTGACTCATTGTTTTTGCCTTGACAAGGTTTTTTCCTGGTAGAACGCTGCTCATCTCAAAATCAGTAACACTATTAGACATTTTTAAGGAAACCTCTTTGTTATCTCAATATAAAACTACTATCTGAAATTtc                              >  1:120246/1‑140 (MQ=255)
                                                                                                          ccTGTGACTCATTGTTTTTGCCTTGACAAGGTTTTTTCCTGGTAGAACGCTGCTCATCTCAAAATTAGTAACACTATTAGACATTTTTAAGGAAACCTTTTTGTTATCTCAATATAAAACTACTATCTGAAATTTCTCAg                          >  1:713687/1‑140 (MQ=255)
                                                                                                          ccTGTGACTCATTGTTTTTGCCTTGACAAGGTTTTTTCCTGGTAGAACGCTGCTCATCTCAAAATTAGTAACACTATTAGACATTTTTAAGGAAACCTCTTTGTTATCTCAATATAAAACTACTATCTGAAATTTCTCAg                          >  1:635416/1‑140 (MQ=255)
                                                                                                          ccTGTGACTCATTGTTTTTGCCTTGACAAGGTTTTTTCCTGGTAGAACGCTGCTCATCTCAAAATTAGTAACACTATTAGACATTTTTAAGAAAACCTCTTTGTTAACGCAATATAAAACTACTATCTGAAATTTCTCAg                          >  2:380643/1‑140 (MQ=255)
                                                                                                           cTGTGACTCATTGTTTTTGCCTTGACAAGGTTTTTTCCTGGTAGAACGCTGCTCATCTCAAAATTAGTAACACTATTAGACATTTTTAAGGAAACCTCTTTGTTATCTCAGTATGAAGCTACTAGCTGAAATTTCTCAgt                         >  1:258503/1‑140 (MQ=255)
                                                                                                             gtgACTCATTGTTTTTGCCTTGACAAGGTTTTTTCCTGGTAGAACGCTGCTCATCTCAAAATTAGTAACACTATTAGACATTTTTAAGGAAACCTCTTTGTTATCTCAATATAAAACTACTATCTGAAATTTCTCAgtgt                       >  2:512596/1‑140 (MQ=255)
                                                                                                               gACTCATTGTTTTTGCCTTGACAAGGTTTTTTCCTGGTAGAACGCTGCTCATCTCAAAATTAGTAACACTATTAGACATTTTTAAGGAAACCTCTTTGTTATCTCAATATAAAACTACTATCTGAAATTTCTCAgtgt                       >  2:611850/1‑138 (MQ=255)
                                                                                                               gACTCATTGTTTTTGCCTTGACAAGGTTTTTTCCTGGTAGAACGCTGCTCATCTCAAAATTAGTAACACTATTAGACATTTTTAAGGAAACCTCTTTGTTATCTCAATATAAAACTACTATCTGAAATTTCTCAgtgt                       <  1:611850/138‑1 (MQ=255)
                                                                                                               gACTCATTGTTTTTGCCTTGACAAGGTTTTTTCCTGGTAGAACGCTGCTCATCTCAAAATTAGTAACACTATTAGACATTTTTAAGGAAACCTCTTTGTTATCTCAATATAAAACTACTATCTGAAATTTCTCAGTGTTg                     >  1:786560/1‑140 (MQ=255)
                                                                                                                    aTTGTTTTTGCCTTGACAAGGTTTTTTCCTGGTAGAACGCTGCTCATCTCAAAATTAGTAACACTATTAGACATTTTTAAGGAAACCTCTTTGTTATCTCAATATAAAACTACTATCTGAAATTTCTCAGTGTTGCCCCa                >  1:736211/1‑140 (MQ=255)
                                                                                                                                ttGACAAGGTTTTTTCCTGGTAGAACGCTGCTCATCTCAAAATTAGTAACACTATTAGACATTTTTAAGGAAACCTCTTTGTTATCTCAATATAAAACTACTATCTGAAATTTCTCAGTGTTGCCCCATTTTTCTCgggg    <  1:135290/140‑1 (MQ=255)
                                                                                                                                  gACAAGGTTTTTTCCTGGTAGAACGCTGCTCATCTCAAAATTAGTAACACTATTAGACATTTTTAAGGAAACCTCTTTGTTATCTCAATATAAAACTACTATTTGAAATTTCTCAGTATTGCCCCATTTTTCTCGGGGCa  >  1:301040/1‑140 (MQ=255)
                                                                                                                                  gACAAGGTTTTTTCCTGGTAGAACGCTGCTCATCTCAAAATTAGTAACACTATTAGACATTTTTAAGGAAACCTCTTTGTTATCTCAATATAAAACTACTATCTGAAATTTCTCAGTGTTGCCCCATTTTTCTCGGGGCa  >  2:31102/1‑140 (MQ=255)
                                                                                                                                         |                                                                                                                                     
GAAGAGTGAGATTTATTTGTTAGTCTGGCATAATGATCAATTACTCCTAATGAAGATTGATTTTTGAACGTGACTGAGTGTTTTTTTATGGGAGTAGTCTGTACTACCTGTGACTCATTGTTTTTGCCTTGACAAGGTTTTTTTCCTGGTAGAACGCTGCTCATCTCAAAATTAGTAACACTATTAGACATTTTTAAGGAAACCTCTTTGTTATCTCAATATAAAACTACTATCTGAAATTTCTCAGTGTTGCCCCATTTTTCTCGGGGCA  >  NC_000913/294609‑294879

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