Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,329,708 A→G 5.4% G30G (GGT→GGC) ‡ yfaQ ← tandem DUF2300 domain protein, putative host defense protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,329,7080AG 5.4% 217.9 / 2.5 75G30G (GGT→GGC) ‡yfaQtandem DUF2300 domain protein, putative host defense protein
Reads supporting (aligned to +/- strand):  ref base A (32/38);  new base G (2/2);  total (34/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.82e-01

GGGCTGCTCCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCGCAACACTTTGCTGCGCAG  >  NC_000913/2329561‑2329853
                                                                                                                                                   |                                                                                                                                                 
gggCTGCTCCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCAcc                                                                                                                                                 <  2:1177944/150‑1 (MQ=255)
   ctgctCCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACcgcg                                                                                                                                              >  1:1480427/1‑150 (MQ=255)
    tgctCCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCgcgc                                                                                                                                             <  1:795026/150‑1 (MQ=255)
      ctCCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCgcag                                                                                                                                           <  1:1240156/150‑1 (MQ=255)
             tGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGCCTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGt                                                                                                                                    >  1:281607/1‑150 (MQ=255)
             tGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCGGTGCGCAGCAcgatc                                                                                                                                    <  2:453284/150‑5 (MQ=255)
              gATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTAGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGtt                                                                                                                                   >  1:211837/1‑150 (MQ=255)
               aTGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTg                                                                                                                                  <  2:191726/150‑1 (MQ=255)
                tGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGt                                                                                                                                 <  1:960513/150‑1 (MQ=255)
                  ggTATCCTCCAGCCTGGCGTAAACATACAGTTTCCACAAACGCCCTATCGGCGTGGTGAGCGAGTCAGGTAAGGCTGTGACTTTGGGCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTaa                                                                                                                               <  2:1344288/150‑1 (MQ=255)
                        ctccAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGt                                                                                                                         >  1:1630585/1‑150 (MQ=255)
                          ccagccagGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTTGTCACGCCTGATGACGACAGCTGATACAGTTTATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGtct                                                                                                                       <  1:1516931/150‑1 (MQ=255)
                             gccagGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTctcct                                                                                                                    >  1:940656/1‑150 (MQ=255)
                                   cGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCaa                                                                                                              >  2:420692/1‑150 (MQ=255)
                                    gTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAg                                                                                                             <  1:1732170/150‑1 (MQ=255)
                                         cATACAGTTTCCACAAACTCCCTAAGGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTg                                                                                                        >  2:1339142/1‑150 (MQ=255)
                                         cATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTg                                                                                                        >  2:1587382/1‑150 (MQ=255)
                                              aGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCaga                                                                                                   >  2:227617/1‑150 (MQ=255)
                                              aGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCaga                                                                                                   >  2:166681/1‑150 (MQ=255)
                                                 ttCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTa                                                                                                >  2:1508560/1‑150 (MQ=255)
                                                 ttCCACAAACTCCCTAACGGCGGGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTa                                                                                                <  1:1045424/150‑1 (MQ=255)
                                                   ccACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAAc                                                                                              >  1:55369/1‑150 (MQ=255)
                                                        aaCTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACg                                                                                         >  1:887821/1‑150 (MQ=255)
                                                                  cGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCACCGCCAATAGCCa                                                                               >  2:330110/1‑150 (MQ=255)
                                                                  cGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCa                                                                               >  2:1682095/1‑150 (MQ=255)
                                                                   ggCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGACAGAGTAACCAACGCCAAAAGCCac                                                                              >  2:1533722/1‑149 (MQ=255)
                                                                    gCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGa                                                                             >  2:1311345/1‑150 (MQ=255)
                                                                     cGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGGTGATCGTGCTGCGCACCGCGCAGCACCAGTTGAAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGAc                                                                            <  1:897800/150‑1 (MQ=255)
                                                                        ggtgAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGAt                                                                         >  2:506903/1‑150 (MQ=255)
                                                                        ggtgAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGAGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGAt                                                                         >  2:118285/1‑150 (MQ=255)
                                                                            aGCGAGTCAGTTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTg                                                                     <  2:656325/150‑1 (MQ=255)
                                                                                aGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGcc                                                                 >  1:863822/1‑150 (MQ=255)
                                                                                 gTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCa                                                                >  2:724536/1‑150 (MQ=255)
                                                                                  tCAGGTAATGCTGAGACTTTGGGCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAg                                                               <  1:1196684/150‑1 (MQ=255)
                                                                                     ggaaaTGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGATGTAACGCCGTCTCCGCCGCAAGAGTTGGCACAGTAACCAACACCAAAAGCCCGACGATTCTGCGCCAGTTc                                                            >  2:1447788/4‑150 (MQ=255)
                                                                                       tAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCAc                                                          <  2:968077/150‑1 (MQ=255)
                                                                                        aaTGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACt                                                         <  1:571992/150‑1 (MQ=255)
                                                                                        aaTGCTGAGACTTTGGTCACGCCGGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGGGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACt                                                         <  2:267021/150‑1 (MQ=255)
                                                                                           gCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTAt                                                      <  1:1145402/150‑1 (MQ=255)
                                                                                           gCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTAt                                                      >  1:921761/1‑150 (MQ=255)
                                                                                             tGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTAttt                                                    <  1:1374416/150‑1 (MQ=255)
                                                                                              gagaCTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTc                                                   >  2:1704152/1‑150 (MQ=255)
                                                                                              gagaCTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTc                                                   <  1:1205348/150‑1 (MQ=255)
                                                                                                 aCTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCAcc                                                <  2:1664694/150‑1 (MQ=255)
                                                                                                 aCTTTGGGCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGGCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCAcc                                                <  2:1161446/150‑1 (MQ=255)
                                                                                                   tttGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTg                                              <  2:211837/150‑1 (MQ=255)
                                                                                                    ttGGGCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGc                                             <  2:990648/150‑1 (MQ=255)
                                                                                                        tCACGCCTGATGACGACAGCTGATACAGGTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATcc                                         <  2:916811/150‑1 (MQ=255)
                                                                                                         cACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGACAGTTCACTTATTTCACCTGCATcca                                        >  2:126335/1‑149 (MQ=255)
                                                                                                            gCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGt                                     <  2:584500/150‑1 (MQ=255)
                                                                                                              cTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTcc                                   <  1:420692/150‑1 (MQ=255)
                                                                                                                 ttgaCGACAGCTGATCCTGTTTATCGTGCTGCGCTCCGCGCAGCACCAGTTGTATCGGCGTCTCCTCCGCAAGTGTTGGCAGAGTTACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATATCACCTGCATCCCGGTCCAtt                                <  1:1504395/149‑1 (MQ=255)
                                                                                                                  tgaCGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTc                               <  1:227617/150‑1 (MQ=255)
                                                                                                                  tgaCGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTc                               <  1:166681/150‑1 (MQ=255)
                                                                                                                    acgacAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCAc                             >  1:1200910/1‑150 (MQ=255)
                                                                                                                    acgacAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCAc                             >  1:1415552/1‑150 (MQ=255)
                                                                                                                     cgacAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACt                            <  1:372345/150‑1 (MQ=255)
                                                                                                                     caacAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCCCCTGCATCCCGGTCCATTCACt                            >  2:673638/3‑150 (MQ=255)
                                                                                                                       acAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTcc                          <  1:1682095/150‑1 (MQ=255)
                                                                                                                          gCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACAAACGCCAATAGCCAGACGATCCTGCGCCAGTTCACTTATTTCACCTGCCTCCCGGTCCATTCCCTCCCTg                       >  1:125191/1‑150 (MQ=255)
                                                                                                                            tGATACAGTTGATCGTGCTGCGCACCGCTCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGcc                     <  1:1508560/150‑1 (MQ=255)
                                                                                                                            tGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGcc                     <  1:330110/150‑1 (MQ=255)
                                                                                                                               tACAGTTGCTCGGGCTGCGCACCGCGCAGCACCAGTTGTACCGGCGGCTCCTCCGCAAGTGTTGGCAGAGTAACCCACGCCAATAGCCCGACGATTCTGCGCCCGTTCACTTATTTCACCTGCATCCGGGTCCATTCACTCCCTGCCGCa                  <  2:948532/150‑1 (MQ=255)
                                                                                                                               tACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGGGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGGCCATTCACTCCCTGCCGCa                  <  1:1374117/150‑1 (MQ=255)
                                                                                                                                 cAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCGCAac                >  2:269483/1‑150 (MQ=255)
                                                                                                                                    ttGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCGCAACACt             <  1:1122594/150‑1 (MQ=255)
                                                                                                                                      gATCGTGCTGCGCACCGCGCAGCACCAGTTGTATCGGCGTCTCATCCGAATGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCGCAACACttt           <  2:863822/150‑1 (MQ=255)
                                                                                                                                       ctggTGCTGCGCGGTGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCGCAACACtt            <  1:719956/145‑1 (MQ=255)
                                                                                                                                       ctggTGCTGCGCGGTGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCGCAACACtt            >  2:719956/4‑148 (MQ=255)
                                                                                                                                          gTGCTGCGCGGTGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCGCAACACTTtgctg       >  1:1212786/1‑150 (MQ=255)
                                                                                                                                          gTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATGCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCGCAACACTTtgctg       <  1:1117202/150‑1 (MQ=255)
                                                                                                                                            gctgcGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCTCTCCCTGCCGCAACACTTTGCTgcg     >  2:631690/1‑150 (MQ=255)
                                                                                                                                            gctgcGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCTCTCCCTGCCGCAACACTTTGCTgcg     >  2:1206984/1‑150 (MQ=255)
                                                                                                                                             ctgcGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCGCAACACTTTGCTgcgc    <  1:724536/150‑1 (MQ=255)
                                                                                                                                               gcgcACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCGCAACACTTTGCTGCGCAg  <  1:444537/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GGGCTGCTCCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCGCAACACTTTGCTGCGCAG  >  NC_000913/2329561‑2329853

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: