Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,162,153 C→A 5.9% intergenic (+23/+17) 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,1530CA 5.9% 395.8 / 7.9 116intergenic (+23/+17)yijD/trmAconserved inner membrane protein YijD/tRNA m(5)U54 methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (52/57);  new base A (4/3);  total (56/60)
Fisher's exact test for biased strand distribution p-value = 7.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 125 total aligned reads displayed.

CAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTCCAGATTCTTGCATAACGTTTCCGGGTTACAG  >  NC_000913/4162004‑4162302
                                                                                                                                                     |                                                                                                                                                     
cAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGgc                                                                                                                                                       >  2:664740/1‑150 (MQ=255)
cAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGgc                                                                                                                                                       >  1:490718/1‑150 (MQ=255)
 aGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGgcg                                                                                                                                                      <  2:1395684/150‑1 (MQ=255)
  gCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGgcgc                                                                                                                                                     >  1:219886/1‑150 (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)
             cGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATgg                                                                                                                                          >  2:95978/1‑150 (MQ=255)
             cGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATgg                                                                                                                                          >  1:749366/1‑150 (MQ=255)
              gTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATggg                                                                                                                                         >  2:640916/1‑150 (MQ=255)
              gTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATggg                                                                                                                                         >  2:1264727/1‑150 (MQ=255)
                  ggAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGtttt                                                                                                                                     <  1:1503232/150‑1 (MQ=255)
                       atCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTc                                                                                                                                <  2:490718/150‑1 (MQ=255)
                          cGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTcgcg                                                                                                                             <  2:264992/150‑1 (MQ=255)
                          cGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTcgcg                                                                                                                             <  1:845962/150‑1 (MQ=255)
                           ggATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCgg                                                                                                                            <  1:1340879/150‑1 (MQ=255)
                                tCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGt                                                                                                                       >  2:926033/1‑150 (MQ=255)
                                   gCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAAt                                                                                                                    >  1:1287041/1‑150 (MQ=255)
                                   gCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAAt                                                                                                                    >  2:14376/1‑150 (MQ=255)
                                   gCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAAt                                                                                                                    >  2:1315535/1‑150 (MQ=255)
                                      ctAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATAcg                                                                                                                 >  1:696396/1‑150 (MQ=255)
                                       tAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATAcgc                                                                                                                >  2:89755/1‑150 (MQ=255)
                                       tAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATAcgc                                                                                                                <  2:912605/150‑1 (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)
                                              ttcCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTc                                                                                                         >  1:639683/1‑150 (MQ=255)
                                               tcCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTcc                                                                                                        <  2:1399224/150‑1 (MQ=255)
                                                  cGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCata                                                                                                     <  2:1213642/150‑1 (MQ=255)
                                                      aGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGt                                                                                                 >  2:319941/1‑150 (MQ=255)
                                                       gTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTg                                                                                                <  1:1281070/150‑1 (MQ=255)
                                                         aCTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCg                                                                                              <  1:523218/150‑1 (MQ=255)
                                                           tctcCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCgtg                                                                                            <  1:1529295/150‑1 (MQ=255)
                                                            ctcCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCgtgt                                                                                           >  1:777793/1‑150 (MQ=255)
                                                              ccGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAg                                                                                         >  1:380269/1‑150 (MQ=255)
                                                               cGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGAATGGGTTTTACTTCGCGGTCATTAATATTCCGCACTCCATATGGTGCGTGTAgt                                                                                        <  2:390849/150‑2 (MQ=255)
                                                                gTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAggg                                                                                       >  2:1328390/1‑150 (MQ=255)
                                                                 tCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAgggg                                                                                      <  1:1190670/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)
                                                                            ttCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCaaa                                                                           <  2:219886/150‑1 (MQ=255)
                                                                             tCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAAc                                                                          <  1:297694/150‑1 (MQ=255)
                                                                             tCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAAc                                                                          <  1:133141/150‑1 (MQ=255)
                                                                              cTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACa                                                                         >  1:621232/1‑150 (MQ=255)
                                                                              cTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACa                                                                         >  2:43446/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)
                                                                                    tGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCa                                                                   >  1:69752/1‑150 (MQ=255)
                                                                                     ggCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAg                                                                  >  1:65508/1‑150 (MQ=255)
                                                                                      gcgcGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGa                                                                 >  2:1626850/1‑150 (MQ=255)
                                                                                        gcgAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGTGGAACTGATCAAACAGAGCCAGACg                                                               <  1:1264727/150‑1 (MQ=255)
                                                                                             gATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGa                                                          <  1:664740/150‑1 (MQ=255)
                                                                                               tGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGAGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGAcc                                                        >  2:1069082/1‑150 (MQ=255)
                                                                                                 cGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCtt                                                      >  2:194086/1‑150 (MQ=255)
                                                                                                       cGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGt                                                >  1:1574805/1‑150 (MQ=255)
                                                                                                        gTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTc                                               <  1:616913/150‑1 (MQ=255)
                                                                                                          aaGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTg                                             <  2:611242/150‑1 (MQ=255)
                                                                                                               ggAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTa                                        >  2:515063/1‑150 (MQ=255)
                                                                                                                    ttGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGtt                                   >  2:297269/1‑150 (MQ=255)
                                                                                                                    ttGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGtt                                   >  2:621439/1‑150 (MQ=255)
                                                                                                                     tGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGttt                                  >  2:1294119/1‑150 (MQ=255)
                                                                                                                      gCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTc                                 >  1:1206104/1‑150 (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:319941/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:1590518/1‑150 (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:1085129/1‑150 (MQ=255)
                                                                                                                                    attatGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTCCAGATTCTTGCATa                   <  2:612440/150‑1 (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:746460/1‑150 (MQ=255)
                                                                                                                                          ccGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTCCAGATTCTTGCATAACGttt             >  2:1357667/1‑150 (MQ=255)
                                                                                                                                            ggATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTCCAGATTCTTGCATAACGTTTcc           >  2:1113466/1‑150 (MQ=255)
                                                                                                                                             gATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTCCAGATTCTTGCATAACGTTTCCg          <  1:297269/150‑1 (MQ=255)
                                                                                                                                             gATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTCCAGATTCTTGCATAACGTTTCCg          <  1:643196/150‑1 (MQ=255)
                                                                                                                                                   ggCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTCCAGATTCTTGCATAACGTTTCCGGGTTAc    <  1:1139731/150‑1 (MQ=255)
                                                                                                                                                   ggCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTCCAGATTCTTGCATAACGTTTCCGGGTTAc    >  2:1557334/1‑150 (MQ=255)
                                                                                                                                                     cgcATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTCCAGATTCTTGCATAACGTTTCCGGGTTACAg  >  1:1088621/1‑150 (MQ=255)
                                                                                                                                                     cgcATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTCCAGATTCTTGCATAACGTTTCCGGGTTACAg  <  1:1069082/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTCCAGATTCTTGCATAACGTTTCCGGGTTACAG  >  NC_000913/4162004‑4162302

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