Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,147,987 G→C 7.3% E123Q (GAG→CAG)  plsX → putative phosphate acyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,147,9870GC 7.3% 246.3 / 10.3 92E123Q (GAG→CAG) plsXputative phosphate acyltransferase
Reads supporting (aligned to +/- strand):  ref base G (38/47);  new base C (2/5);  total (40/52)
Fisher's exact test for biased strand distribution p-value = 6.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.73e-01

CCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTCTGGCTGAAGAGGT  >  NC_000913/1147843‑1148135
                                                                                                                                                |                                                                                                                                                    
ccTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTTGAggg                                                                                                                                                  >  2:191117/1‑149 (MQ=255)
ccTTCGCAAGCTATCCGCGCCAGGCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAggg                                                                                                                                                  <  2:503206/149‑1 (MQ=255)
  ttCGCAATCTATCCGCGGCAGTCGTGGGAGTTCACTGCGCGGGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGAAATACAGGGGCGCTGTTGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGa                                                                                                                                                <  2:818760/149‑1 (MQ=255)
  ttCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGa                                                                                                                                                >  2:810064/1‑149 (MQ=255)
  ttCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGa                                                                                                                                                >  1:1289057/1‑149 (MQ=255)
  ttCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGa                                                                                                                                                >  2:1328243/1‑149 (MQ=255)
    cGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAATCCCCTCCAGGGGCtt                                                                                                                                              <  2:715961/149‑1 (MQ=255)
     ccAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTg                                                                                                                                             <  1:756270/148‑1 (MQ=255)
      cAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTGCTCAAGCCCCTGGAGGGGATTGa                                                                                                                                            >  1:1422731/1‑149 (MQ=255)
      cAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAATCCCCTCCAGGGGCTTGa                                                                                                                                            <  1:288039/149‑1 (MQ=255)
       aaGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAATCCCCTCCAGGGGCTTGAg                                                                                                                                           <  1:609779/149‑1 (MQ=255)
       aaGCTATCCGCGCCAGGCGGGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAATCCCCTCCAGGGGCTTGAg                                                                                                                                           <  1:534331/149‑1 (MQ=255)
       aaGCTATCCCCGCCAGTCGTGTGAGTTCACTGCGCTTGGCCCTCGAGCTGGTGAAAGAAGGGCGAGCGCAAGCCTGTGGCAGTTCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAATCCCCTCCAGGGGCTTGAg                                                                                                                                           <  2:1002562/149‑1 (MQ=255)
          cTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGAGGCCCTCGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGTTTGAGCGt                                                                                                                                        >  1:1517887/1‑149 (MQ=255)
             tCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCg                                                                                                                                     >  2:1074337/1‑149 (MQ=255)
             tCCGCGCCAGTCATGGGAGTTCAATGCGCGTGGCCCTGGAGCTTGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGGGCGGATTGGGCTGGGAAAATTTTATCTCAAGCCAGTGGAGCGGATTGAGAGACCg                                                                                                                                     >  2:852894/1‑149 (MQ=255)
              ccGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGGCGAGCGCAAGCCTGTGTCAGGGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCgg                                                                                                                                    <  2:444511/149‑1 (MQ=255)
                 cgcCAGTCGTGGGAGTTCAAAGCGCGTGGCCCTGGAGCGGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGgcgc                                                                                                                                 <  2:1126974/149‑1 (MQ=255)
                 cgcCACCCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGATAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGgcgc                                                                                                                                 <  1:1131774/149‑1 (MQ=255)
                         gTGGGAGTTCAATGCGCGGGGCCCTGGAGCTGGTGAAAGACGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTTGCAAAATTATTCCTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACg                                                                                                                         <  1:1042086/149‑1 (MQ=255)
                          tGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAATGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACgg                                                                                                                        <  2:200864/149‑1 (MQ=255)
                          tGGGAGGTCAATGCGCGGGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACgg                                                                                                                        <  1:698467/149‑1 (MQ=255)
                          ggggAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACgg                                                                                                                        <  1:698470/148‑1 (MQ=255)
                            ggAGTTCAATGCGCGTGGCCCGGGAGCTGGGGAAAGAAGGTCGAGGGCAAGCCTGTGTCAGTGCCGGTAATAGCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTa                                                                                                                      <  2:1414589/149‑1 (MQ=255)
                               gTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGACGGTAATACCGGGGCGATGATGGGGCTGGCAAAATTATTACACAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTa                                                                                                                   >  2:952572/1‑149 (MQ=255)
                                     tGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGTGGGGATGGAGCTTCCGGCGCTGGTGACGTTATTACcacaa                                                                                                             >  2:347290/1‑148 (MQ=255)
                                      gcgcGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATc                                                                                                            >  1:1249918/1‑149 (MQ=255)
                                      gcgcGTGGCACTGGAGCTGGTGAAAGAAGGTCGAGCGCAACCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGTGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGTTTGAGCGTCCCGCGCTGGTGACGGTATTACcacaaa                                                                                                            >  1:232398/1‑147 (MQ=255)
                                       cgcgTGGCCCTGGAGCTGGTGAAAGAAGGGCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGCGACGGTATTACCACATCa                                                                                                           <  1:851358/149‑1 (MQ=255)
                                       cgcgTGGCCCTGGAGCTGGTGAAAGAACGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCa                                                                                                           <  1:480322/149‑1 (MQ=255)
                                       agcgTGTCCCTGGATCTGGTGAAAGATGGTCGAGCTCAAGACTGTGTCATTGCCGGTAATTCGGGGGCGCTTATTGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTACGGCGCTGGTTACTGTATTACCACATCa                                                                                                           <  2:1352545/148‑1 (MQ=255)
                                         cgTGGCCCTGGACCAGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGc                                                                                                         <  2:1459890/149‑1 (MQ=255)
                                           ggggccTGGAGCTGGTGAACGAAGGTCGAGCGCAAGCCTGTGGGAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCaa                                                                                                       <  2:164781/145‑1 (MQ=255)
                                              cccTGGAGCTGGTGTAAGACGGGCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAg                                                                                                    <  2:472275/149‑1 (MQ=255)
                                                cTGGAGCTGGTGAAAGAAGGGCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCTCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAggg                                                                                                  <  1:382325/149‑1 (MQ=255)
                                                   gAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCaa                                                                                               <  2:77610/149‑1 (MQ=255)
                                                        ggTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGTGCGCAGATGGGGCTGGCAAAATTATTACTCAAGCACCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGcaaaacggcaaaacgg                                                                                          >  2:128008/1‑149 (MQ=255)
                                                         gTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACggt                                                                                         <  1:60596/149‑1 (MQ=255)
                                                          tGACAGATGGTCGAGCGCAAGCCTGTGTCAGTGCCCGTAATTCCGGGGCGCTGATGGGGCTGGCAAATATATTACTCTAGCGCCTGGAGGGGATTGAGAGTGCTGCGATGGTGACGGTATTACCACATCAACAAAAGGGCAAAACggtg                                                                                        <  2:19582/149‑1 (MQ=255)
                                                          tGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAACAAAAGGGCAAAACggtg                                                                                        <  2:1502023/149‑1 (MQ=255)
                                                          tGAAAGAAGGTCGAGCGCAAGCCTGTGGCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACggtg                                                                                        <  1:1074337/149‑1 (MQ=255)
                                                            aaaGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGTGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACggtggt                                                                                      >  2:714865/1‑149 (MQ=255)
                                                            aaaGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTTACGGTATTACCACATCAGCAAAAGGGCAAAACggtggt                                                                                      >  1:318716/1‑149 (MQ=255)
                                                              agaagGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCTCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCCAAACGGTGGTcc                                                                                    >  2:1314546/1‑149 (MQ=255)
                                                              agaagGGCGAGCTCAAGCCGGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTcc                                                                                    <  2:543172/149‑1 (MQ=255)
                                                                aagGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAATTTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCtt                                                                                  <  2:1228010/149‑1 (MQ=255)
                                                                 agGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTg                                                                                 <  1:1369545/149‑1 (MQ=255)
                                                                  gggcGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGa                                                                                <  2:1562767/146‑1 (MQ=255)
                                                                  ggTCGAGCGCTAGCGTGGGTCAGTGCCGGTATTACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGATGGGATTGAGCGTCCGGCGCTGGTGACTGTATTACAACATCAGCAAAAGGGCAAAACGGTGTTCCTTGa                                                                                <  2:130798/149‑1 (MQ=255)
                                                                  ggTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGa                                                                                <  2:1289057/149‑1 (MQ=255)
                                                                   gTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGAc                                                                               >  1:1359706/1‑149 (MQ=255)
                                                                       aGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAg                                                                           <  1:395857/149‑1 (MQ=255)
                                                                       aGCGCAAGCCTGTGTCAGGGCCGGTAATACCGGGGGGATGATGGGGCTGGCAAAATTATAACTCAAGCCACTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAg                                                                           <  2:668510/149‑1 (MQ=255)
                                                                         cgcAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAggg                                                                         >  1:931352/1‑149 (MQ=255)
                                                                            aaGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGCCTTAGGGGcc                                                                      >  2:84367/1‑149 (MQ=255)
                                                                              gCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGTAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGGCCTTGACTGAGGGGCCac                                                                    >  2:1325471/1‑148 (MQ=255)
                                                                              gCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCGTTGACTTAGGGGCCaa                                                                    >  2:938967/1‑149 (MQ=255)
                                                                                 tgtgTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGt                                                                 <  1:1556795/149‑1 (MQ=255)
                                                                                 ggtgtCAGTGCCGGTAATACCGGGGCGCGGATGGGGCGGGCAAAATTCTTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTCGGGGCCAACGt                                                                 <  2:1212899/148‑1 (MQ=255)
                                                                                 ggtgtCAGTGCCCGTAATCCCGGGGCGCTGAGGGGGCTGGCAAAAATATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGt                                                                 <  1:1528222/148‑1 (MQ=255)
                                                                                  gtgtCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTc                                                                >  1:540953/1‑149 (MQ=255)
                                                                                  gtgtCAGTGCCGGTAATACCGGGGCGCGGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTc                                                                <  1:600131/149‑1 (MQ=255)
                                                                                    gtCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGa                                                              >  1:404798/1‑149 (MQ=255)
                                                                                     tCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGAt                                                             <  2:553747/149‑1 (MQ=255)
                                                                                     tCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGAt                                                             >  1:975898/1‑149 (MQ=255)
                                                                                       aGTGCCTGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATAGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATtg                                                           >  1:902725/1‑149 (MQ=255)
                                                                                        gTGCCGGTAATCCCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGAGTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATtgt                                                          <  1:658767/149‑1 (MQ=255)
                                                                                         tGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATtgtg                                                         >  1:349956/1‑149 (MQ=255)
                                                                                          gCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGa                                                        <  2:1239764/149‑1 (MQ=255)
                                                                                            cGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACa                                                      <  1:973490/149‑1 (MQ=255)
                                                                                                   aCCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAAt                                               >  1:792222/1‑149 (MQ=255)
                                                                                                     cGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTTCCACATCAGCAAAAGGGCAAAACGGTGGTGCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGc                                             <  1:574806/149‑1 (MQ=255)
                                                                                                       gggCGCTGAAGAGGCTGGCATAATTATTAATAAAGCCACTGGAGGGGATAGAGAGTCCGCAGCTGGTGACGGTATTAACACATCAGCAAAAGGGAAAAACGTTGGTCCATGACTTAGGGTCCAACGTCGATTGTGACAGCACAAtcatg                                           >  2:736614/1‑145 (MQ=255)
                                                                                                          cgcTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTg                                        >  2:835761/1‑149 (MQ=255)
                                                                                                           gcTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAACGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGc                                       >  1:1029567/1‑149 (MQ=255)
                                                                                                           gcTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGAGGTCGTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGc                                       >  1:624542/1‑149 (MQ=255)
                                                                                                               aTGGGGCTGGCAAAATTATAACTCATGCCCCTGGAGTGGATTGAGCTTCCGGCGCTGGTGACGGTATTACCACATCAGCTAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAAtt                                   <  2:1424933/149‑1 (MQ=255)
                                                                                                                 ggggCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACTTCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTg                                 >  1:751138/1‑149 (MQ=255)
                                                                                                                  gggCTGGCAAAATTATTACTCAAGGCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGTCGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGc                                <  2:739729/149‑1 (MQ=255)
                                                                                                                       ggCAAAATTATTACTCAAGCCCCTGGAGGGGAGTGAGCGTCCGGCGCGGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTa                           <  1:810064/149‑1 (MQ=255)
                                                                                                                        gCAAAATTATTACTCAAGCCGCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGGCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTAt                          <  1:1267593/149‑1 (MQ=255)
                                                                                                                                attaCTCAAGCCCCTGGAGCGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATGGTGTCAGCACAATGCTGGTGCTAGTTGCCATGATGGGCTCAg                  >  1:1083504/1‑149 (MQ=255)
                                                                                                                                  taCTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCGGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGtt                <  2:1107176/149‑1 (MQ=255)
                                                                                                                                    cTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGCCTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTCt              >  2:1224847/1‑149 (MQ=255)
                                                                                                                                    cTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCAGTATGGGCTCAGTTCt              >  2:407631/1‑149 (MQ=255)
                                                                                                                                    cTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGAATGTGACACCACAATGCTGGTGCAATTTTCCATTATGGGCTCAGTTCt              >  2:634462/1‑149 (MQ=255)
                                                                                                                                     tCAATCCCCTCCAGGGGCTTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCTCAATGCTGGTGCAAGTTGCCATTATGGGCTCAGTTCTg             >  1:457344/1‑149 (MQ=255)
                                                                                                                                     tCAATCCCCTCCAGGGGCTTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTCTg             >  2:1557410/1‑149 (MQ=255)
                                                                                                                                     tCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTACTg             >  1:746441/1‑149 (MQ=255)
                                                                                                                                      aaaGCCCCTGGAGGGGATTGAGCGGCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTCTgg            <  1:808297/148‑1 (MQ=255)
                                                                                                                                         gCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTCTGGCTg         >  2:263090/1‑149 (MQ=255)
                                                                                                                                           cccTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTCTGGCTGaa       <  2:404798/149‑1 (MQ=255)
                                                                                                                                              gggAGGGTATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGGCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTCTGGCTGAagag    <  1:745357/148‑1 (MQ=255)
                                                                                                                                                gAGGGGATTGAGCGTCCGGCGCTGGTGCCGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGGTGGTGCAATTTGCCATTATGGGCTCAGTTCTGGCTGAAGAcgt  >  1:939607/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
CCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTCTGGCTGAAGAGGT  >  NC_000913/1147843‑1148135

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