Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,074,844 G→A 8.5% G59D (GGT→GAT)  yihW → putative transcriptional regulator for sulphoquinovose utilization

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,074,8440GA 8.5% 315.8 / 12.3 105G59D (GGT→GAT) yihWputative transcriptional regulator for sulphoquinovose utilization
Reads supporting (aligned to +/- strand):  ref base G (42/54);  new base A (3/6);  total (45/60)
Fisher's exact test for biased strand distribution p-value = 7.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.43e-01

Only 100 of 107 total aligned reads displayed.

CCCGCGGCACGACCAACTCCTCATGCTGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTCAACAATATTTATCACCATTGGTACGACTGTTG  >  NC_000913/4074695‑4074993
                                                                                                                                                     |                                                                                                                                                     
cccGCGGCACGACCAACTCCTCATGCTGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCAGCACgg                                                                                                                                                       >  2:611721/1‑150 (MQ=255)
 ccGCGGCACGACCAACTCCTCATGCTGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGt                                                                                                                                                      >  2:1696323/1‑150 (MQ=255)
  cgcgGCACGACCAACTCCTCATGCTGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTg                                                                                                                                                     >  2:1443121/1‑150 (MQ=255)
  cgcgGCACGACCAACTCCTCATGCTGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTg                                                                                                                                                     <  1:973259/150‑1 (MQ=255)
   gcgGCACGACCAACTCCTCATGCTGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTgg                                                                                                                                                    <  1:1139525/150‑1 (MQ=255)
    cgGCACGACCAACTCCTCATGCTGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGgc                                                                                                                                                   <  1:1282348/150‑1 (MQ=255)
       cACGACCAACTCCTCATGCTGATCGCCGATCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGcgcg                                                                                                                                                <  2:763430/150‑1 (MQ=255)
        aCGACCAACTCCTCATGCTGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAACTTAAGCGAGCAAGGCCTGATTACGCGCCACCGTGATGGcgcgt                                                                                                                                               <  1:1088793/150‑2 (MQ=255)
        aCGACCAACTCCTCATGCTGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCACCGTGATGGcgcgt                                                                                                                                               <  1:787732/150‑2 (MQ=255)
         cGACCAACTCCTCATGCTGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCggg                                                                                                                                              <  1:1127946/150‑1 (MQ=255)
         cGACCAACTCCTCATGCTGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCACCGTGATGGcgcgta                                                                                                                                              <  2:1659138/150‑3 (MQ=255)
          gACCAACTCCTCATGCTGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGt                                                                                                                                             >  1:1342710/1‑150 (MQ=255)
          gACCAACTCCTCATGCTGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCACCGTGATGGcgcgtaa                                                                                                                                             <  1:1377094/150‑4 (MQ=255)
              aaCTCCTCATGCTGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCggg                                                                                                                                         <  1:607969/150‑1 (MQ=255)
              aaCTCCTCATGCTGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCggg                                                                                                                                         <  2:1452306/150‑1 (MQ=255)
              aaCTCCTCATGCTGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCACCGTGATGGcgcgtaatcag                                                                                                                                         <  1:1293588/150‑8 (MQ=255)
               aCTCCTCATGCTGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGc                                                                                                                                        <  1:1160180/150‑1 (MQ=255)
                  cctcATGCTGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTc                                                                                                                                     <  2:358766/150‑1 (MQ=255)
                   ctcATGCTGATCGCCGAGCGTGGGTATATGAATATTGCTGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTcc                                                                                                                                    <  1:14590/150‑1 (MQ=255)
                     cATGCTTATCGCCGAGCGGGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGACTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAg                                                                                                                                  <  1:117392/150‑1 (MQ=255)
                     cATGCTGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAg                                                                                                                                  <  1:495159/150‑1 (MQ=255)
                       tGCTGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGcg                                                                                                                                <  1:293691/150‑1 (MQ=255)
                       tGCTGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCACCGTGATGGcgcgtaatcag                                                                                                                                         <  2:171665/141‑8 (MQ=255)
                       tGCTGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCACCGTGATGGcgcgtaatcag                                                                                                                                         >  1:171665/1‑134 (MQ=255)
                         cTGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGcgtc                                                                                                                              >  1:1630471/1‑150 (MQ=255)
                          tGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGcgtcg                                                                                                                             >  2:359600/1‑150 (MQ=255)
                          tGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGcgtcg                                                                                                                             <  1:1161043/150‑1 (MQ=255)
                          tGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGcgtcg                                                                                                                             >  2:252671/1‑150 (MQ=255)
                               gCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAAt                                                                                                                        >  2:1013964/1‑150 (MQ=255)
                                 cGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATAc                                                                                                                      >  2:1531933/1‑150 (MQ=255)
                                  gAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACg                                                                                                                     <  2:762830/150‑1 (MQ=255)
                                   aGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCCCGCAGACGGACCGCCGGGATATTCGTAAATTAAGCGATCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGGCGGGCTTCCAGCGTCGTTATTACgg                                                                                                                    >  2:1467294/1‑150 (MQ=255)
                                       tGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGtt                                                                                                                >  2:208126/1‑150 (MQ=255)
                                       tGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGAAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGtt                                                                                                                >  1:1303855/1‑150 (MQ=255)
                                        gggTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTc                                                                                                               <  1:1013964/150‑1 (MQ=255)
                                          gTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGCCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGa                                                                                                             >  2:1396748/1‑150 (MQ=255)
                                           tataTGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGag                                                                                                            >  2:745902/1‑150 (MQ=255)
                                           taaatGAATATTGACGAGCTGGCAAAACTGCTGGATGGCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCCCGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGag                                                                                                            <  1:1578442/147‑1 (MQ=255)
                                             tatGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGagca                                                                                                          >  2:637897/1‑150 (MQ=255)
                                               tGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGagcagc                                                                                                        >  1:457126/1‑150 (MQ=255)
                                                gAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCg                                                                                                       <  1:745902/150‑1 (MQ=255)
                                                  atatTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGGTAATACGGCGTTCGAGCAGCGTg                                                                                                     >  2:153512/1‑150 (MQ=255)
                                                   tatTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCgaga                                                                                                    >  2:770864/1‑150 (MQ=255)
                                                   tatTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGa                                                                                                    >  2:1187110/1‑150 (MQ=255)
                                                    tttGATTAGCTGGCAAATATGCTGGAGGTCTCCAGGCAGACGGTCCTCCGGGATACTCGTAAATTAAGCGAGCAAGGCCTGTTTACGCTCCAGCAGTGTTGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAg                                                                                                   <  1:544455/149‑1 (MQ=255)
                                                      tgatgaGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCAGCACGGGGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGt                                                                                                 <  1:1581662/150‑1 (MQ=255)
                                                      tgatgaGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGCGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCCCGGTGGCGCGGGTCGGGCTTCCAGAGTCGTTAAAACGGCGTTCGAGCCGCGTGAGGt                                                                                                 >  1:298705/1‑150 (MQ=255)
                                                       gatgaGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGtt                                                                                                >  1:609743/1‑150 (MQ=255)
                                                         tgaGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGTAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTc                                                                                              >  1:836378/1‑150 (MQ=255)
                                                          gaGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCg                                                                                             <  2:69552/150‑1 (MQ=255)
                                                          gaGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCg                                                                                             >  2:297936/1‑150 (MQ=255)
                                                           cgCTGGCAACTGTGCTGGCTCTCTCCCCGCAGACGGTCCGCCGGGATATTCGTAACTTACGCGAGCAATGCCTGATTCCGCGCCATCACGGCGGCGCGTGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGc                                                                                            <  2:289583/149‑1 (MQ=255)
                                                           aGCTGGCAAATCTGCTGGATGGCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTACGCGAGCAAGGCCTGATTACGCGCCATCACGGGGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGc                                                                                            <  1:1292811/150‑1 (MQ=255)
                                                            gCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCa                                                                                           <  2:1472050/150‑1 (MQ=255)
                                                               ggCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCAATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAAc                                                                                        >  2:307333/1‑150 (MQ=255)
                                                                     tCTGCGGGATGGCGCCACGCAGACGCTCCTCCGGGATATTCGTAAAATCAGCGAGCCAGGCCTGATTACGCGCCATCTCGGTGGCGCGGGTCGGGCTTCCAGCGTCGTAAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGa                                                                                  <  1:431882/150‑1 (MQ=255)
                                                                       tgctgGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGaaa                                                                                <  1:917368/150‑1 (MQ=255)
                                                                             aTGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAGAAAAAAg                                                                          <  1:1696323/150‑1 (MQ=255)
                                                                                   ccACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTg                                                                    <  2:1342710/150‑1 (MQ=255)
                                                                                     aCGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGcc                                                                  <  1:1458310/150‑1 (MQ=255)
                                                                                       gCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGa                                                                <  1:208126/150‑1 (MQ=255)
                                                                                       gCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGa                                                                <  2:530097/150‑1 (MQ=255)
                                                                                       gCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGACCAAGGCCTGAGTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGGCGTTAATACGGCGTGCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGa                                                                <  2:134248/150‑1 (MQ=255)
                                                                                          gACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGc                                                             <  1:1443121/150‑1 (MQ=255)
                                                                                             ggTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGt                                                          <  1:297936/150‑1 (MQ=255)
                                                                                               tCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTgg                                                        <  1:1531933/150‑1 (MQ=255)
                                                                                                 cgccgGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCa                                                      <  2:1630471/150‑1 (MQ=255)
                                                                                                    ctggATTTTCTTAAATTAAGCGATCAAGGCCTGTTTACGCGCCATCACTGTGGCGCGGGTCGGGCTTCCAGCGTGGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACGGAGGAATAAAAATCGTTTGCCGAAGCGGTGGCAGAc                                                   <  2:694515/148‑1 (MQ=255)
                                                                                                       gATATTCGTAAATAAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACtat                                                <  2:867962/150‑1 (MQ=255)
                                                                                                          atTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGGCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATAtt                                             <  2:363007/150‑1 (MQ=255)
                                                                                                              gTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTg                                         <  2:425374/150‑1 (MQ=255)
                                                                                                               tAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGa                                        >  1:1582219/1‑150 (MQ=255)
                                                                                                               tAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGa                                        >  1:274686/1‑150 (MQ=255)
                                                                                                                aaaTTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGAt                                       >  1:607138/1‑150 (MQ=255)
                                                                                                                 aaTTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATg                                      >  2:1211583/1‑150 (MQ=255)
                                                                                                                 aaTTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATg                                      >  2:1205585/1‑150 (MQ=255)
                                                                                                                      tgCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTca                                 <  1:1638065/149‑1 (MQ=255)
                                                                                                                        cGAGCAAGGCCTGATTACGCGCCATCAGGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTcaac                               <  2:1640708/150‑1 (MQ=255)
                                                                                                                        cGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTcaac                               <  1:829732/150‑1 (MQ=255)
                                                                                                                        cGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTcaac                               >  1:525212/1‑150 (MQ=255)
                                                                                                                         gAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTcaaca                              >  1:701742/1‑150 (MQ=255)
                                                                                                                          aGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTcaacaa                             >  2:195958/1‑150 (MQ=255)
                                                                                                                            cAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTCAACAata                           <  1:252671/150‑1 (MQ=255)
                                                                                                                            cAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTCAACAata                           <  1:359600/150‑1 (MQ=255)
                                                                                                                               ggCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCACACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTCAACAATAttt                        >  2:162225/1‑150 (MQ=255)
                                                                                                                               ggCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTCAACAATAttt                        >  2:941250/1‑150 (MQ=255)
                                                                                                                                cccTGATTACGCGCCATCACGGGGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTCAACAATATTTa                       <  1:1533384/149‑1 (MQ=255)
                                                                                                                                 ccTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTCAACAATATTTAt                      >  1:1303390/1‑150 (MQ=255)
                                                                                                                                  cTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTCAACAATATTTATc                     >  2:587297/1‑150 (MQ=255)
                                                                                                                                  cTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTCAACAATATTTATc                     >  2:760177/1‑150 (MQ=255)
                                                                                                                                   tGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTCAACAATATTTATCa                    >  2:483368/1‑150 (MQ=255)
                                                                                                                                   ccgacccgcgCCACCGTGATGGCGCGGGTCGGGCTTCCAGCGACGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTCAACAATATTTATCa                    >  1:1355332/7‑150 (MQ=255)
                                                                                                                                        aCGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTCAACAATATTTATCACCAtt               <  1:143186/150‑1 (MQ=255)
                                                                                                                                         cgcgCCACCGTGATGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTCAACAATATTTATCACCATTg              >  1:1672749/1‑150 (MQ=255)
                                                                                                                                          gcgcCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTCAACAATATTTATCACCATTgg             <  2:701742/150‑1 (MQ=255)
                                                                                                                                               aTCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCCATCCGAGGGAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTCAACAATATTTATCACCATTGGTACGa        >  1:279495/1‑150 (MQ=255)
                                                                                                                                                tCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTCAACAATATTTATCACCATTGGTACGAc       <  1:637897/150‑1 (MQ=255)
                                                                                                                                                   cGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTCAACAATATTTATCACCATTGGTACGACTGt    <  1:1689231/150‑1 (MQ=255)
                                                                                                                                                   cGGTGGCGCGGGTCGCGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTCAACAATATTTATCACCATTGGTACGACTGt    >  2:663438/1‑150 (MQ=255)
                                                                                                                                                     gTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTCAACAATATTTATCACCATTGGTACGACTGTTg  <  2:609743/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCCGCGGCACGACCAACTCCTCATGCTGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTCAACAATATTTATCACCATTGGTACGACTGTTG  >  NC_000913/4074695‑4074993

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