Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,612,330 G→C 5.7% intergenic (+41/‑81) ytjA → / → yjjU uncharacterized protein/putative patatin‑like family phospholipase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,612,3300GC 5.7% 364.1 / 8.2 104intergenic (+41/‑81)ytjA/yjjUuncharacterized protein/putative patatin‑like family phospholipase
Reads supporting (aligned to +/- strand):  ref base G (49/49);  new base C (4/2);  total (53/51)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 104 total aligned reads displayed.

TTTGGTGGTCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGCCGTTCCAGCCTG  >  NC_000913/4612185‑4612477
                                                                                                                                                 |                                                                                                                                                    
    gtggtCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTg                                                                                                                                              >  2:102720/1‑149 (MQ=255)
         cTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCGCT‑GGACTGGCttt                                                                                                                                         <  3:384586/149‑1 (MQ=255)
         cTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCGCT‑GGACTGGCttt                                                                                                                                         <  3:531366/149‑1 (MQ=255)
         cTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCC‑GCTGGACTGGCttt                                                                                                                                         <  1:46578/149‑1 (MQ=255)
            gCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCC‑GCTGGACTGGCtttgga                                                                                                                                      <  2:661648/149‑4 (MQ=255)
             ccGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCg                                                                                                                                     >  1:943640/1‑149 (MQ=255)
             ccGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCg                                                                                                                                     <  3:815225/149‑1 (MQ=255)
              cGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCgg                                                                                                                                    >  1:310774/1‑149 (MQ=255)
              cGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCgg                                                                                                                                    <  1:312733/149‑1 (MQ=255)
               ggTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGt                                                                                                                                   >  1:299346/1‑149 (MQ=255)
                gTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGtt                                                                                                                                  <  3:301145/149‑1 (MQ=255)
                  aCCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGtttt                                                                                                                                >  2:450413/1‑149 (MQ=255)
                  aCCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGtttt                                                                                                                                <  1:450413/149‑1 (MQ=255)
                   ccGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTa                                                                                                                               >  2:499486/1‑149 (MQ=255)
                     gCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGc                                                                                                                             <  1:684924/149‑1 (MQ=255)
                     gCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGc                                                                                                                             <  4:444864/149‑1 (MQ=255)
                     gCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGc                                                                                                                             <  2:511471/149‑1 (MQ=255)
                        gCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGaa                                                                                                                          <  3:655906/149‑1 (MQ=255)
                           ggCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTgt                                                                                                                                      >  3:532118/1‑133 (MQ=255)
                           ggCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGc                                                                                                                                      <  4:532118/134‑1 (MQ=255)
                                 gCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAAtt                                                                                                                 >  2:93638/1‑149 (MQ=255)
                                   tAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGc                                                                                                               <  2:850541/149‑1 (MQ=255)
                                   tAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGc                                                                                                               <  4:62051/149‑1 (MQ=255)
                                        ttGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTa                                                                                                                               >  3:755294/1‑128 (MQ=255)
                                        ttGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTa                                                                                                                               <  4:755294/128‑1 (MQ=255)
                                        ttGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTAc                                                                                                          >  1:777071/1‑149 (MQ=255)
                                          gTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACtt                                                                                                        >  2:688259/1‑149 (MQ=255)
                                           tCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAGAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACttt                                                                                                       >  2:317569/1‑149 (MQ=255)
                                             tttGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGt                                                                                                     <  1:499486/149‑1 (MQ=255)
                                             tttGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGt                                                                                                     <  3:880548/149‑1 (MQ=255)
                                             tttGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATCTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGt                                                                                                     >  1:943920/1‑149 (MQ=255)
                                              ttGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTc                                                                                                    <  3:339724/149‑1 (MQ=255)
                                              ttGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTc                                                                                                    <  2:632527/149‑1 (MQ=255)
                                              ttGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTc                                                                                                    <  2:519633/149‑1 (MQ=255)
                                                    tcgGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCtt                                                                                              >  4:1045210/1‑149 (MQ=255)
                                                     cgGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCttt                                                                                             <  2:310774/149‑1 (MQ=255)
                                                      gggATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTc                                                                                            >  1:132155/1‑149 (MQ=255)
                                                      gggATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTc                                                                                            <  1:979808/149‑1 (MQ=255)
                                                        gATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTcac                                                                                          >  1:926838/1‑149 (MQ=255)
                                                        gATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTcac                                                                                          >  1:738567/1‑149 (MQ=255)
                                                         attattCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTcaca                                                                                         >  2:739747/1‑149 (MQ=255)
                                                          ttattCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCacaa                                                                                        >  2:1037589/1‑149 (MQ=255)
                                                           tattCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCacaac                                                                                       <  1:475796/149‑1 (MQ=255)
                                                             ttCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCAcaacaa                                                                                     >  1:1045361/1‑149 (MQ=255)
                                                              tCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACaaa                                                                                    >  1:27486/1‑149 (MQ=255)
                                                                tGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACaaaaa                                                                                  >  4:744974/1‑149 (MQ=255)
                                                                  ttCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACa                                                                                >  2:990428/1‑149 (MQ=255)
                                                                   tCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATCAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAg                                                                               <  1:968127/149‑1 (MQ=255)
                                                                   tCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAg                                                                               >  2:77399/1‑149 (MQ=255)
                                                                     cTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGa                                                                             <  2:777071/149‑1 (MQ=255)
                                                                      tGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGaa                                                                            >  4:301338/1‑149 (MQ=255)
                                                                             tttGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAg                                                                           <  2:255672/143‑1 (MQ=255)
                                                                             tttGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAg                                                                           >  1:255672/1‑143 (MQ=255)
                                                                             tttGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCagag                                                                     >  2:814455/1‑149 (MQ=255)
                                                                                    aTGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGc                                                              >  2:993584/1‑149 (MQ=255)
                                                                                    aTGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGATACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGc                                                              >  2:112313/1‑149 (MQ=255)
                                                                                       ggCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGc                                                           <  2:1045361/149‑1 (MQ=255)
                                                                                            aaaaCGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGTCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATa                                                      <  3:7324/149‑1 (MQ=255)
                                                                                            aaaaCGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACttt                                                                                                       >  2:463609/1‑100 (MQ=255)
                                                                                            aaaaCGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACttt                                                                                                       <  1:463609/100‑1 (MQ=255)
                                                                                            aaaaCGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATa                                                      <  1:551940/149‑1 (MQ=255)
                                                                                            aaaaCGACCCTAGATTTCAGACCATTAAGAAAATCTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATa                                                      <  2:943920/149‑1 (MQ=255)
                                                                                             aaaCGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATAc                                                     <  3:744974/149‑1 (MQ=255)
                                                                                             aaaCGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATAc                                                     <  1:124159/149‑1 (MQ=255)
                                                                                             aaaCGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATAc                                                     <  4:725046/149‑1 (MQ=255)
                                                                                               aCGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCt                                                   <  1:814455/149‑1 (MQ=255)
                                                                                                  aCCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTa                                                <  1:411992/149‑1 (MQ=255)
                                                                                                     cTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTgg                                                                 <  3:732082/129‑1 (MQ=255)
                                                                                                     cTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTgg                                                                 >  4:732082/1‑129 (MQ=255)
                                                                                                      tAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGc                                            >  2:862385/1‑149 (MQ=255)
                                                                                                        gATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCtt                                          <  1:1003966/149‑1 (MQ=255)
                                                                                                         aTTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTg                                         <  1:1037589/149‑1 (MQ=255)
                                                                                                              aGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAat                                    >  3:177811/1‑149 (MQ=255)
                                                                                                              aGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAat                                    >  1:44370/1‑149 (MQ=255)
                                                                                                              aGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGAGGGGCAGCGAATACCTGTAACGCGTGATAat                                    >  4:22376/1‑149 (MQ=255)
                                                                                                               gACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAata                                   >  4:377226/1‑149 (MQ=255)
                                                                                                               gACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAata                                   <  2:525454/149‑1 (MQ=255)
                                                                                                                aCCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAatat                                  <  3:175972/149‑1 (MQ=255)
                                                                                                                aCCATTAAGAAAATTTTAACCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAGATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAatat                                  <  1:471345/149‑1 (MQ=255)
                                                                                                                  cATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTgg                                                                 <  3:264717/116‑1 (MQ=255)
                                                                                                                  cATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTgg                                                                 >  4:264717/1‑116 (MQ=255)
                                                                                                                  cATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTg                                >  1:788553/1‑149 (MQ=255)
                                                                                                                   aTTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTgc                               >  4:699194/1‑149 (MQ=255)
                                                                                                                    ttAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTgcg                              >  2:614043/1‑149 (MQ=255)
                                                                                                                     tAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTgcgc                             >  2:549433/1‑149 (MQ=255)
                                                                                                                      aaGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAg                                                                               <  1:325607/98‑1 (MQ=255)
                                                                                                                      aaGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAg                                                                               >  2:325607/1‑98 (MQ=255)
                                                                                                                       aGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGACAGAGGTGGGGCAGCGAATACCTGTAACGATTGGTAATATTGCGCCg                           >  4:438829/1‑149 (MQ=255)
                                                                                                                        gAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGc                                                                                                               <  3:217326/64‑1 (MQ=255)
                                                                                                                        gAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGc                                                                                                               >  4:217326/1‑64 (MQ=255)
                                                                                                                                tATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGtt                         <  1:231366/142‑1 (MQ=255)
                                                                                                                                tATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGtt                         >  2:231366/1‑142 (MQ=255)
                                                                                                                                  gcaaaaGCCAGTCCGCT‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGGGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAg                >  2:441914/4‑149 (MQ=255)
                                                                                                                                    aaaaGCCAGTCCGCT‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAAC‑AAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGcc             >  1:920961/2‑149 (MQ=255)
                                                                                                                                     aaaGCCAGTCCCCT‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGcc             >  1:555803/1‑149 (MQ=255)
                                                                                                                                      aaGCCAGTCCGCT‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGCCg            >  2:358676/1‑149 (MQ=255)
                                                                                                                                      aaGCCAGTCC‑GCTGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGCCg            >  2:152941/1‑149 (MQ=255)
                                                                                                                                        gCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGCCGtt          >  1:154997/1‑149 (MQ=255)
                                                                                                                                           aGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGCCGTTCCa       <  1:990428/149‑1 (MQ=255)
                                                                                                                                               cAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGCCGTTCCAGCCt   <  2:299346/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
TTTGGTGGTCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGCCGTTCCAGCCTG  >  NC_000913/4612185‑4612477

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