Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,777,746 G→T R98S (AGG→AGT ypjF → CP4‑57 prophage; toxin of the YpjF‑YfjZ toxin‑antitoxin system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,777,7460GT99.1% 381.2 / ‑4.1 109R98S (AGG→AGTypjFCP4‑57 prophage; toxin of the YpjF‑YfjZ toxin‑antitoxin system
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base T (66/42);  minor base A (0/1);  total (66/43)
Fisher's exact test for biased strand distribution p-value = 3.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.31e-01

Only 100 of 109 total aligned reads displayed.

ATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTTCCGTAGTAGATATTCTGCGAGCAAGGCGCTCTACCGGCTTGCTAAAAACTAATGTGAAATAAACGCTTAAATACAGAGCAGACTGAAGGAAAGCAAAATGCTAATCTCATAGACGAAGAGACTCCCGCTGTAACCTCCCCTCCCCCGAAAAAACCTGACATT  >  NC_000913/2777613‑2777882
                                                                                                                                     |                                                                                                                                        
aTGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTTCCGTAGTAGATATTCTGCGAGCAAGTCGctct                                                                                                                                    >  2:1680118/1‑140 (MQ=255)
 tGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTTCCGTAGTAGATATTCTGCGAGCAAGTCGCTCTa                                                                                                                                   >  2:71348/1‑140 (MQ=255)
 tGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTTCCGTAGTAGATATTCTGCGAGCAAGTCGCTCTa                                                                                                                                   <  1:267584/140‑1 (MQ=255)
  gCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTTCCGTAGTAGATATTCTGCGAGCAAGTCGCTCTAc                                                                                                                                  >  1:1513556/1‑140 (MQ=255)
    tGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTTCCGTAGTAGATATTCTGCGAGCAAGTCGCTCTACCg                                                                                                                                >  1:1087775/1‑140 (MQ=255)
    tGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTTCCGTAGTAGATATTCTGCGAGCAAGTCGCTCTACCg                                                                                                                                >  1:1396169/1‑140 (MQ=255)
    tGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTTCCGTAGTAGATATTCTGCGAGCAAGTCGCTCTACCg                                                                                                                                >  1:528835/1‑140 (MQ=255)
      gTATCACTCTGGCCGATGCAGTCAATTTTCTCGTGGAAAAGTACGAACTGGTACGCAGCGATCACAGAGGATTTTCCTCGCAACAACAGTCTCCATATATTTCCGTAGTAGATATTCTGCGAGCAAGTCGCTCTACCGGc                                                                                                                              <  2:1300497/140‑1 (MQ=255)
       tATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTTCCGTAGTAGATATTCTGCGAGCAAGTCGCTCTACCGGCt                                                                                                                             <  1:900795/140‑1 (MQ=255)
        aTCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTTCCGTAGTAGATATTCTGCGAGCAAGTCGCTCTACCGGCtt                                                                                                                            >  1:602163/1‑140 (MQ=255)
           aCTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTTCCGTAGTAGATATTCTGCGAGCAAGTCGCTCTACCGGCTTGCt                                                                                                                         <  2:1087775/140‑1 (MQ=255)
           aCTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTTCCGTAGTAGATATTCTGCGAGCAAGTCGCTCTACCGGCTTGCt                                                                                                                         >  2:890007/1‑140 (MQ=255)
           aCTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTTCCGTAGTAGATATTCTGCGAGCAAGTCGCTCTACCGGCTTGCt                                                                                                                         >  2:526767/1‑140 (MQ=255)
           aCTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTTCCGTAGTAGATATTCTGCGAGCAAGTCGCTCTACCGGCTTGCt                                                                                                                         >  2:495385/1‑140 (MQ=255)
           aCTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTTCCGTAGTAGATATTCTGCGAGCAAGTCGCTCTACCGGCTTGCt                                                                                                                         <  1:988221/140‑1 (MQ=255)
              ctGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTTCCGTAGTAGATATTCTGCGAGCAAGTCGCTCTACCGGCTTGCTaaa                                                                                                                      >  1:751823/1‑140 (MQ=255)
               tGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTTCCGTAGTAGATATTCTGCGAGCAAGTCGCTCTACCGGCTTGCTaaaa                                                                                                                     >  1:221606/1‑140 (MQ=255)
                ggCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTTCCGTAGTAGATATTCTGCGAGCAAGTCGCTCTACCGGCTTGCTaaaaa                                                                                                                    <  2:1308000/140‑1 (MQ=255)
                ggCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTTCCGTAGTAGATATTCTGCGAGCAAGTCGCTCTACCGGCTTGCTaaaaa                                                                                                                    >  1:1300497/1‑140 (MQ=255)
                 gCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTTCCGTAGTAGATATTCTGCGAGCAAGTCGCTCTACCGGCTTGCTAAAAAc                                                                                                                   >  2:1571518/1‑140 (MQ=255)
                 gCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTTCCGTAGTAGATATTCTGCGAGCAAGTCGCTCTACCGGCTTGCTAAAAAc                                                                                                                   <  2:420857/140‑1 (MQ=255)
                 gCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTTCCGTAGTAGATATTCTGCGAGCAAGTCGCTCTACCGGCTTGCTAAAAAc                                                                                                                   >  2:416683/1‑140 (MQ=255)
                      tGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTTCCGTAGTAGATATTCTGCGAGCAAGTCGCTCTACCGGCTTGCTAAAAACTAAtg                                                                                                              >  2:180118/1‑140 (MQ=255)
                      tGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTTCCGTAGTAGATATTCTGCGAGCAAGTCGCTCTACCGGCTTGCTAAAAACTAAtg                                                                                                              <  1:1680118/140‑1 (MQ=255)
                           tCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTTCCGTAGTAGATATTCTGCGAGCAAGTCGCTCTACCGGCTTGCTAAAAACTAATGTGaaa                                                                                                         >  1:620233/1‑140 (MQ=255)
                           tCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTTCCGTAGTAGATATTCTGCGAGCAAGTCGCTCTACCGGCTTGCTAAAAACTAATGTGaaa                                                                                                         >  1:1474067/1‑140 (MQ=255)
                               ttttCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTTCCGTAGTAGATATTCTGCGAGCAAGTCGCTCTACCGGCTTGCTAAAAACTAATGTGAAATaaa                                                                                                     <  2:401340/140‑1 (MQ=255)
                               ttttCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTTCCGTAGTAGATATTCTGCGAGCAAGTCGCTCTACCGGCTTGCTAAAAACTAATGTGAAATaaa                                                                                                     <  2:1511369/140‑1 (MQ=255)
                               ttttCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTTCCGTAGTAGATATTCTGCGAGCAAGTCGCTCTACCGGCTTGCTAAAAACTAATGTGAAATaaa                                                                                                     <  1:526767/140‑1 (MQ=255)
                               ttttCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTTCCGTAGTAGATATTCTGCGAGCAAGTCGCTCTACCGGCTTGCTAAAAACTAATGTGAAATaaa                                                                                                     <  2:685898/140‑1 (MQ=255)
                               ttttCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTTCCGTAGTAGATATTCTGCGAGCAAGTCGCTCTACCGGCTTGCTAAAAACTAATGTGAAATaaa                                                                                                     <  2:620233/140‑1 (MQ=255)
                               ttttCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTTCCGTAGTAGATATTCTGCGAGCAAGACGCTCTACCGGCTTGCTAAAAACTAATGTGAAATaaa                                                                                                     <  2:410932/140‑1 (MQ=255)
                                    tggtggAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTTCCGTAGTAGATATTCTGCGAGCAAGTCGCTCTACCGGCTTGCTAAAAACTAAt                                                                                                               <  1:1514468/125‑1 (MQ=255)
                                    tggtggAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTTCCGTAGTAGATATTCTGCGAGCAAGTCGCTCTACCGGCTTGCTAAAAACTAAt                                                                                                               >  2:1514468/1‑125 (MQ=255)
                                    tggtggAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTTCCGTAGTAGATATTCTGCGAGCAAGTCGCTCTACCGGCTTGCTAAAAACTAATGTGAAATAAACGCtt                                                                                                <  2:602163/140‑1 (MQ=255)
                                    tggtggAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTTCCGTAGTAGATATTCTGCGAGCAAGTCGCTCTACCGGCTTGCTAAAAACTAATGTGAAATAAACGCtt                                                                                                <  1:248400/140‑1 (MQ=255)
                                           aaaGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTTCCGTAGTAGATATTCTGCGAGCAAGTCGCTCTACCGGCTTGCTAAAAACTAATGTGAAATAAACGCTTAAATACa                                                                                         >  1:866826/1‑140 (MQ=255)
                                           aaaGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTTCCGTAGTAGATATTCTGCGAGCAAGTCGCTCTACCGGCTTGCTAAAAACTAATGTGAAATAAACGCTTAAATACa                                                                                         <  1:973353/140‑1 (MQ=255)
                                              gTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTTCCGTAGTAGATATTCTGCGAGCAAGTCGCTCTACCGGCTTGCTAAAAACTAATGTGAAATAAACGCTTAAATACagag                                                                                      <  1:180118/140‑1 (MQ=255)
                                              gTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTTCCGTAGTAGATATTCTGCGAGCAAGTCGCTCTACCGGCTTGCTAAAAACTAATGTGAAATAAACGCTTAAATACagag                                                                                      <  2:1474067/140‑1 (MQ=255)
                                              gTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTTCCGTAGTAGATATTCTGCGAGCAAGTCGCTCTACCGGCTTGCTAAAAACTAATGTGAAATAAACGCTTAAATACagag                                                                                      <  1:495385/140‑1 (MQ=255)
                                              gTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTTCCGTAGTAGATATTCTGCGAGCAAGTCGCTCTACCGGCTTGCTAAAAACTAATGTGAAATAAACGCTTAAATACagag                                                                                      <  1:751534/140‑1 (MQ=255)
                                                aCGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTTCCGTAGTAGATATTCTGCGAGCAAGTCGCTCTACCGGCTTGCTAAAAACTAATGTGAAATAAACGCTTAAATACAGAGCa                                                                                    <  2:602688/140‑1 (MQ=255)
                                                 cGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTTCCGTAGTAGATATTCTGCGAGCAAGTCGCTCTACCGGCTTGCTAAAAACTAATGTGAAATAAACGCTTAAATACAGAGCAg                                                                                   >  2:540717/1‑140 (MQ=255)
                                                         tACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTTCCGTAGTAGATATTCTGCGAGCAAGTCGCTCTACCGGCTTGCTAAAAACTAATGTGAAATAAACGCTTAAATACAGAGCAGACTGAAgg                                                                           >  2:253926/1‑140 (MQ=255)
                                                         tACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTTCCGTAGTAGATATTCTGCGAGCAAGTCGCTCTACCGGCTTGCTAAAAACTAATGTGAAATAAACGCTTAAATACAGAGCAGACTGAAgg                                                                           >  1:1039847/1‑140 (MQ=255)
                                                         tACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTTCCGTAGTAGATATTCTGCGAGCAAGTCGCTCTACCGGCTTGCTAAAAACTAATGTGAAATAAACGCTTAAATACAGAGCAGACTGAAgg                                                                           >  2:718505/1‑140 (MQ=255)
                                                         tACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTTCCGTAGTAGATATTCTGCGAGCAAGTCGCTCTACCGGCTTGCTAAAAACTAATGTGAAATAAACGCTTAAATACAGAGCAGACTGAAgg                                                                           >  1:528893/1‑140 (MQ=255)
                                                             tATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTTCCGTAGTAGATATTCTGCGAGCAAGTCGCTCTACCGGCTTGCTAAAAACTAATGTGAAATAAACGCTTAAATACAGAGCAGACTGAAGGAAAg                                                                       >  1:1595294/1‑140 (MQ=255)
                                                                  aTCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTTCCGTAGTAGATATTCTGCGAGCAAGTCGCTCTACCGGCTTGCTAAAAACTAATGTGAAATAAACGCTTAAATACAGAGCAGACTGAAGGa                                                                          <  2:501957/132‑1 (MQ=255)
                                                                  aTCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTTCCGTAGTAGATATTCTGCGAGCAAGTCGCTCTACCGGCTTGCTAAAAACTAATGTGAAATAAACGCTTAAATACAGAGCAGACTGAAGGa                                                                          >  1:501957/1‑132 (MQ=255)
                                                                  aTCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTTCCGTAGTAGATATTCTGCGAGCAAGTCGCTCTACCGGCTTGCTAAAAACTAATGTGAAATAAACGCTTAAATACAGAGCAGACTGAAGGAAAGCaaaa                                                                  >  1:150187/1‑140 (MQ=255)
                                                                  aTCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTTCCGTAGTAGATATTCTGCGAGCAAGTCGCTCTACCGGCTTGCTAAAAACTAATGTGAAATAAACGCTTAAATACAGAGCAGACTGAAGGAAAGCaaaa                                                                  >  2:67439/1‑140 (MQ=255)
                                                                  aTCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTTCCGTAGTAGATATTCTGCGAGCAAGTCGCTCTACCGGCTTGCTAAAAACTAATGTGAAATAAACGCTTAAATACAGAGCAGACTGAAGGAAAGCaaaa                                                                  >  1:1610385/1‑140 (MQ=255)
                                                                  aTCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTTCCGTAGTAGATATTCTGCGAGCAAGTCGCTCTACCGGCTTGCTAAAAACTAATGTGAAATAAACGCTTAAATACAGAGCAGACTGAAGGAAAGCaaaa                                                                  >  2:1670683/1‑140 (MQ=255)
                                                                   tCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTTCCGTAGTAGATATTCTGCGAGCAAGTCGCTCTACCGGCTTGCTAAAAACTAATGTGAAATAAACGCTTAAATACAGAGCAGACTGAAGGAAAGCAAAAt                                                                 >  2:1245289/1‑140 (MQ=255)
                                                                       agagGATTTTCGTGGCAACAACAGTCTCCATATATTTCCGTAGTAGATATTCTGCGAGCAAGTCGCTCTACCGGCTTGCTAAAAACTAATGTGAAATAAACGCTTAAATACAGAGCAGACTGAAGGAAAGCAAAATGCTa                                                             <  1:1669524/140‑1 (MQ=255)
                                                                             ttttCGTGGCAACAACAGTCTCCATATATTTCCGTAGTAGATATTCTGCGAGCAAGTCGCTCTACCGGCTTGCTAAAAACTAATGTGAAATAAACGCTTAAATACAGAGCAGACTGAAGGAAAGCAAAATGCTAATCTCa                                                       <  2:1039847/140‑1 (MQ=255)
                                                                             ttttCGTGGCAACAACAGTCTCCATATATTTCCGTAGTAGATATTCTGCGAGCAAGTCGCTCTACCGGCTTGCTAAAAACTAATGTGAAATAAACGCTTAAATACAGAGCAGACTGAAGGAAAGCAAAATGCTAATCTCa                                                       >  2:838385/1‑140 (MQ=255)
                                                                             ttttCGTGGCAACAACAGTCTCCATATATTTCCGTAGTAGATATTCTGCGAGCAAGTCGCTCTACCGGCTTGCTAAAAACTAATGTGAAATAAACGCTTAAATACAGAGCAGACTGAAGGAAAGCAAAATGCTAATCTCa                                                       >  1:223892/1‑140 (MQ=255)
                                                                             ttttCGTGGCAACAACAGTCTCCATATATTTCCGTAGTAGATATTCTGCGAGCAAGTCGCTCTACCGGCTTGCTAAAAACTAATGTGAAATAAACGCTTAAATACAGAGCAGACTGAAGGAAAGCAAAATGCTAATCTCa                                                       >  1:466264/1‑140 (MQ=255)
                                                                             ttttCGTGGCAACAACAGTCTCCATATATTTCCGTAGTAGATATTCTGCGAGCAAGTCGCTCTACCGGCTTGCTAAAAACTAATGTGAAATAAACGCTTAAATACAGAGCAGACTGAAGGAAAGCAAAATGCTAATCTCa                                                       >  2:1163761/1‑140 (MQ=255)
                                                                               ttCGTGGCAACAACAGTCTCCATATATTTCCGTAGTAGATATTCTGCGAGCAAGTCGCTCTACCGGCTTGCTAAAAACTAATGTGAAATAAACGCTTAAATACAGAGCAGACTGAAGGAAAGCAAAATGCTAATCTCATa                                                     <  1:113680/140‑1 (MQ=255)
                                                                               ttCGTGGCAACAACAGTCTCCATATATTTCCGTAGTAGATATTCTGCGAGCAAGTCGCTCTACCGGCTTGCTAAAAACTAATGTGAAATAAACGCTTAAATACAGAGCAGACTGAAGGAAAGCAAAATGCTAATCTCATa                                                     <  2:751823/140‑1 (MQ=255)
                                                                                   tGGCAACAACAGTCTCCATATATTTCCGTAGTAGATATTCTGCGAGCAAGTCGCTCTACCGGCTTGCTAAAAACTAATGTGAAATAAACGCTTAAATACAGAGCAGACTGAAGGAAAGCAAAATGCTAATCTCATAGACg                                                 >  1:1188153/1‑140 (MQ=255)
                                                                                          aacaGTCTCCATATATTTCCGTAGTAGATATTCTGCGAGCAAGTCGCTCTACCGGCTTGCTAAAAACTAATGTGAAATAAACGCTTAAATACAGAGCAGACTGAAGGAAAGCAAAATGCTAATCTCATAGACGAAGAGAc                                          <  1:1571518/140‑1 (MQ=255)
                                                                                            caGTCTCCATATATTTCCGTAGTAGATATTCTGCGAGCAAGTCGCTCTACCGGCTTGCTAAAAACTAATGTGAAATAAACGCTTAAATACAGAGCAGACTGAAGGAAAGCAAAATGCTAATCTCATAGACGAAGAGACTc                                        <  1:540717/140‑1 (MQ=255)
                                                                                             aGTCTCCATATATTTCCGTAGTAGATATTCTGCGAGCAAGTCGCTCTACCGGCTTGCTAAAAACTAATGTGAAATAAACGCTTAAATACAGAGCAGACTGAAGGAAAGCAAAATGCTAATCTCATAGACGAAGAGACTcc                                       <  1:1163761/140‑1 (MQ=255)
                                                                                                ctcCATATATTTCCGTAGTAGATATTCTGCGAGCAAGTCGCTCTACCGGCTTGCTAAAAACTAATGTGAAATAAACGCTTAAATACAGAGCAGACTGAAGGAAAGCAAAATGCTAATCTCATAGACGAAGAGACTCCCGc                                    <  2:1595294/140‑1 (MQ=255)
                                                                                                    atatatTTCCGTAGTAGATATTCTGCGAGCAAGTCGCTCTACCGGCTTGCTAAAAACTAATGTGAAATAAACGCTTAAATACAGAGCAGACTGAAGGAAAGCAAAATGCTAATCTCATAGACGAAGAGACTCCCGCTGTa                                <  2:1610385/140‑1 (MQ=255)
                                                                                                    atatatTTCCGTAGTAGATATTCTGCGAGCAAGTCGCTCTACCGGCTTGCTAAAAACTAATGTGAAATAAACGCTTAAATACAGAGCAGACTGAAGGAAAGCAAAATGCTAATCTCATAGACGAAGAGACTCCCGCTGTa                                >  2:1646369/1‑140 (MQ=255)
                                                                                                    atatatTTCCGTAGTAGATATTCTGCGAGCAAGTCGCTCTACCGGCTTGCTAAAAACTAATGTGAAATAAACGCTTAAATACAGAGCAGACTGAAGGAAAGCAAAATGCTAATCTCATAGACGAAGAGACTCCCGCTGTa                                <  2:866826/140‑1 (MQ=255)
                                                                                                    atatatTTCCGTAGTAGATATTCTGCGAGCAAGTCGCTCTACCGGCTTGCTAAAAACTAATGTGAAATAAACGCTTAAATACAGAGCAGACTGAAGGAAAGCAAAATGCTAATCTCATAGACGAAGAGACTCCCGCTGTa                                >  1:1528441/1‑140 (MQ=255)
                                                                                                    atatatTTCCGTAGTAGATATTCTGCGAGCAAGTCGCTCTACCGGCTTGCTAAAAACTAATGTGAAATAAACGCTTAAATACAGAGCAGACTGAAGGAAAGCAAAATGCTAATCTCATAGACGAAGAGACCCCCGCt                                   >  1:1380036/1‑137 (MQ=255)
                                                                                                      atatTTCCGTAGTAGATATTCTGCGAGCAAGTCGCTCTACCGGCTTGCTAAAAACTAATGTGAAATAAACGCTTAAATACAGAGCAGACTGAAGGAAAGCAAAATGCTAATCTCATAGACGAAGAGACTCCCGCTGTAAc                              <  2:223892/140‑1 (MQ=255)
                                                                                                      atatTTCCGTAGTAGATATTCTGCGAGCAAGTCGCTCTACCGGCTTGCTAAAAACTAATGTGAAATAAACGCTTAAATACAGAGCAGACTGAAGGAAAGCAAAATGCTAATCTCATAGACGAAGAGACTCCCGCTGTAAc                              <  2:466264/140‑1 (MQ=255)
                                                                                                         tttCCGTAGTAGATATTCTGCGAGCAAGTCGCTCTACCGGCTTGCTAAAAACTAATGTGAAATAAACGCTTAAATACAGAGCAGACTGAAGGAAAGCAAAATGCTAATCTCATAGACGAAGAGACTCCCGCTGTAAcctc                           >  1:429010/1‑140 (MQ=255)
                                                                                                           tCCGTAGTAGATATTCTGCGAGCAAGTCGCTCTACCGGCTTGCTAAAAACTAATGTGAAATAAACGCTTAAATACAGAGCAGACTGAAGGAAAGCAAAATGCTAATCTCATAGACGAAGAGACTCCCGCTGTAAcctccc                         >  2:695523/1‑140 (MQ=255)
                                                                                                             cGTAGTAGATATTCTGCGAGCAAGTCGCTCTACCGGCTTGCTAAAAACTAATGTGAAATAAACGCTTAAATACAGAGCAGACTGAAGGAAAGCAAAATGCTAATCTCATAGACGAAGAGACTCCCGCTGTAAcctcccct                       >  2:145950/1‑140 (MQ=255)
                                                                                                             cGTAGTAGATATTCTGCGAGCAAGTCGCTCTACCGGCTTGCTAAAAACTAATGTGAAATAAACGCTTAAATACAGAGCAGACTGAAGGAAAGCAAAATGCTAATCTCATAGACGAAGAGACTCCCGCTGTAAcctccccc                       >  2:849724/1‑139 (MQ=255)
                                                                                                               tagtagATATTCTGCGAGCAAGTCGCTCTACCGGCTTGCTAAAAACTAATGTGAAATAAACGCTTAAATACAGAGCAGACTGAAGGAAAGCAAAATGCTAATCTCATAGACGAAGAGACTcc                                       <  2:610428/122‑1 (MQ=255)
                                                                                                               tagtagATATTCTGCGAGCAAGTCGCTCTACCGGCTTGCTAAAAACTAATGTGAAATAAACGCTTAAATACAGAGCAGACTGAAGGAAAGCAAAATGCTAATCTCATAGACGAAGAGACTcc                                       >  1:610428/1‑122 (MQ=255)
                                                                                                               tagtagATATTCTGCGAGCAAGTCGCTCTACCGGCTTGCTAAAAACTAATGTGAAATAAACGCTTAAATACAGAGCAGACTGAAGGAAAGCAAAATGCTAATCTCATAGACGAAGAGACTCCCGCTGTAAcctcccctcc                     >  1:1002262/1‑140 (MQ=255)
                                                                                                                 gtagATATTCTGCGAGCAAGTCGCTCTACCGGCTTGCTAAAAACTAATGTGAAATAAACGCTTAAATACAGAGCAGACTGAAGGAAAGCAAAATGCTAATCTCATAGACGAAGAGACTCCCGCTGTAACCtcccctgccc                   >  1:1530708/1‑140 (MQ=255)
                                                                                                                    gATATTCTGCGAGCAAGTCGCTCTACCGGCTTGCTAAAAACTAATGTGAAATAAACGCTTAAATACAGAGCAGACTGAAGGAAAGCAAAATGCTAATCTCATAGACGAAGAGACTCCCGCTGTAACctcccctccc                    <  2:1431584/136‑1 (MQ=255)
                                                                                                                    gATATTCTGCGAGCAAGTCGCTCTACCGGCTTGCTAAAAACTAATGTGAAATAAACGCTTAAATACAGAGCAGACTGAAGGAAAGCAAAATGCTAATCTCATAGACGAAGAGACTCCCGCTGTAACCTCCCCTCCCCCGa                >  2:524138/1‑140 (MQ=255)
                                                                                                                       atTCTGCGAGCAAGTCGCTCTACCGGCTTGCTAAAAACTAATGTGAAATAAACGCTTAAATACAGAGCAGACTGAAGGAAAGCAAAATGCTAATCTCATAGACGAAGAGACTCCCGCTGTAACCTCCCCTCCCCCGaaaa             <  1:1245289/140‑1 (MQ=255)
                                                                                                                         tCTGCGAGCAAGTCGCTCTACCGGCTTGCTAAAAACTAATGTGAAATAAACGCTTAAATACAGAGCAGACTGAAGGAAAGCAAAATGCTAATCTCATAGACGAAGAGACTCCCGCTGTAACCTCCCCTCCCCCGaaaaaa           >  1:1058668/1‑140 (MQ=255)
                                                                                                                            gCGAGCAAGTCGCTCTACCGGCTTGCTAAAAACTAATGTGAGATAAACGCTTAAATACAGAGCAGACTGAAGGAAAGCAAAATGCTAATCTCATAGACGAAGAGACTCCCGCTGTAACCTCCCCTCCCCCGAAAAAACCt        >  2:1463138/1‑140 (MQ=255)
                                                                                                                            gCGAGCAAGTCGCTCTACCGGCTTGCTAAAAACTAATGTGAAATAAACGCTTAAATACAGAGCAGACTGAAGGAAAGCAAAATGCTAATCTCATAGACGAAGAGACTCCCGCTGTAACCTCCCCTCCCCCGAAAAAACCt        >  2:797087/1‑140 (MQ=255)
                                                                                                                            gCGAGCAAGTCGCTCTACCGGCTTGCTAAAAACTAATGTGAAATAAACGCTTAAATACAGAGCAGACTGAAGGAAAGCAAAATGCTAATCTCATAGACGAAGAGACTCCCGCTGTAACCTCCCCTCCCCCGAAAAAACCt        >  1:99994/1‑140 (MQ=255)
                                                                                                                            gCGAGCAAGTCGCTCTACCGGCTTGCTAAAAACTAATGTGAAATAAACGCTTAAATACAGAGCAGACTGAAGGAAAGCAAAATGCTAATCTCATAGACGAAGAGACTCCCGCTGTAACCTCCCCTCCCCCGAAAAAACCt        >  2:1064547/1‑140 (MQ=255)
                                                                                                                            gCGAGCAAGTCGCTCTACCGGCTTGCTAAAAACTAATGTGAAATAAACGCTTAAATACAGAGCAGACTGAAGGAAAGCAAAATGCTAATCTCATAGACGAAGAGACTCCCGCTGTAACCTCCCCTCCCCCGAAAAAACCt        >  1:455167/1‑140 (MQ=255)
                                                                                                                            gCGAGCAAGTCGCTCTACCGGCTTGCTAAAAACTAATGTGAAATAAACGCTTAAATACAGAGCAGACTGAAGGAAAGCAAAATGCTAATCTCATAGACGAAGAGACTCCCGCTGTAACCTCCCCTCCCCCGAAAAAACCt        >  2:357708/1‑140 (MQ=255)
                                                                                                                            gCGAGCAAGTCGCTCTACCGGCTTGCTAAAAACTAATGTGAAATAAACGCTTAAATACAGAGCAGACTGAAGGAAAGCAAAATGCTAATCTCATAGACGAAGAGACTCCCGCTGTAACCTCCCCTCCCCCGAAAAAACCt        >  2:234041/1‑140 (MQ=255)
                                                                                                                              gAGCAAGTCGCTCTACCGGCTTGCTAAAAACTAATGTGAAATAAACGCTTAAATACAGAGCAGACTGAAGGAAAGCAAAATGCTAATCTCATAGACGAAGAGACTCCCGCTGTAACCTCCCCTCCCCCGAAAAAACCTGa      >  2:612329/1‑140 (MQ=255)
                                                                                                                               aGCAAGTCGCTCTACCGGCTTGCTAAAAACTAATGTGAAATAAACGCTTAAATACAGAGCAGACTGAAGGAAAGCAAAATGCTAATCTCATAGACGAAGAGACTCCCGCTGTAACCTCCCCTCCCCCGAAAAAACCTGAc     >  2:1244897/1‑140 (MQ=255)
                                                                                                                               aGCAAGTCGCTCTACCGGCTTGCTAAAAACTAATGTGAAATAAACGCTTAAATACAGAGCAGACTGAAGGAAAGCAAAATGCTAATCTCATAGACGAAGAGACTCCCGCTGTAACCTCCCCTCCCCCGAAAAAACCTGAc     >  2:1393320/1‑140 (MQ=255)
                                                                                                                               aGCAAGTCGCTCTACCGGCTTGCTAAAAACTAATGTGAAATAAACGCTTAAATACAGAGCAGACTGAAGGAAAGCAAAATGCTAATCTCATAGACGAAGAGACTCCCGCTGTAACCTCCCCTCCCCCGAAAAAACCTGAc     >  1:1451652/1‑140 (MQ=255)
                                                                                                                                  aaGTCGCTCTACCGGCTTGCTAAAAACTAATGTGAAATAAACGCTTAAATACAGAGCAGACTGAAGGAAAGCAAAATGCTAATCTCATAGACGAAGAGACTCCCGCTGTAACCTCCCCTCCCCCGAAAAAACCTGACAtt  >  1:1264127/1‑140 (MQ=255)
                                                                                                                                     |                                                                                                                                        
ATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTTCCGTAGTAGATATTCTGCGAGCAAGGCGCTCTACCGGCTTGCTAAAAACTAATGTGAAATAAACGCTTAAATACAGAGCAGACTGAAGGAAAGCAAAATGCTAATCTCATAGACGAAGAGACTCCCGCTGTAACCTCCCCTCCCCCGAAAAAACCTGACATT  >  NC_000913/2777613‑2777882

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