Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,162,152 G→C 6.1% intergenic (+22/+18) yijD → / ← trmA conserved inner membrane protein YijD/tRNA m(5)U54 methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,162,1520GC 6.1% 383.8 / 12.1 114intergenic (+22/+18)yijD/trmAconserved inner membrane protein YijD/tRNA m(5)U54 methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (51/56);  new base C (4/3);  total (55/59)
Fisher's exact test for biased strand distribution p-value = 7.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.24e-01

Only 100 of 124 total aligned reads displayed.

CAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTCCAGATTCTTGCATAACGTTTCCGGGTTAC  >  NC_000913/4162004‑4162300
                                                                                                                                                    |                                                                                                                                                    
cAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGgc                                                                                                                                                     >  1:490718/1‑150 (MQ=255)
cAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGgc                                                                                                                                                     >  2:664740/1‑150 (MQ=255)
 aGCACTGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTggtg                                                                                                                                                    <  2:1045521/150‑3 (MQ=255)
 aGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGgcg                                                                                                                                                    <  1:1420184/150‑1 (MQ=255)
 aGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGgcg                                                                                                                                                    <  2:1395684/150‑1 (MQ=255)
  gCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGgcgc                                                                                                                                                   >  2:1238391/1‑150 (MQ=255)
  acACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGAAGTAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGgcgc                                                                                                                                                   >  1:204443/2‑150 (MQ=255)
   cACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCa                                                                                                                                                  >  1:612440/1‑150 (MQ=255)
          gTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCa                                                                                                                                           >  1:250257/1‑150 (MQ=255)
          gTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCGCCAGCATCCGGCa                                                                                                                                           <  1:1211678/150‑1 (MQ=255)
             cGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATgg                                                                                                                                        >  1:749366/1‑150 (MQ=255)
             cGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATgg                                                                                                                                        >  2:95978/1‑150 (MQ=255)
              gTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATggg                                                                                                                                       >  1:733929/1‑150 (MQ=255)
              gTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATggg                                                                                                                                       >  2:640916/1‑150 (MQ=255)
                  ggAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGtttt                                                                                                                                   <  1:1503232/150‑1 (MQ=255)
                       atCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTc                                                                                                                              <  2:490718/150‑1 (MQ=255)
                        tCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTACGCAGGAAGTGGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTcg                                                                                                                             <  2:1332530/150‑1 (MQ=255)
                          cGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTcgcg                                                                                                                           <  2:264992/150‑1 (MQ=255)
                            gATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGt                                                                                                                         <  1:95978/150‑1 (MQ=255)
                                tCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGt                                                                                                                     >  2:926033/1‑150 (MQ=255)
                                   gCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAAt                                                                                                                  >  2:1315535/1‑150 (MQ=255)
                                   gCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAAt                                                                                                                  >  2:14376/1‑150 (MQ=255)
                                   gCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAAt                                                                                                                  >  1:1287041/1‑150 (MQ=255)
                                      ctAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATAcg                                                                                                               >  1:696396/1‑150 (MQ=255)
                                       tAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATAcgc                                                                                                              <  2:912605/150‑1 (MQ=255)
                                       tAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATAcgc                                                                                                              >  2:89755/1‑150 (MQ=255)
                                       tAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATAcgc                                                                                                              >  2:1289260/1‑150 (MQ=255)
                                        aaCTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCgg                                                                                                                                             <  1:626129/118‑1 (MQ=255)
                                        aaCTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCgg                                                                                                                                             >  2:626129/1‑118 (MQ=255)
                                        aaCTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATAcgcc                                                                                                             <  2:103481/150‑1 (MQ=255)
                                          cttcttCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATAcgccgc                                                                                                           <  1:1562829/150‑1 (MQ=255)
                                          cttcttCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATAcgccgc                                                                                                           <  2:1447606/150‑1 (MQ=255)
                                          cttcttCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATCATGCCGGATGCGCCAGCATCCGGCAtaat                                                                                                                                       >  2:1322999/1‑119 (MQ=39)
                                          cttcttCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATCATGCCGGATGCGCCAGCATCCGGCAtaat                                                                                                                                       <  1:1322999/122‑4 (MQ=39)
                                             cttcCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACt                                                                                                        >  1:242800/1‑150 (MQ=255)
                                              ttcCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTc                                                                                                       >  1:639683/1‑150 (MQ=255)
                                                  cGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCata                                                                                                   <  2:1213642/150‑1 (MQ=255)
                                                      aGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGt                                                                                               >  2:319941/1‑150 (MQ=255)
                                                           cctcCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCgtg                                                                                          <  1:1112574/149‑1 (MQ=255)
                                                            ctcCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCgtgt                                                                                         >  1:777793/1‑150 (MQ=255)
                                                              ccGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAg                                                                                       >  1:380269/1‑150 (MQ=255)
                                                               cGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGAATGGGTTTTACTTCGCGGTCATTAATATTCCGCACTCCATATGGTGCGTGTAgt                                                                                      <  2:390849/150‑2 (MQ=255)
                                                                 tCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAgggg                                                                                    <  2:1610254/150‑1 (MQ=255)
                                                                    ttATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAAc                                                                                 <  2:1207972/150‑1 (MQ=255)
                                                                     tATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATTGTGCTTGTAGTGGAACt                                                                                <  2:722818/150‑1 (MQ=255)
                                                                      ctGGTGTGCTGGCCTGGCGCGAAGATGCGTAACCGAACGCAGGCAGGTGCTGAGCAATCGGTATTCTGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCCGTCCGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTg                                                                               <  2:473094/149‑1 (MQ=255)
                                                                         gtgtTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATc                                                                            <  1:1238391/150‑1 (MQ=255)
                                                                          tgtTCTGGATTGGCGCGAAGATGCTTAACCGTAAGCAGGAAGTTGCTGAGTAATGGGTATTATGCCGGATGCTGGCGCTTGCGGCATGTGTTTTACTTCGCTGTCAGTAATAGGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCa                                                                           <  2:204443/150‑1 (MQ=255)
                                                                           gtTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCaa                                                                          <  1:1405760/150‑1 (MQ=255)
                                                                             tCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAAc                                                                        <  1:133141/150‑1 (MQ=255)
                                                                             tCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAAc                                                                        >  1:1534761/1‑150 (MQ=255)
                                                                             tCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAAc                                                                        <  1:297694/150‑1 (MQ=255)
                                                                              cTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACa                                                                       >  2:43446/1‑150 (MQ=255)
                                                                                ggATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACaga                                                                     >  2:1139731/1‑150 (MQ=255)
                                                                                ggATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACaga                                                                     >  2:1515340/1‑150 (MQ=255)
                                                                                ggATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACaga                                                                     >  1:22062/1‑150 (MQ=255)
                                                                                   ttGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGcc                                                                  <  1:1530828/150‑1 (MQ=255)
                                                                                      gcgcGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGa                                                               >  2:1626850/1‑150 (MQ=255)
                                                                                        gcgAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGTGGAACTGATCAAACAGAGCCAGACg                                                             <  1:1264727/150‑1 (MQ=255)
                                                                                             gATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGa                                                        <  1:664740/150‑1 (MQ=255)
                                                                                              aTGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGAc                                                       <  2:22062/150‑1 (MQ=255)
                                                                                               tGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGAGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGAcc                                                      >  2:1069082/1‑150 (MQ=255)
                                                                                                 cGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCtt                                                    >  2:194086/1‑150 (MQ=255)
                                                                                                    aaCCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTtgtg                                                 <  1:372999/150‑1 (MQ=255)
                                                                                                       cGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGt                                              >  1:1574805/1‑150 (MQ=255)
                                                                                                        gTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTc                                             <  1:616913/150‑1 (MQ=255)
                                                                                                          aaGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTg                                           <  2:611242/150‑1 (MQ=255)
                                                                                                           aGCAGGAAGTGGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTgg                                          <  1:664738/150‑1 (MQ=255)
                                                                                                               ggAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTa                                      >  2:515063/1‑150 (MQ=255)
                                                                                                                    ttGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGtt                                 >  2:621439/1‑150 (MQ=255)
                                                                                                                    ttGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGtt                                 >  2:297269/1‑150 (MQ=255)
                                                                                                                     tGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGttt                                >  2:1294119/1‑150 (MQ=255)
                                                                                                                      gCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTc                               >  1:182841/1‑150 (MQ=255)
                                                                                                                      gCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTc                               >  1:1206104/1‑150 (MQ=255)
                                                                                                                       cTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTcc                              <  2:1426118/150‑1 (MQ=255)
                                                                                                                        tGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGGCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTCCa                             <  1:193990/150‑1 (MQ=255)
                                                                                                                         gAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTCCAg                            >  1:104166/1‑150 (MQ=255)
                                                                                                                          aGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTCCAGa                           <  1:620943/150‑1 (MQ=255)
                                                                                                                           gTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTCCAGAt                          <  1:1289260/150‑1 (MQ=255)
                                                                                                                            tAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTCCAGAtt                         >  1:665498/1‑150 (MQ=255)
                                                                                                                             aaTCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTCCAGATTc                        <  2:243002/150‑1 (MQ=255)
                                                                                                                              aTCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTCCAGATTCt                       <  2:733929/150‑1 (MQ=255)
                                                                                                                              aTCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTCCAGATTCt                       <  1:1159654/150‑1 (MQ=255)
                                                                                                                               tCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTCCAGATTCtt                      >  2:671739/1‑150 (MQ=255)
                                                                                                                                cGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTCCAGATTCTTg                     <  2:639574/150‑1 (MQ=255)
                                                                                                                                 ggTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTCCAGATTCTTGc                    >  1:1089845/1‑150 (MQ=255)
                                                                                                                                 aaacccatGCCGGATGCGCCAGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTCCAGATTCTTGc                    <  1:402406/144‑1 (MQ=255)
                                                                                                                                    attatGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTCCAGATTCTTGCATa                 <  2:612440/150‑1 (MQ=255)
                                                                                                                                    attatGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTCCAGATTCTTGCATa                 >  2:1085129/1‑150 (MQ=255)
                                                                                                                                      tatGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTCCAGATTCTTGCATAAc               <  1:14376/150‑1 (MQ=255)
                                                                                                                                      tatGCCGGATGCTGGCGCATCCGGCATCGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTCCAGATTCTTGCATAAc               <  1:1596722/150‑1 (MQ=255)
                                                                                                                                       atGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTCCAGATTCTTGCATAACg              >  2:1353900/1‑150 (MQ=255)
                                                                                                                                        tGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTCCAGATTCTTGCATAACGt             <  2:242476/150‑1 (MQ=255)
                                                                                                                                         gCCGGATGCGCCAGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTCCAGATTCTTGCATAACGtt            >  1:572312/1‑150 (MQ=255)
                                                                                                                                         gCCGGATGCGCCAGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTCCAGATTCTTGCATAACGtt            >  1:746460/1‑150 (MQ=255)
                                                                                                                                          ccGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTCCAGATTCTTGCATAACGttt           >  2:1357667/1‑150 (MQ=255)
                                                                                                                                           cGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTCCAGATTCTTGCATAACGTTTc          <  1:621439/150‑1 (MQ=255)
                                                                                                                                             gATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTCCAGATTCTTGCATAACGTTTCCg        <  1:643196/150‑1 (MQ=255)
                                                                                                                                             gATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTCCAGATTCTTGCATAACGTTTCCg        <  1:297269/150‑1 (MQ=255)
                                                                                                                                                   ggCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTCCAGATTCTTGCATAACGTTTCCGGGTTAc  <  1:1139731/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTCCAGATTCTTGCATAACGTTTCCGGGTTAC  >  NC_000913/4162004‑4162300

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